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Several formulae are derived expressing the intensity of electron scattering for a variety of 
orientations of hydrocarbon films on a solid surface. Intensity contour maps may be constructed 
from these formulae for arbitrarily chosen models. These maps are helpful for determining the 
structure of a film and are being prepared for future publication. 


ONSIDERABLE attention has been given 
to the study of the unusual properties of 
long chain hydrocarbon films on various solid 
surfaces. An example is the field of boundary 
lubrication where proper significance .has been 
given to the importance of the strong adhesive 
forces between an oil film and the metal surface 
being lubricated.! These films are at times perhaps 
only one molecule thick, oriented so that their 
long molecular axes are vertical or almost so. 
Electrons are ideal for studying such films be- 
cause of their low penetrating power. As a result 
the major portion of the scattering arises from 
the outermost layers of atoms when the beam 
strikes at a glancing angle. Murison,’ using fast 
electrons, performed the first experiments on 
long chain organic compounds. A study of 
lubricating oils was subsequently made using 
electron diffraction by Andrew* who obtained a 
* This work was performed in the Chemistry Department 
of the University of Michigan as part of a cooperative 
project between the University and the Naval Research 
Laboratory under Contract Number N173s—10452. 
_ | Fora general discussion containing many references one 
is referred to the book by N. K. Adam, The Physics and 
Chemistry of Surfaces (Oxford University Press, New York, 
1941). See also The Chemical Background for Engine Re- 
search (Interscience Publishers, Inc., New York, 1943), 
‘Lubrication,’ by O. Beeck, page 265. 


*C. A. Murison, Phil. Mag. 17, 201 (1934). 
§L. T. Andrew, Trans. Faraday Soc. 32, 607 (1936). 


correlation between the lubricating qualities of 
the oils and their corresponding patterns. Another 
interesting study of surface films has been made 
by Germer and Storks.‘ They investigated the 
structure of monomolecular and multimolecular 
films of barium stearate and stearic acid prepared 
by the Langmuir-Blodgett technique. 

Thus far an expression for the intensity of 
scattering has been given only for the case where 
the molecules are distributed randomly about 
their axes and these axes are normal to the 
surface of attachment.” The expressions to be 
derived below give the intensity of scattering for 
more general spatial arrangements and by means 
of these, more detailed information concerning 
the structure of surface films may be obtained. 


MATHEMATICAL DEVELOPMENT 


We wish to obtain the intensity of coherent 
scattering of fast electrons (ca. 30,000 volt 
electrons) by long hydrocarbon chains in arbi- 
trary orientations relative to the surface of 
attachment. It is assumed that a chain;is com- 
posed of two parallel rows of carbon atoms in a 
staggered configuration, and that the contribu- 
tion from hydrogen atoms and end groups, such as 


4L. H. Germer and K. H. Storks, J. Chem. Phys. 6, 280 
(1938). 
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Fic. 1. End of hydrocarbon chain rests on flat surface. 
The direct beam (mo) and the scattered beam (n) reach the 
photographic plate (x, y).° 


carboxyl, may be neglected. The distribution of 
the chains is considered to be sufficiently dilute so 
that each one scatters independently. 

In order to discuss the various spatial arrange- 
ments which will be considered, we define a set of 
parameters in terms of which the orientation of 
the carbon chain may be completely determined. 
We place the carbon atom closest to the sup- 
porting surface at the origin and define the row in 
which it occurs as the axis of the molecule, vector 
ra (see Fig. 1). The other row of atoms may be 
free to rotate about this axis and the angle of 
rotation, 6, is measured relative to a vertical 
plane containing axis rq. By fixing the position of 
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r, in space and the value of angle 6, the orienta- 
tion of the carbon chain is determined. 

Expressions for the intensity of scattering will 
be obtained for the following cases: 

I. The molecules are fixed in space. Thus, 
angles fq and ¢aq, the angles which determine the 
position of rz, and angle 6 which determines the 
orientation of the plane of the carbon chain 
assume specific arbitrary values.°® 

II. The axes of the molecules are arbitrarily 
fixed in space but the planes of the carbon chains 
are considered to have random orientations about 
these axes. °® 

III. The tilt of the molecular axis and the 
orientation of the plane of the chain are arbitrarily 
chosen. However, the molecules are randomly 
oriented about axis z. 

IV. The tilt of the molecular axis is arbitrarily 
fixed. However, the molecules are randomly 
oriented about the vertical, z, and also about 
their own axes. 

Our next step is to obtain an expression for the 
amplitude of electron scattering for a generalized 
molecule consisting of 7 atoms, bearing in mind 
that the intensity is equal to the square of the 
absolute magnitude of the amplitude. The general 
expression is then applied to the above listed 
cases for the staggered carbon chain. 

We describe by r; the vector from the origin to 
the atom 7 (see Fig. 2). As an electron beam 


z 





| Ne 





Fic. 2. Vectors defining path difference, n-r;—mMo- fi. 


5 Angle ya is the tilt of ra from the vertical axis, z, and 
angle ¢. is the clockwise rotation (when looking down on 
the surface) of the raz plane from the plane bounded by 
and the incident beam after traversing the origin, moz. 

6 The case of free rotation about an axis cannot be dis- 
tinguished from fixed random orientations among the 
many molecules. 
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strikes the molecule, each atom scatters a wavelet ; 
and the path difference between the wavelets 
scattered by the atom at the origin and the atom 
iis (M—Mo)-Fs, where n and ny are unit vectors in 
the direction of the scattered wave and the 
incident wave, respectively. If the phase of the 
wavelet scattered at the origin is considered to be 
zero, the amplitude of the wave scattered by the 
carbon chain in the direction n is proportional to 


s > exp [ik(n—no)-r;], (1) 


Y i=1 


where f, is the scattering factor for carbon, a 
function of n and ng, r is the radial coordinate of 
the scattered beam measured from the origin to 
any reference point in space (in our case points on 
the photographic plate), k=27/\ where X is the 
electron wave-length and N is the number of 
carbon atoms. The regularity of the structure of 
the hydrocarbon chain allows one to sum expres- 
sion (1) as a power series before squaring its 
absolute value to obtain the intensity of scat- 
tering. Cases I and II are then obtained in a 
compact form. However, for Cases III and IV 
there is no gain in summing the series and so 
instead, a general expression for the intensity is 
obtained from (1) by multiplying it by its com- 
plex conjugate. The result is 


ferN N 
— > > exp [tk(n—np) - 135], (2) 


y? i=l j=1 


where r;; is a vector connecting the 7th and jth 
atoms. Expression (2) is used for deriving 
formulae III and IV, resulting in less compact 
expressions than those of the first two cases. 


Case I 


The plan of this section is (a) to write expres- 
sion (1) in the form of a power series and sum it, 
(b) to obtain the path differences included in (a) 
in terms of the internal structural parameters of 
the chain ya, 7a, and r; (defined below), functions 
of the polar coordinates of the scattered beam, 
and the orientation parameters Wa, ¢a, and 6, and 


(c) to insert the results of (6) into (a) and arrive. 


at an expression for the intensity. Finally the 
relations between the polar coordinates of the 
scattered beam and the functions thereof are 








Fic. 3. Trigonometric details at point of scattering. 


given. We thus may express the intensity in 
terms of the internal structural parameters, the 
orientation parameters, and the coordinates of 
the scattered beam. 

(a) We may write expression (1), for a chain 
containing an even number of carbon atoms, JN, 
in the form of a power series as follows: 


fe N/2-1 
— > exp [ikg(n—no)-r; | 


Yr g=0 
X(1+exp [7k(n—mno)-ra]), (3) 


where rz is the vector which connects the carbon 
atom at the origin to the nearest one and rf; is the 
vector which connects the carbon atom at the 
origin to the nearest one in the same row. The 
sum of (3) is 


1 —exp [iNk(n—np) ‘r,/2] 





1—exp [ik(n—no) -r, ] 
X(1+exp [7k(m—mo)-ta]). (4) 


(6) The angle between n—mnp and r% is a and 
the angle between n—n» and fy is ag (see Fig. 3). 
We proceed to relate these angles to the structural 
and orientation parameters and functions of the 
coordinates of the scattered beam. Let us define 
the following quantities: 


r,=a vector perpendicular to r, in the plane of 
the carbon chain and placed at the origin. 

r.=a vector perpendicular to r, in the raz plane 
and placed at the origin, 
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The angle of rotation of the chain, 6, is clearly the 
angle between r, and rp. 


B=angle between n—n» and r.. 
7t=angle between n—np and r,. 
va=angle between r, and rq. 
p=angle between planes defined by r, and 
rz and rz, and n—np. 
¢=angle between planes defined by r,. and 
r, and r, and n—np». 
yY,=angle of tilt of n—no from the vertical 
axis, 2. 
¢,=angle of rotation of n—ny from the noz 
plane.’ 
y.=angle of tilt of r, from the vertical axis, z. 
¢.=angle of rotation of r, from the moz 
plane.’ 


From the sine and cosine laws of spherical 
trigonometry, the following relations are ob- 
tained: 


COS Gq =COS Gq COS Ya+SiN a SiN Ya COS p, 
COS T=SiN @q COS p, 
cos t=cos 6 cos B+sin 6 sin B sin ¢, 
sin {=sin a, sin (6—p)/sin B, 

COS Q=COs Wx, COS Wa 


+sin y, sin , cos (¢a—¢n), 


(5) 


and 
cos B=cosy, cos ¥.+sin y, sin yz cos (¢:—¢n), 
where 


g:=t+¢o. and y,=90—ya. 


By combining the first four expressions in (5) we 
get® 
cos ag=ko+k1 cos 6k: sin 4, (6) 
where 
kyo =COS a COS Ya, ki =sin ya cos B, 
and 


ko=sin ya sin vp, sin (da— on). 


(c) We multiply (4) by its complex conjugate, 
introduce the expressions for cos aq and cos aq 
obtained in (5) and (6), and find the intensity 


7 These angles are measured by the same convention as 
applies to angle ¢a. 

8 The ambiguity regarding the sign may be resolved by 
defining 6 as increasing when the rotation of 4 about the 
axis of the chain, fj, appears counterclockwise to one 
observing the tip, of r. and facing the origin. Then a 
— check shows that the negative sign must be 
taken. 
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proportional to® 


f.? sin? (Nsr; cos aq/4) 





r? sin? (sr; COS a/2) 


[1-++cos (sra.cos aa) |, (7) 





where s=k|n—no 

We denote the angular coordinates of the 
scattered beam, n, by y, the tilt from the vertical 
axis, z, and by ¢ the rotation from the noz plane 
using the same conventions as apply to angle du. 
The following readily obtained relations enable 
us to express (7) and the subsequent expressions 
in II, III, and IV, in terms of these coordinates. 


s=4n[1—cos ¢ sin y/2 }}/A, 
s cos ~,=22r cos y/X, 
s sin ¥,=22[2(1—cos ¢ sin yp) 
—cos?y]!/A, (8) 
ssin ~, cos $, = —2mr(1—cos ¢ sin W) /X, 
s sin y, sin ¢,= +22[2(1—cos (+¢) sin y) 
—cos? y—(1—cos (4+¢) sin wy)? ]3/A.1 


The scattering factor, f., may also be expressed in 
terms of the angular coordinates of the scattered 
beam since it is equal to 6— F./s* where F, is the 
x-ray scattering factor for carbon atoms." 

Cartesian coordinates may be defined for the 
photographic plate in terms of the main beam as 
the origin, an x axis parallel to the plane of the 
supporting surface and a y axis parallel to the 
surface normal, z. Since 

sin ¢=x/(L?+x?)} (9) 
and 
cos y=y/(L?+x?+y’)}, 


where L is the distance from the point of scatter- 
ing to the origin of the photographic plate, the 
relations in (8) may be expressed as accurately as 
desired in terms of x and y. The fact that L is 
much larger than x or y allows one to make 


® Note that Nis even. When JN is odd an expression much 
less compact than (7) is obtained. Because of the rarity 
of occurrence of odd chains, it does not seem worth while 
to discuss this case further at present. 

10 We are interested in the scattered intensities occurring 
in the solid angle above the undeviated beam included by 
a photographic plate. Thus, only small values of ¢ are 
involved. If we mean by +@ a small clockwise rotation 
from the moz plane and —@ a small counterclockwise 
rotation, the method for choosing the correct sign of the 
square root is clear. 

See Internationale Tabellen Zur Bestimmung von Kris- 
tallsstrukturen (J. W. Edwards, Ann Arbor, Michigan, 
1944), revised edition, Vol. II, p. 571, 
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greatly simplifying approximations. Also in 
practice the change in r may be neglected over 
the area. of a photographic plate ordinarily 
considered. 


Case II 


For Case II, an average of the contributions 
from all positions of the molecule about axis fq is 
required. This is accomplished by integrating (7) 
with respect to 6, from 0 to 27, the probability of 
the chain occurring in any interval 6 and 6+d6 is 
d6/2m. We note in (6) that cos ag is a function of 6 
and rewrite it as 


cos ag=ky+A cos (6+B), (10) 


where 


9 


A=+(ki2+hk22)! and B=tan-!—. 


ky 


Expression (10) replaces the expression for cos ag 
in (7) and we find that the contribution from all 
positions about axis rz may be obtained from 
considering the real part of the integral, 

exp [tsrako | 


Qn 
xf exp [isrzA cos (6+B) ]d6/2xr (11) 
0 
which has the value, 
cos (srako)Jo(sraA), (12) 


where Jo represents the Bessel function of zero 
order. We thus have from (7) and (12) the desired 
formula II, 


f.? sin? (Nsr; cos aq/4) 
2 





~ 2 sin? (sr, COS a@a/2) 
[1+ cos (sra cos ya COS aa) 
XJo(sra sin ya sin aa) ], (13) 


where cos a, is defined in (5) and N is even. The 
relations required for writing (13) in terms of the 
angular coordinates of the scattered beam, ¥ and 
¢ are found in (8). The term s sin ag can be ex- 
pressed in terms of y and ¢ using s and $s cos ag. 
Since Jo(x) =Jo(—x), no difficulty results from 
the ambiguity of sign arising from deriving the 
sine of an angle from its cosine. 


Case III 


In order to obtain formulae III and IV we 
work with expression (2). The vector rg in the 
above analysis can be replaced by the general 
interatomic vector r;; and then aq is replaced by 
a:; and ya by yi;. We may write instead of 
expression (6), 


cos a;;= Hy+H, cos 6— Hz sin 6, (14) 
where 


Hy=cos a, cos ¥:;, H,=sin y;; cos B, 
and 
H2=sin yi; sin Yn sin (da—Gn). 


We substitute (14) into (2) and find the intensity 
proportional to!” 


fe 


9 


4 


i) 


M= 


N 


> 


1 j=1 





i] 


exp [isr:;(Ho+Hi cos 6— He sin 6)]}. (15) 


To obtain formula III an average is taken of the 
contribution from all values of ¢, from 0 to 27. 
The expressions for cos a, and cos 8 given in (5) 
are substituted into (15) and we find that we 
must consider the integral 


f exp [isr:;C cos (¢’+D) ]do’/2m, (16) 


0 
where 


¢’ =d¢a—n; 
C=+([cos yi; sin PY, sin Wa 
—sin yi; sin ¥, cos Ya cos 6 |? 
+[sin y;;sin y, sin 6 }*)! 


and D is an arc tangent function which disap- 
pears, as previously, in the integration. The value 
of (16) is 

Jo(sri;C). (17) 


Next let us break up the summation over the 
interatomic distances into two groups, one over 
the distances along the two rows of atoms and the 
other over the distances across the rows. If we 
consider N even, it can be shown that there are 

2% Note that since rij=—rji, expression (2) is real and 
thus so is (15). It will be convenient henceforth, however, 
to consider only the contribution from r;; (when 1#7) and 


obtain our final result by multiplying the real part of this 
contribution by 2. 
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2v distances of length (V/2—y)r; along the two from 1 to N/2 and also X distances of projected 
rows where v goes from 1 to (N—2)/2. Acrossthe length —({[N—1]/2—A)r; as X goes from 1 to 
rows there are yw distances projected by the (N-—2)/2. All the distances across the rows, 
cosines of their respective angles y;;,and have the multiplied bysin y;;, give 1 .e inter-row distance." 
projected length ([V+1]/2—y)r: where u goes We get for Case II] 





oe (N—2)/2 


f 
Nf.2/1°+2— > 2vcos[(N/2—»)ri cos Pas cos ~x |Jof (N/2—v)ri sin Pas sin Pn | 
7 v=1 


2 N/2 
+o > wcos [(LN+1 ]/2—pu)r: cos Was cos Hn +7o sin Pas COS Pz COs 6 | 
7? y= 
XJo{s sin ¥n(L(LN +1 ]/2—u)ri sin Wa—70 cos Wa cos 6 |?+[ 79 sin 6 ]?)3} 
Phi 2 (N—2)/2 
+o Pp \ cos [ —([N—1]/2—A)r: cos Was CoS Hn +7o Sin Was COS Pp COs 6 | 
7? 


XJo{s sin ¥.([—(LN—1]/2—A)r: sin Wa—710 cos Wa cos 6 |?+[ 70 sin 6 ]?)'}, (18) 


where 79 is the inter-row distance. The factor 2 appears since all interatomic distances are counted 
twice in the double summation. By using (8), expression (18) may be given in terms of the angular 
coordinates of the scattered beam, ¢ and y. 


Case IV 


It would seem natural to obtain formula IV from formula III directly by integrating with respect 
to 6 to get the average contribution from random positions about rz. However, a formula which 
appears more simple to calculate is obtained if one reverses the order of integrating, that is to derive 
the average contribution from random positions about r, first and then the average contribution from 
random positions about z. We note that (14) may be written 


cos a;;=H)+E cos (6+ F), (19) 
where 
E=+(H,?+H,’)! and F=tan- (AH2/M). 
We have by substituting this into (2) 


2N N 
al > > exp [isri(Ho+E cos [6+ F)]) ]. (20) 


v2 i=1 j=1 


The average over all values of 6 is obtained from considering the integral, 


2r 


f exp [isri;E cos (6+ F) ]d5/2z, (21) 
0 


whose value is Jo(s7;;E). The contribution from the rotation of the chain about axis r, is then 


Jes 
r? 


N WN 
X be exp [tsri; COS Yi; COS a |Jo(s7;; sin yi; Sin aa). (22) 


8 The case when N is odd is readily considered. There are 2y—1 distances of length ([N+1]/2—y)r: along the two 
rows where » goes from 1 to (N—1)/2. Across the rows there are yu distances projected by the cosines of their respective 
angles, y:;, and have the projected length (V/2—y)r: as » goes from 1 to (N—1)/2 and also \ distances of projected 
length —(N/2—\A)r: as \ goes from 1 to (V—1)/2. 
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Let us now split the above summation among distances along the two rows of atoms and those across 
the rows (neglecting for the present terms for which 17=j, the atomic scattering). We get 


f2w—mi fe2nie 
- r? u=1 


fe? (N-2)/2 


2/2 
2— > 2vexp [i(N/2—v)ris cos aq |+2— ¥ wexp [(4(LN+1]/2—u)ris cos a |Jo(sro sin aa) 
“ v=1 


+2— > rexp[—i(LN—1]/2—A)ris cos aa |Jo(sro sin aa). (23) 
A=1 


9 


We now wish to obtain the average contribution 
from all possible values of ¢,. The expression for 
COS @ is substituted into the first summation in 
(23) and it is found that we must consider the 
integral, 


f exp [7(N/2—y)ms sin y, 


0 


Xsin Pa COS (Ga—Gn) \dda/27. (24) 
Its value is 
Ji(LN/2—v]ri sin Pas sin p,). (25) 


For the second and third summations in (23) we 
must obtain a useful expansion for Jo(s7o sin aa). 
We have 


sin a= (1—cos? ag)?! 
=(P-—Q cos (da—Gn) (26) 
—R cos” (da —¢,)), 


where 


P=1-cos? y, cos? Wa, 
QO=2 cos yf, cos Ha SiN Wp SiN Wa, 
R=sin? y,, sin? Pq. 


The above expression may be approximated by 
neglecting R. This imposes a strict limitation on 
the applicability of the final formula. An estimate 
may be made of this limitation by setting ¢a—@n 
equal to 0° or 180° and calculating the values of 
P,Q,and R for some typical cases. It is found that 
if one assumes a value of 20° for angle a, the 
area enclosed in —1.5<x<1.5 cm and 1<y<4 
cm may be covered with an error in the intensity 
probably not greater than 10 percent in the 
vicinity of the main maxima, and usually much 
less. In the region mentioned the approximation 
improves with larger y and smaller x since sin y, 
is very nearly x«/(x?+y?)}. 





We may write 
sin a, =[a?+b?—2ab cos (¢a—@n) |!, (27) 
where 
a’ =(P+(P?—Q?)*]/2 
and 
b?=Q?/[2P+2[ P?—Q?]}*]. 


From the addition theorem for Jo, we get 


Jo(sro sin a) 


=P émJm(S70)Im(srob) cos m(da—on), (28) 


m=0 


where €,, is equal to unity when m=0 and is 
equal to 2 for all other values of m.'* By substi- 
tuting (28) into the second summation of (23) the 
following integral arises from considering the 
contribution from all values of ¢z, 


f ” eip LN 41Y2~pritateaed, 

" Xc0s (Gen) ] 608 m(ba—b)dGa/2. (29) 
Its value is 
(i)"InC(CN+1]/2—n)risin Yes sin Ya]. (30) 


Clearly, a similar integral arises for the third 
summation whose value is 


(4)"Im{—(LN—1]/2—A)ri sin pas sin pp]. (31) 


From the results, (25), (30), and (31) and the 
terms in (23) which were constant with respect to 
the integration over ¢4,, we may now write down 
formula IV. (The contributions from the second 
and third summations in (23) may be combined 
since J,(—z)=—J,(z) if v is odd and J,(—z) 


= J,(z) if v is even.) 


14 My thanks to Professor Otto Laporte for suggesting 
this expansion. 
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fer? f-? W-2/2 
N ~+2 9 


'™ r° 








2f.2 N/2 
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~ 2vcos [(N/2—v)r, cos Pas COS Hn |JoL(N/2—v)r: sin Pas sin Wp |] 
v=1 


+— } (2u—1) cos [(([N+1]/2—p)ri cos Was cos Wn | 
u=l1 


r 


XD (—1)"€mJ 2m(S704) Jam (Sob) Joml (IN +1]/2—pn)r: sin Pos sin yp | 
m=0 


2f 2 N/2 





+- 


r: 


p> (2u—1) sin (((N+1]/2—y)r. cos Yas cos Wp | 


x ¥ (- 1)"*12Toms1(S100) J om41(S00) Jomsil (LN +1 ]/2 —p)r; sin WS sin Vn |. (32) 


m=0 


Expression IV, though quite complex, permits 
calculation of the intensity distribution from a 
chosen structure. Examination shows that, within 
the range of validity mentioned above, the argu- 
ments s7oa and srob are sufficiently small to cause 
the infinite series to converge very rapidly. 

Because of the small value of sin y, for points 
along the y axis, formulae III and IV will be 
indistinguishable from formulae | and I], re- 
spectively. In other words the intensity along the 
y axis is virtually independent of the variable ga. 


Additional Remarks 


The intensity from a molecule with r, vertical 
has been previously obtained.2* For the case 
when the plane of the molecule is rotated 4 
degrees from the noz plane (measured according 
to previous conventions) we have 


f.? sin? (Nris cos f,/4) 





r? sin? (75 COS Wn/2) 
X(1+cos [ras cos Pn COS Ya 
+ras sin Wy, sin ya COS (@n—¢a) }). (33) 


By integrating (33) with respect to ¢a from 0 to 
2m we get 


fe? sin? (Nrs cos W,/4) 





r® sin? (r,s cos W,,/2) 
X (1+ cos [ra cos yas cos Pp | 
XJolrasin yas siny,]) (34) 


for the intensity when there is free rotation about 
r,. Expression (34) agrees with the result of 
Germer and Storks when given in terms of the 





coordinates of the photographic plate, making 
suitable approximations. 


DISCUSSION 


For the purpose of interpreting diffraction 
patterns produced by chains tilted away from the 
normal, the effect on scattering of only one of the 
two rows of carbon atoms in a hydrocarbon chain 
has been used. This has been successful in locating 
the main maxima. However, subtle effects on the 
intensity of the interaction of the two rows of the 
staggered chain could not be anticipated. The 
formulae obtained above which take into account 
the effect of both rows give the scattered in- 
tensity for some of the more interesting spatial 
orientations. They are, of course, not exhaustive. 
For instance, expressions were not obtained for 
the average effect of a variation in angle ya, the 
tilt from the vertical, through a small angle range. 

It was previously noted that fixed random 
positions of the molecules about an axis is 
mathematically indistinguishable from the case 
of free rotation of each molecule about this axis. 
Thus the same pattern can be expected from both 
cases. Another situation which appears at first 
glance indistinguishable from the previous two is 
the case in which there are small patches on a 
surface in which all the molecules have the same 
orientation, but the orientation varies from patch 
to patch. This case, however, would probably 
arise when the density of packing approaches 
that present in a crystal. Thus, such an arrange- 
ment of the molecules could affect the intensity 
of scattering because of their repeated intermo- 
lecular distances. This effect is not included in 
the above formulae. Instead, the latter are 








Phi 
Phi 
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derived on the assumption that the distribution 
of molecules is sufficiently dilute to make 
intermolecular interaction negligible. A lack of 
agreement, therefore, between the scattered in- 
tensity calculated from the theoretical formulae 
and that actually found may be caused by a close 
packing of the molecules. This effect has already 
been suggested.‘ 

When making calculations for contour maps of 
the intensity, it is very helpful to know the region 
in which the main maxima occur. The formulae 
for the intensity, conveniently expressed in terms 
of the plate coordinates, x and y, may then be 
considered in fairly restricted areas.'® Such calcu- 
lations are being made at present in this labora- 
tory for future publication. 





's The following useful tables are available: (a) Tables 
of Sines and Cosines for Radian Arguments, Federal Works 
Agency Work Projects Administration, Project No. 765-97- 
3-10, New York, 1940. (6) Hague, Proc. Phys. Soc. 29, 211 
(1917). (c) Gray, Mathews, and MacRobert, A Treatise on 
Bessel Functions (The Macmillan Company, New York, 
1931), p. 286. 


An improvement in interpretation can be ex- 
pected from a method of photographing which 
would allow for quantitative measurements of the 
intensity. P. P. Debye'® has suggested the use of a 
rotating sector in gas diffraction to eliminate the 
effect of the scattering factors. In our case only 
one occurs, f,. It would appear that hydrocarbon 
films are unusually suited for use by the sector 
method since one sector would suffice for all 
films, the radiation allowed to pass through being 
inversely proportional to (6— F,)?/s*. In addition 
the difficulty of vapors diffusing into the camera 
on long exposure, as occurs in gas diffraetion, 
does not arise. 

For his interest I wish to thank Professor 
Lawrence O. Brockway who has directed the 
work performed in the Chemistry Department of 
the University of Michigan under contract with 
the Naval Research Laboratory. 





16 P, P. Debye, Physik. Zeits. 40, 66 and 404 (1939). 
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A relation of the form 
k=aN(xaxp/d2)*+b 


has been found to hold accurately for a large number of diatomic and simple polyatomic 
molecules in their ground states. (The average deviation of k calculated from k observed for 
seventy-one cases is 1.84 percent.) Here k is the bond-stretching force constant, d the bond 
length, N the bond order, and x4 and x z are the electronegativities of the bonded atoms. 
If k is measured in dynes/cm X10 and d in Angstrom units, a and b have the values 1.67 
and 0.30, respectively, for stable molecules exhibiting their normal covalencies, except those 
in which both bonded atoms have only one electron in the valence shell; for diatomic molecules 
of the alkali metals, Nas, NaK, etc., a and 6 are 1.180 and —0.013, respectively; for hydrides 
of elements having a single electron in the valence shell, 1.180 and 0.040, respectively; and for 
diatomic hydrides of elements having two to four electrons in the valence shell, 1.42 and 0.08, 
respectively. Numerous applications of the relation are made and certain exceptions are 
pointed out. 


DESCRIPTION OF THE RELATION have been found to hold approximately for a 


EVERAL relations! ? between force constants 
or vibrational frequencies and bond distances 


‘P. M. Morse, Phys. Rev. 34, 57 (1929); D. Clark, 
Phil. Mag. 18, 459 (1934); H. S. Allen and A. K. Longair, 
Phil. Mag. 19, 1032 (1935); M. L. Huggins, J. Chem. 
Phys. 3, 473 (1933); 4, 308 (1936); J. J. Fox and A. E. 


limited number of molecules. They, in general, 


Martin, J. Chem. Soc. 884 (1939); G. B. B. M. 
Sutherland, J. Chem. Phys. 8, 161 (1940); J. W. Linnett, 
Trans. Faraday Soc. 36, 1123 (1940); C. K. Wu and 
Chang-Tsing Yang, J. Phys. Chem. 48, 295 (1944). 

(1938) M. Badger, J. Chem. Phys. 2, 128 (1934); 3, 710 
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have two or more parameters that depend on the 
atoms forming the bond, or their position in the 
periodic table. The most widely used of these 
seems to be that of Badger,’ 

ke=(C:5/re)'+di;, 
where c;; and d;; are the parameters and 7, and k, 
the equilibrium bond distance and bond-stretch- 
ing force constants, respectively. 

Remick? has correctly ascribed the increasing 
force constants of the series CH < NH <OH<FH 
partly to the increasing electronegativities of the 
series C<N<O<F, on the basis that the dis- 
placement of H would in turn cause a displace- 
ment of the valence electrons in the potential 
field of the heavier atom. In a vibratory link 
A-—B both A and B are displaced with respect to 
the center of gravity of the bonding electron 
cloud, and it seems that the bond-stretching 
force constant k would vary with the electro- 
negativities of both A and B and could be ex- 
pressed as some function of their product, 


k=f(xaxp). 


It is well known that k increases with bond 
order, that, other things being equal, a double 
bond is stronger than a single bond, a triple bond 
stronger than a double bond; and the bond order, 
which I shall designate as N, may be regarded 
as a relative measure of the electronic cloud 
effective in holding the two atoms together. 

Thus k varies with all the quantities—JN, x,xp, 
and d—and it should be possible therefore to 
express it as a combined function, 


k=f(N, xaxp, d). 


Though the assumption of such a relationship 
may be justified by qualitative considerations, 
a theoretical derivation of the function may not 
be possible. In the present study the function 
is evaluated empirically by the correlation of the 
quantities involved for a large number of mole- 
cules in their ground states. It is not possible to 
check it for the different excited states of one 
and the same molecule because values for NV and x 
are available for ground states only. The form 
of the function is: 


k=aN(xaxp/d*)' +0, (1) 


3A. E. Remick, Electronic Interpretations of Organic 
“ee (John Wiley & Sons, Inc., New York, 1943), 
p. : 


where a and 0 are constants for certain broad 
classes of molecules; N, k, and d are the bond 
order, bond-stretching force constant, and bond 
distance, respectively; and x, and xg are the 
electronegativities of the atoms A and B. 

The constants @ and 6} in the above equation 
have the same values for stable bonds of atoms 
exhibiting their normal covalencies, except those 
in which both atoms have a single electron in 
the valence shell. Most of the simple molecules 
with bonds of this type, for which k, N, and d 
are known with reasonable accuracy, are listed 
in Table I. Their experimentally observed bond- 
stretching force constants are compared with 
those computed from the equation, 


k=1.67N(xaxp/d*)*+0.30. (IA) 


Here k is measured in dynes/cm X 1075 and d in 
Angstrom units. Specific values of N, k, and d 
which were used for the calculations are listed in 
Tables I and VI. The electronegativity values are 
those found in column 3 of Table VIII. 

An examination of the types of molecules in 
Table I will be helpful in judging for what other 
molecules Eq. (IA) may be expected to apply. 
It may not be applied to polyatomic molecules 
in which there are appreciable interactions be- 
tween non-bonded atoms, as it assumes a charac- 
teristic bond-stretching constant which depends 
wholly on the properties of the two atoms form- 
ing the bond. Polyatomic molecules for which 
the valence force field function is inadequate 
may not, in general, be expected to fit the rule. 
Odd molecules or molecular links of atoms with 
unsatisfied valencies also seem to be exceptions. 
Certain classes of these may, nevertheless, be 
made to fit approximately the more general 
relation (1) with a proper choice of values for 
a and b. 

For diatomic molecules of the alkali metals, 
Naz, Nak, etc., relation (1) has the form, 


k =1.180(xaxp/d?2)?—0.013. (IB) 


Table II will illustrate how accurately it holds. 

Diatomic hydrides having one electron in the 
valence shell fit accurately relation (I) in the 
form, 


k=1.180(xaxp/d?)?+0.040. (IC) 
See Table III. 
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TABLE I. Observed and calculated bond-stretching force constants. 








Bond-stretching 

force constant in 

dynes/cm X10-5 
Molecule N k obs. 


Gordy Badger 

k calc. from k calc. from 
k obs. k obs. 
% dev. % dev. 





2.81° 
3.98 
7.508 
2.468 
3.135 
15.73 
9.793 
4.503 
15.6 
5.18* 
7.58! 
3.645 
5.965 
5.45™ 
6.20" 
6.248 
4.86° 
18.1? 
13.4! 
5.004 
3.05° 
7.478 
2.59" 
3.298 
2.78 
4.118 
5.15¢ 
9.678 
3.14 
2.078 
2.398 
1.72 
(.83)s 
(.85 )e 
798 
1.048 
1.18 
948 
22.98 
6.86° 
8.234 
7.67% 
4.82” 
3.33¢ 
4.25* 
3.379 
2.$2* 
2.95? 
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2.59 
2.92 


SOmMANWN OR UCONN SOWUPROSSOWH NOON WODL 


m DO bo DO be tO & 


Average % deviation 17.8 








* The effects of triple bond character on the electronegativity of C 
are taken into account. See discussion in section on electronegativity. 

® Bond distances used in calculations are listed in Table VI, column 
3; electronegativities in Table VIII, column 3; bond orders in column 3 
of this table. 

b From R.S. Mulliken, C. A. Rieke, and W. G. Brown, J. Am. Chem. 
Soc. 63, 41 (1941). 

¢ Omitted in computing average deviation. 

4 Badger gives no constants for this combination of period and row of 
the atomic table. 

e J. Howard, J. Chem. Phys. 3, 208 (1935). The value given is in- 
creased 8 percent to correct for anharmonicity. 

f See reference 29 in text. 

_& Computed from the we values selected by G. Herzberg, Molecular 
ey and Molecular Structure (Prentice-Hall, Inc., New York, 1939), 
p. 483. 

5 B. L. Crawford and S. R. Brinkley, Jr., J. Chem. Phys. 9, 69 (1941). 

iG. B. B. M. Sutherland and D. M. Dennison, Proc. Roy. Soc. 148, 
250 (1935). 

iF. Stitt, J. Chem. Phys. 7, 297 (1939). 

k B. L. Crawford, J. Chem. Phys. 8, 526 (1940). 

'R. C. Lord and D. H. Andrews, J. Phys. Chem. 41, 149 (1938). 

™ Computed from the wi given by D. M. Dennison (Rev. Mod. 
Phys. 12, 175 (1940)) using the formula w1 =(1/27)(k/My). 

= T. A. You Wu and A. T. Kiang, J. Chem. Phys. 7, 178 (1938). 

° J. T. Edsall, J. Chem. Phys. 5, 225 (1937). 


pe P. F. Bartunek and E. F. Barker, Phys. Rev. 48, 516 (1935). 

a See reference 10 in text. 

rK. W. F. Kohlrausch, Raman Spectra (Akademische Verlags- 
gesellschaft Becker and Erler, Leipzig, 1943), p. 184. 

8s G. Glockler and T. Y. Tung, J. Chem. Phys. 13, 388 (1945). 

t J. W. Straley, C. H. Tindal, and H. H. Nielsen, Phys. Rev. 62, 161 
(1942). The value given is increased 8 percent to correct for anhar- 
monicity. 

u L. G. Bonner, Phys. Rev. 46, 458 (1934). 

v Computed from the valence force field formula, 


vs =(1/2m) [k(1/Ma) +((1 —cos a) /Mg))}, 


taking v3 as 1280 cm~! and the vertex angle a measured by electron 
diffraction, 100°. 

~ Computed as described in footnote v of Table I, taking vs as 973 
cm~! and the vertex angle a measured by electron diffraction, 115°. 
( xC. R. Bailey, J. W. Thompson, and J. B. Hale, Phys. Rev. 49, 777 
1936). 

y D. M. Cameron, W. C. Sears, and H. H. Nielsen, J. Chem. Phys. 
7, 994 (1942). The value given is increased 8 percent to correct for 
anharmonicity. 

 C. H. Tindal, W. H. Straley, and H. H. Nielsen, Phys. Rev. 62, 151 
(1942). The value given is increased 8 percent to correct for anhar- 
monicity. 
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TABLE II. Observed and calculated force constants 
and bond distances of alkali metals. 











Force 
con- 
stant in Bond k cal. k calc. d calc. d calc. 
dynes/cm distance by from by from 
Mole- X10-5 d obs. Eq. k obs. Eq. d obs. 
cule k obs.* in A (1B) % dev. (IB) % dev. 
Cse 0.069 4.50> 0.067 2.9 4.44 ) 
Ke .099 = 3.923¢ .096 3.0 3.845 2.0 
Liz .254 2.672° 238 (6.3)¢ 2.558 (4.3)e¢ 
LiCs 109 3.754 Als aa 3.83 2A 
LiK 148 3.304 .147 7 3.29 Pe 
LiRb 129 3.504 131 1.6 3.53 9 
Nae 172 3.079 .173 6 3.096 6 
NaCs 108 3.764 .108 0 3.75 P| 
Nak 130 8=©3.504 128 1.5 3.47 9 
NaRb 120 83.594 .120 0 3.59 0 
Rbz .081 4.22% .081 0 4.22 0 
RbCs .075 4.364 .074 1.3 4.31 1.1 


| 
| 


Average % Average % 
deviation 1.39 deviation 0.86 








* Computed from the we values of G. Herzberg, reference 5 in text. 

b Added radii. See V. Schomaker and D. P. Stevenson, J. Am. Chem. 
Soc. 63, 37 (1941). 

¢ From G. Herzberg, reference 5 in text. 

4 Determined by Schomaker and Stevenson (J. Am. Chem. Soc. 63, 
37 (1941)) from Badger’s rule with its parameters fitted to the similar 
molecules—Liz, Naz, and Ko. 

¢ Omitted in computing average deviation. 


Diatomic hydrides of elements having two to 
four electrons in the valence shell may be 
grouped together and fitted approximately to 


k=1.42(xaxp/d?)*+0.08. (ID) 


Results are shown in Table IV. 

In all the above forms of relation (1), & is in 
dynes/cm X10-5, d in Angstroms, and the x 
values those of column 3, Table VIII. These 
units are used throughout the paper. 


FORCE CONSTANTS 


One use of relation (1) is for the determination 
of force constants and thus for the estimation of 
fundamental frequencies associated with the 
stretching of bonds. Frequently it is a difficult 
or impossible task to decide which of the many 
frequencies of a polyatomic molecule are the 
fundamentals associated with the stretching of 
bonds. From x-ray and electron diffraction 
measurements the bond distances in a large 
number of polyatomic molecules are now known 
with an accuracy of 1 percent to 3 percent. The 
bond distances of many others can be estimated 
with a fair degree of certainty from the additivity 
rule.* Electronegativity values of nearly all the 


‘L. Pauling, The Nature of the Chemical Bond (Cornell 
University Press, Ithaca, New York, 1939), p. 153; 
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common elements are available. Perhaps the 
chief difficulty in using the relation to estimate k 
is the uncertainty of the bond order N for many 
bonds which are of fractional order. 


Accuracy of the Relation 


Table | gives a comparison of the force con- 
stants calculated from Eq. (IA) with those 
obtained directly from the observed spectra. 
Tables II to 1V show similar comparisons for k's 
computed from Eqs. (IB)—-(ID), respectively. 
The differences of the two values, k computed 
and k observed, are for most of the cases listed 
in Tables I-IV no greater than the probable 
error in the obseryed k’s, the average deviation 
for the seventy-one bonds being only 1.84 per- 
cent. Only three of the seventy-five computed k’s, 
those for KBr, KCl, and SiH, deviate more than 
6.4 percent from the observed values, and it 
seems probable that the wide deviation in these 
may be due to errors in the observed k’s. For 
example, the ratios of the observed k’s of KI, 
KBr, and KCl are not of the proper magnitude 
when compared to those of Nal, NaBr, and 
NaCl; whereas the computed values have the 
expected ratios. The w, value used to compute 
SiH is listed as uncertain by Herzberg. 

Though Badger’s rule? is of advantage in the 
estimation of bond distances, force constants 
computed from it differ so much from the ob- 


TABLE III. Observed and calculated force constants and 
bond distances of diatomic hydrides of elements with 
one electron in the valence shell. 











Force 
con- 
stant in Bond k calc k cale d calc. d calc. 
dynes/cm distance by from by from 
Mole 10-5 d obs.» Eq. k obs. Eq. d obs. 
cule k obs.* in A (IC) % dev (IC) % dev. 
CsH 0.467 2.494 0.466 0.2 2.488 0.2 
RbH 514 2.368 514 0 2.368 0 
He 5.759 741 5.791 6 .744 A 
KH .560 2.244 .564 Pe 2.254 4 
LiH 1.026 1.596 1.033 a 1.602 A 
NaH 779 1.888 781 a 1.892 p 


Average % 


Average % ge 
deviation 0.27 


deviation 0.42 








* Computed from we values of G. Herzberg, reference 5 in text. 
b From G. Herzberg, reference 5 in text. 


V. Schomaker and D. P. Stevenson, J. Am. Chem. Soc. 
63, 37 (1941). 

5G. Herzberg, Molecular Spectra and Molecular Structure 
(Prentice-Hall, Inc., New York, 1939), p. 483. 
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served values as to make it of questionable value 
in the computation of unknown k’s. This fact 
has been pointed out previously by Penney and 
Sutherland. Other rules, such as the Clark 
relation,’ in which k is expressed as an inverse 
5th or 6th power of d cannot be accurate for the 
computation of k because the errors in d are 
multiplied several times in the evaluation of k. 
The advantage of the present relation for the 
computation of k is that none of the quantities 
involved in the expression for k are raised to 
high powers. 


Observed Force Constants Used for Comparison 


Some discussion of the observed force con- 
stants used for testing the relation is justified. 
Those of all diatomic molecules and of all M—H 
bonds in polyatomic molecules are for vibrations 
of infinitesimal amplitudes. Since vibration con- 
stants involving the light hydrogen atoms are 
appreciably affected by anharmonicity, it was 
necessary to estimate a correction term for these 
when the values found in the literature were 
not for vibrations of infinitesimal amplitude. 
The correction of +8 percent was determined by 
averaging the anharmonic effects in the M—H 
bonds of the molecules for which the correct 
constants are known. 

The simplest molecule containing a given bond 
of known constants was chosen as representative 
of that bond, and the same link with the same 
bond order in other molecules was not repeated 
except for the case of acetylenic CH. Single bonds 
to acetylenic carbon, =C—, were regarded as 
different from other single bonds to carbon 
because of the effective difference in the electro- 
negativity of the carbon in the two cases. 

The force constants of the diatomic molecules 
were computed with the simple harmonic for- 
mula, 


we = (1/27) (k/u)', 


from the w, values selected by Herzberg.® The k’s 
of polyatomic molecules were chosen from the 
literature referred to in the table and were 
computed by the various authors on the assump- 
tion of a valence force field function, or a slightly 


6 W. G. Penney and G. B. B. M. Sutherland, Proc. Roy. 
Soc. 156, 654 (1936). 


TABLE IV. Observed and calculated force constants and 
bond distances of diatomic hydrides of elements with 
two to four electrons in the valence shell. 














Force 
con- 
stants in Bond k cale. k cale. d calc. d calc. 
dynes/cm_ distance by from by from 
Mole- x10-5 d obs.» Eq. k obs. Ka. d obs. 
cule k obs.* in A (ID) % dev. (ID) % dev. 
AlH 1.62 1.646 1.69 4.3 1.693 2.9 
BH 3.04 1.233 3.04 0 1.233 0 
BaH 81 2.232 78 3.7 2.184 ‘s 
BeH 2.26 1.343 2.20 2.7 1.321 1.6 
CH 4.37 1.120 4.34 Ps 1.115 4 
CaH .98 2.002 .97 1.0 1.978 1.1 
MgH 1.28 1.731 1.34 4.7 1.789 3.3 
SiH (2.48) 1.521 2.15 (13.3) (1.380) (9.3) 
SrH 85 2.146 88 3.5 2.195 2.3 


Average % 
deviation 1.71 


Average % 
deviation 2.58 











®* Computed from we values of G. Herzberg, reference 5 in text. 
b From G. Herzberg, reference 5 in text. 


modified form of this function. This treatment 
assumes that the fundamental vibrations of the 
molecules are controlled entirely by the stretch- 
ing and bending of the chemical bonds, that the 
interactions between non-bonded atoms are neg- 
ligible. Constants for several of the polyatomic 
molecules have been computed by more than one 
researcher. In general the agreement of the 
different values is sufficiently close for the 
accuracy required here. Those for such complex 
molecules as CH;OH, CH3;NHo, CoH6S, CoHeSe 
cannot be regarded as very accurate. Rather 
wide disagreement exists in the literature for 
OC], and OF», because of the fact that different 
workers have assigned the fundamentals differ- 
ently. Sutherland and Penney’ have criticized 
the assignments of Hettner, Pohlman, and 
Schumacher® and have suggested an alternative 
assignment. Relation (IA) agrees with the 
Sutherland and Penney values for OCl, but 
with the Hettner, Pohlman, and Schumacher 
assignment for v3 of OF 2. There is sufficient inter- 
action between the hydrogens of the two CH; 
groups of ethane that a modified form of the 
valence force field is required to compute the 
observed frequencies satisfactorily.* This inter- 
action no doubt causes some deviation from 


7G. B. B. M. Sutherland and W. G. Penney, Proc. Roy. 
Soc. 156, 678 (1936). 

8 G. Hettner, R. Pohlman, and H. I. Schumacher, Zeits. 
f. Physik. 96, 203 (1935). 

9F. Stitt, J. Chem. Phys. 7, 297 (1939). 
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TABLE V. Predicted force constants and vibrational 











frequencies. 
Bond-stretching 

Bond force constant 

length calculated* Expected 

in A from Eq. (IA) in frequency 
Molecule Link d obs. dynes/cm X1075 cm"! 
CsBr CsBr 3.14¢ 0.82 168 
CsCl CsCl 3.06° .87 226 
CsI CsI 3.41¢ 12 137 
F, FF 1.4354 7.93 1190 
Ge2He GeGe 2.41° 1.29 238 
RbBr RbBr 2.94¢ .89 192 
RbCl RbCl 2.89¢ .94 248 
RbI RbI 3.26° 76 159 
SCl. SCl 2.00° 2.98 564 
SbH; SbH 1.71¢ 2.35 
SiH;Cl SiCl 2.06° 2.29 480 
SnH, SnH 1.70¢ 2.28 1960° 
TeBre TeBr 2.49¢ 1.88 
TeCl. TeCl 2.36° pM F ag 
HoTe TeH 1.67¢ 2.68 2300» 








* After this paper was submitted for publication Dr. H. Sponer 
called to my attention a recent measurement of the spectrum of TeCle 
by W. Spinnler (Helv. Phys. Acta 18, 297 (1945)), who obtains 2.18 105 
dynes/cm for the bond-stretching constant. 

* A normal single bond is assumed for each line. 

b The observed frequency should be about 4 percent lower because 
of anharmonicity. 

¢L. R. Maxwell, J. Opt. Soc. Am. 30, 374 (1940). 

4M. T. Rogers, V. Schomaker, and D. P. Stevenson, J. Am. Chem. 
Soc. 63, 2610 (1941). 

e Added covalent radii, see reference 4 in text. 


relation (IA). It is interesting that k of 4.30 
computed from the —C—-C— valence vibration 


(1034 cm) alone!® agrees closely with that 
computed from relation (IA). 


Prediction of Some Unknown Force Constants 


Table V lists a few molecules to which relation 
(IA) may be expected to apply. To my knowledge 
the force constants of these have not been deter- 
mined from their spectra. Their bond-stretching 
constants are computed from relation (IA) and 
are listed in column 5. Vibrational frequencies 
obtained from these are given in the last column. 

In computing the vibrational frequencies the 
molecules GesHg and H;SiCl were treated as 
diatomic molecules with the masses of the 
hydrogen atoms added to the heavier atoms to 
which they are bonded. The valence force field 
formula for the bond-stretching frequency, 1, 
of tetrahedral molecules, AB,, is the simple 
harmonic formula with the reduced mass, uaz, 
being replaced by mz, the mass of B. This allows 
a direct calculation of the symmetrical frequency 


10. G. Bonner, J. Chem. Phys. 5, 293 (1937). 


of SnH, from &snu alone. The parallel frequencies, 
v3, Of the triatomic molecules, SCle and H.Te, 
were computed from the valence force formula" 


1 1 1-—cosa\} 
het oem 
Qn ma Mp 


assuming the vertex angle a to be 100° in each 
case. Here m4 and mg are the masses of the 
central and end atoms, respectively. 


BOND LENGTHS 


Table VI compares bond distances calculated 
from relation (IA) and from Badger’s rule? with 
the experimentally observed values. Lengths of 
other type bonds, computed with Eqs. (1B)—(ID), 
are listed in Tables II-IV. 


Observed Values Used for Comparison 


The experimentally observed bond lengths 
taken from the literature have, on the average, 
an accuracy of the order of 1 percent. Those for 
diatomic molecules and for bonds of the poly- 
atomic molecules which are between hydrogen 
and a heavier atom were with few exceptions 
determined from measurements of spectra. The 
probable error for most of them is well below 
1 percent. Other bond lengths, determined by 
electron diffraction, may deviate as much as 
1 percent to 3 percent.from the correct value. 

The values used for IBr and for the C—Se 
bond of C.H,Se are assumed to be the sum of 
the covalent radii of the bonded atoms. Some of 
the ‘‘observed”’ d’s listed in Table III were, as 
indicated, determined from the vibrational fre- 
quencies by Schomaker and Stevenson” using a 
special form of Badger’s rule fitted to the similar 
molecules Liz, Nae, and Ks. These, therefore, do 
not provide an independent test but merely 
show that Eq. (IB) agrees closely with Badger’s 
rule for this type of molecule. The fact that the 
Schomaker-Stevenson rule also agrees closely, 
except for LiCs, lends credence to these values. 
As pointed out by Schomaker and Stevenson, 
the disagreement with their rule by LiCs may 
indicate an error in w, for this molecule. 


11 J. H. Hibben, The Raman Effect and Its Chemical 
Applications (Reinhold Publishing Corporation, New 
York, 1939), p. 94. 

12,V. Schomaker and D. P. Stevenson, J. Am. Chem. 
Soc. 63, 37 (1941). 
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BOND FORCE 


Accuracy of the Relation for Calculation of 
Bond Lengths 


The agreement of the bond lengths calculated 
from relation (1) with the observed values for 
the bonds listed in Tables II-IV and VI is almost 
as good as the accuracy of the observed values. 
The average deviation of d calculated from d 
observed for the seventy-one bonds is 1.34 
percent. If CoHs, CH;0H, CH;N He, and C2H¢Se 
are omitted, the average deviation becomes 1.29 
percent for the remaining forty-two molecules of 
Table IV and 1.18 percent for the sixty-seven 
remaining ones in all four tables. It will be 
noticed that the agreement for the molecules 
of Table IV is hardly as good as that for the 
molecules in the other tables. 


Comparison with Badger’s Rule 


Though the agreement with the experimentally 
observed bond lengths, on the average, is 
markedly better for relation (IA) than for 
Badger’s rule, the advantage is not nearly so 
great for the computation of bond lengths as for 
the prediction of force constants. A comparison 
of the forms of the two relations will reveal 
why this is true. 

The bond lengths and force constants predicted 
by Badger’s rule in Tables I and VI are pre- 
dominately too small. The average deviation, 
d calculated from d observed, can be reduced by 
approximately 1 percent through a different 
choice of the parameters, c;; and d;;. Apparently 
Badger selected these parameters to give the 
best over-all fit for molecules in certain excited 
states as well as in ground states. The com- 
parisons of Table IV imply that whenever 
Badger’s rule is used with the original c;; and d;; 
values, the predicted bond lengths should be 
increased by about 1 percent. 


BOND ORDERS 


The bond order N as used here represents the 
effective number of covalent (electron pair) or 
electrovalent (ionic) bonds acting between the 
two atoms considered. Bonds like HCl and HF 
which have partial ionic and partial covalent 
character!* are regarded as normal single bonds. 


_8L. Pauling, The Nature of the Chemical Bond (Cornell 
University Press, Ithaca, New York, 1939), Ch. IT. 
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TaBLeE VI. Observed and calculated bond lengths. 











Gordy Badger 
Bond dcalc.* d calc. dcalc. dcalc. 
lengths by from by from 
in A Eq. dobs. Badger’s d obs. 
Molecule Link d obs. (IA) % dev. rule % dev. 
AsHs AsH 1.564 1.57 0.6 1.52 2.6 
BeoHe BH 1.18¢ 1.21 2.5 1.12 S.3 
BeO BeO 1.331f 1.340 a 1.308 1.7 
Bre BrBr 2.284f 2.317 1.4 2.286 on 
CH3Br CBr 1.91¢ 1.86 2.6 1.80 5.8 
C2He ee 1.2045 1.206 Pe 1.171 2.7 
CoH4 CC 1.330! 1.322 6 1.255 5.6 
C2H¢ cS 1.553 1.48 4.5 1.43 Pf 
CH3;:CCH cx; 1.21% 1.18 2.$ 1.17 3.3 
CH3:CCH c*c 1.46* 1.48 1.4 1.39 4.8 
CesHe cx. 1.391 1.40 7 1.31 5.8 
CH;Cl CCl 1.77¢ 1.73 2.3 1.69 4.5 
CH3F CF 1.396™ 1.407 .6 1.359 2.7 
CHa CH 1.093™ 1.100 6 1.033 5.5 
HCN C*H 1.0575 1.058 .l 1.004 5.0 
C2He C*H 1.0575 1.054 oS 1.003 5.1 
CH3NHe CN 1.47" 1.41 4.1 1.41 4.1 
HCN CN 1.1545 1.184 2.6 1.149 4 
H2CO co 1.21° 1.19 1.6 1.20 8 
CH;0H co 1.43P 1.49 4.2 1.40 2.1 
C2H6S CS 1.82¢ 1.82 0 1.74 4.4 
CSe CS 1.548 1.53 6 1.53 6 
CoH¢Se CSe 1.944 2.00 3.1 1.84 5.1 
Cle CICI 1.989! 2.014 1.3 1.960 1.5 
GeH, 7eH 1.48° 1.46 1.3 1.52 aoe 
HBr HBr 1.4108 1.397 09 1.419 0.6 
HCl HCl 1.275! 1.236 5.2 1.299 1.9 
HF HF 9178 919 2 914 3 
HI HI 1.600* 1.603 i ¢ 
IBr [Br 2.474 2.50 1.2 e 
ICI ICI 2.321f 2.323 ott 2.294 1.2 
I2 II 2.667! 2.730 2.4 2.641 1.0 
KBr KBr 2.944 3.18 (8.2)> 2.64 (10.2)> 
KCl KCl 2.794 3.23 (15.7)> e 
KI KI 3.23" 3.17 1.9 © 
NaBr NaBr 2.644 2.71 2.6 ¢ . 
NaCl NaCl 2.514 2.51 0 2.25 10.4 
Nal Nal 2.90" 2.81 3.1 2.59 10.7 
Ne NN 1.095! 1.082 1.2 1.113 1.6 
NH; NH 1.014™ 1.012 2 .983 3.1 
H20 OH .958™ .959 l .944 1.5 
OF 2 OF 1.41¢ 1.37 2.8 1.31 7.1 
OCle Oc! 1.68¢ 1.65 1.8 1.62 3.6 
PH; PH 1.45 1.42 2.1 1.41 2.8 
HS SH 1.334” 1.308 1.9 1.344 a 
H2Se SeH 1.53* 1.51 1.3 1.47 3.9 
H2Se2 SS 2.05@ 2.09 1.9 2.03 1.0 
SiH.« SiH 1.45¥ 1.44 of 1.44 P 


Average % 


Average % ge 
deviation 3.38 


deviation 1.53 








* x value for ethylenic C is taken as 2.83. 

* Force constants used for the calculation are listed in Table I, 
column 4; bond orders in Table I, column 3; electronegativities in 
Table VIII, column 3. 

b Omitted in computing the average deviation. 

¢ Badger gives no constants for this combination of period and row 
of the atomic table. 

4G. B. B. M. Sutherland, E. Lee, and C. K. Wu, Trans. Faraday 
Soc. 35, 1374 (1939). 

© See reference 43 in text. 

f See reference 5 in text. 

& See reference c in Table V. 

bh H. Verleger, Physik. Zeits. 38, 83 (1937). 

iE. H. Eyster, J. Chem. Phys. 6, 580 (1939). 

iL. Pauling and L. O. Brockway, J. Am. Chem. Soc. 59, 1223 (1937). 

k L. Pauling, H. D. Springhill, and K. J. Palmer, J. Am. Chem. Soc. 
61, 927 (1939). 

!'V. Schomaker and L. Pauling, J. Am. Chem. Soc. 61, 1769 (1939). 

m D. M. Dennison, Rev. Mod. Phys. 12, 175 (1940). 

=" R. G. Owens and E. F. Barker, J. Chem. Phys. 8, 229 (1940). 

© D. P. Stevenson, J. E. LuValle, and V. Schomaker, J. Am. Chem. 
Soc. 61, 2922 (1939). 

Pp A. Borden and E. F. Barker, J. Chem. Phys. 6, 553 (1938). 

9 See reference e in Table V. 

( tJ. W. Straley, C. H. Tindal, and H. H. Nielsen, Phys. Rev. 62, 161 
1942). 

’ E. K. Plyler and E. F. Barker, Phys. Rev. 44, 984 (1933). 

tA. H. Nielsen and H. H. Nielsen, Phys. Rev. 47, 585 (1935). 

uL. R. Maxwell, S. B. Hendricks, and V. M. Mosley, Phys. Rev. 52 
968 (1937). 

v C. C. Stephenson and W. F. Giauque, J. Chem. Phys. 5, 149 (1937). 

~ G. Herzberg, Infrared and Raman Spectra of Polyatomic Molecules 
(D. van Nostrand Company, Inc., New York, 1945), p. 489. 

=H. H. Nielsen, privately communicated. 

y G. Herzberg, reference w of Table VI, p. 456. 
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TABLE VII. Predicted bond orders. 








Bond 
order 


Bond-stretching Bond predicted 


force constants in length 
dynes/cm X10-5 in A 
Link k obs. d obs. 


BrC* 4.28 1.798 
CN 16.98 1.A5* 
a §.3° 1.49¢ 
CN 17.5» 1.15° 
CN 5.45» 1.43! 
NC 16.3» 1.18! 
CIC* 5.2 1.67 
CN 16.78 1.15¢ 
C—O 5.4¢ 1.428 


Molecule 
BrCN 
CH;CN 
CH;NC 
CICN 





re Syne Sone Ce or 
SD sw Ye Gh Ah 


11.8° 


1.74 


C=O 


SizHs SiSi 








* x value for acetylenic carbon is taken as 2.83. 
® See reference 6 in text. 

b J. W. Linnett, J. Chem. Phys. 8, 91 (1940). 

¢ L. G. Bonner and R. Hofstadter, J. Chem. Phys. 6, 531 (1938). 
4 T. F. Anderson and A. B. Burg, J. Chem. Phys. 6, 586 (1938). 

¢ See reference h of Table VI. 

fW. Gordy and L. Pauling, J. Am. Chem. Soc. 64, 2952 (1942). 
«J. Karle and L. O. Brockway, J. Am. Chem. Soc. 66, 574 (1944). 


Bonds of higher order having some ionic char- 
acter are treated in a similar manner—i.e., 
electrovalent and covalent bonds are considered 
equivalent and additive in assigning the bond 
orders. 

Considering their chemical properties and the 
valencies of their constituent atoms, the bond 
orders of most of the molecular links listed in 
Table I, except those of fractional order, can be 
assigned with a fair degree of certainty. 


Bonds of Fractional Order 


Pauling, Brockway, and Beach" treat mole- 
cules having bonds of fractional order as made 
up of different resonating structures whose bonds 
are all of integral order, the bond character 
being then determined by averaging the integral 
orders of the bonds between the same pair of 
atoms in the different contributing structures. 
Thus the CC bond of benzene would be halfway 
between a single and a double bond, assuming 
the two equivalent resonating Kékulé structures 
to be the only significant ones. Penney! has 
proposed a treatment which, though still based 
upon the valence-bond method, gives fractional 
orders which are linearly related to bond energies. 

“L. Pauling, L. O. Brockway, and J. Y. Beach, J. Am. 


Chem. Soc. 57, 2705 (1935). 
16 W. G. Penney, Proc. Roy. Soc. 158A, 318 (1937). 


For the CC bond of benzene he obtains an order 
of 1.623. 

Coulson,'* who employs the molecular orbital 
method, defines the order of a bond in an un- 
saturated molecule as 1+, where the unit 
value is contributed by the o electrons which 
are regarded as localized between the two atoms 
forming the bond, and where p is the summation 
of the contribution of all the mobile or unsatura- 
tion electrons. He obtains 1.67 for the bond 
order of the CC benzene bond. 

Mulliken, Rieke, and Brown,!'’? using the 
molecular orbital method, have recently shown 
that hyperconjugation—conjugation with or be- 
tween saturated groups—has an appreciable 
effect upon the orders of bonds in hydrocarbons. 
They predict that normal C—C bonds have 
about 11 percent double bond character and that 
normal C==C bonds such as are found in ethylene 
have about 12 percent triple bond character. 
For the CC bonds of fractional order in Table I, 
I have used the bond order values given by these 
researchers. The agreement of their results with 


TaBLeE VIII. Electronegativities. 
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*®L. Pauling, The Nature of the Chemical Bond (Cornell University 
Press, Ithaca, New York, 1939), p. 64. 
* Single bond to acetylenic carbon. 


14°C A Coulson, Proc. Roy. Soc. 169A, 413 (1939); 
J. E. Lennard-Jones and C. A. Coulson, Trans. Faraday 
Soc. 35, 811 (1939). 

17 R,S. Mulliken, C. A. Rieke, and W. G. Brown, J..Am. 
Chem. Soc. 63, 41 (1941). 








order 


rbital 
1 un- 

unit 
vhich 
toms 
ation 
tura- 
bond 


the 
10WN 
r be- 
iable 
ONS. 
have 
that 
“lene 
crer. 
le I, 
hese 
with 


on 


=— ow we 
wah 


—O~100 | Oo Gn Sate 
oe 00 


a. aes. 2. ah hk OU 


ais | Wo 





rsity 


39); 
day 


Am. 














relation (1A) is pleasing since it lends con- 
firmation to the somewhat surprising hyper- 
conjugation predicted by them. It is unfortunate 
that bond orders computed by their method are 
available only for certain CC bonds. 


Prediction of Some Bond Orders 


Table VII lists several bonds for which stretch- 
ing constants and interatomic distances are 
known, but for which bond orders are regarded 
as uncertain. Relation (IA) is used with the 
known constants to predict the bond orders 
listed in column 7. At present there is no check 
on these other than the method of Pauling, 
Brockway, and Beach" of estimating bond 
character, which is applicable to some. As is 
evident from the above discussion, the fractional 
order predicted by the relation is not strictly 
equal to the bond character as defined by these 
authors but is more nearly equivalent to the bond 
order as computed by Mulliken, Rieke, and 
Brown."" 

From its close similarity to ethane, SigHs may 
be expected to have some hyperconjugation 
which would make the order of the SiSi bond 
somewhat above unity. The order predicted, 
though not accurate, strongly suggests some 
double bond character for SiSi. From the elec- 
tron pair theory the fractional bond order could 
be explained by postulating structures of the 
type 


H H H H 


to contribute significantly to the ground state. 
Resonance between structures of the types 
X—C=N Xt=C=N- 
giving XC some double bond character, and CN 
a bond order somewhat less than 3, could explain 
the bond orders predicted for CH3CN, CICN 
and BrCN. For methyl isocyanide the predicted 
orders indicate some contribution to the ground 
state by structures I] and II] in addition to the 
principal structure I, 
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TABLE IX. Calculation of electronegativities. 








Bond Electronegativity 





Mole- Force constant in length from relation (1) 
cule dynes /cm X107-5 in A Equation 
AB k obs.* d obs.> XA used 
AgCl 1.83 2.42 1.76 (IA) 
AgBr 1.68 2.59 1.89 (IA) 
Agl 1.45 2.81 1.96 (IA) 
AuCl 2.56 2.48 3.10 (IA) 
AICI 2.08 2.423 1.53° (IE) 
AIBr 1.71 2.29 1.53° (IE) 
All 1.31 2.50 1.50° (IE) 
BCl 3.40 1.77 1.91° (IE) 
BBr 2.69 1.94 1.90° (IE) 
BiH 1.70 1.809 1.83 (ID) 
CdH 1.21 1.762 1.07 (ID) 
CuCl 2.29 2.24 2.14 (IA) 
CuBr 2.03 2.42 2.24 (IA) 
Cul 1.72 2.63 2.28 (IA) 
GaCl 1.83 2.11 1.31 (IE) 
GaBr 1.52 2.28 1.33 (IE) 
Gal 1.24 2.50 1.41 (IE) 
HgH 1.14 1.741 97 (ID) 
InCl 1.58 2.29 1.29 (IE) 
InBr 1.36 2.46 1.35 (IE) 
InI 1.11 2.69 1.44 (IE) 
PbH 1.45 1.839 1.51 (ID) 
TICI 1.45 ye 1.19 (IE) 
TIBr 1.26 2.48 1.39 (IE) 
TI 1.07 2.71 1.40 (IE) 
TIH 1.15 1.870 1.13 (ID) 
ZnH 1.51 1.595 1.20 (ID) 











@ Computed from we values of G. Herzberg, reference 5 in text. 

b The bond lengths for hydrides are from G. Herzberg, reference 5 
in text. For the halogens they are determined from the additivity rule 
using the Schomaker-Stevenson corrections. See reference 4 in text. 

¢ These values are included to demonstrate applicability of Eq. (IE). 
Pauling gives x values for Al and B as 1.5 and 1.9. See Table VIII. 


H H 


I. H-C—Nt=C-: I. H- C=Nt=C: 






ELECTRONEGATIVITIES 


Electronegativity represents the attraction of 
an atom for the electrons which bind it in a 
stable molecule. Mulliken'® has shown that the 
average of the first ionization energy and the 
electron affinity of an element is a measure of 
this quantity. Though his is a rather straight- 
forward measure, it has not proved applicable to 
a very large number of elements. 


18 R. S. Mulliken, J. Chem. Phys. 2, 782 (1934); 3, 573 
(1935), 
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Pauling’® has derived a more widely applicable 
method based upon the heats of formation of 
compounds. He shows empirically that 


XA XBR= 0.208A?, 


where xa and xg are the electronegativities of 
the elements A and B; and where A, the energy 
above that expected for a normal covalent bond, 
expressed here in electron volts, is determined 
from the bond energy D of the bonds A—B, 
A—A and B—B, thus, 


A = D(A—B) —3{D(A—A)+D(B—B)}, 


by postulation of the arithmetic mean.?® Mul- 
liken’s values and certain other considerations 
are used to set the origin of the Pauling scale. 

The relation which is the subject of this paper 
offers a means of determining electronegativity 
from molecular spectra. The method is widely 
applicable and appears to be of worthwhile 
accuracy. 


Comparison with Pauling’s Values 


The electronegativity values used to establish 
the relation are those of Pauling.2! When a 
given element was involved in more than one 
bond its x value was adjusted to give the best 
over-all fit to relation (I). The values were 
extended to the third place where the accuracy 
of the data justified it. In no instance was it 
necessary to alter the last significant figure 
given by Pauling—that is, to change the x value 
by more than 0.05. This is remarkable, since 
Pauling’s values were determined in a different 
manner from other bond properties. A com- 
parison is given in Table VIII. 


Determination of Some New Values 


Relation (1) was used to determine a number 
of x values not given by Pauling. These are also 
listed in Table VIII. The molecular constants 
and particular form of relation (1) used for com- 
puting them are shown in Table IX. 

The only molecules involving Ga and In for 
which force constants are available are certain 
diatomic halides. These have unsatisfied valencies 
which make relation (1A) inapplicable, and their 

19 L, Pauling, reference 13, p. 58. 


20. Pauling, reference 13, p. 48. 
*L. Pauling, reference 13, p. 64. 
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interatomic distances have not been measured. 
Consequently, I have chosen the parameters a 
and 6b of relation (I) so as to give the best fit 
for the similar halides of B and Al (elements of 
the same valence) with all interatomic distances 
determined from the Schomaker-Stevenson rule,” 
and have used the resulting equation, 


k=2.25(x4xp/d*)‘—0.20, (IE) 


to determine the x values for Ga, In, and TI. 
The consistency of the results (see Table IX) 
appears to justify the procedure. 

The x values for Ag and Cu and especially 
that for Au seem rather high in comparison to 
those of other metals. The extrapolation of 
Eq. (IA) to these elements is certainly question- 
able. Some double bond character may exist in 
the halides used, and the bond lengths which 
were determined from the additivity rule’ may 
not be sufficiently accurate. Nevertheless, the x 
values for Ag obtained from three different 
halides agree reasonably well. The results are 
likewise consistent for Cu. Values determined 
from the hydrides were not used in assigning x 
for these elements because their polyvalent 
character makes it uncertain which form of 
the relation applies. However, it seems that the 
correct values should be somewhere between 
those given by Eqs. (IC) and (ID): 1.6 and 2.1 
for Ag, 1.7 and 2.3 for Cu, 3.0 and 3.9 for Au. 
The values assigned are in qualitative agreement 
with other properties.” Elements of sub-group | 
—Cu, Ag, and Au—have much higher ionization 
potentials, higher boiling points, form oxides 
that are less basic, etc., than do elements of the 
main group I—Li, Na, etc. 


Effects of Bond Type 


Carbon bound by a triple and a single bond 
shows a significantly higher electronegativity in 
the direction of the single bond than does carbon 
held by bonds of lower order. Possible causes of 
this are the less effective shielding of the nucleus 
by the electrons involved in the triple bond and 
the smaller than normal** single bond radius 
which allows a closer approach of the bonding 

22W. M. Latimer and J. H. Hildebrand, Reference Book 
of Inorganic Chemistry (The Macmillan Company, New 
York, 1940), p. 101. 


3 L. Pauling, H. D. Springhill, and K. J. Palmer, J. Am. 
Chem. Soc. 61, 927 (1939). 
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electrons to the nucleus. This higher value of x 
in ==C— is revealed in other phenomena. Hydro- 
gens joined to acetylenic carbon are known to 
be more acidic™ and to participate in hydrogen 
bonding more effectively than do those of a 
more saturated hydrocarbon.”® In this work the 
x value for =C— was chosen which provided 
the best fit for the CH bond of HCN and of 
CsHe, and this value was used for other single 
bonds to a triple bonded carbon. 

Other variations in bond character—bond 
type, degree of hybridization of the atomic 
orbitals—no doubt influence the effective electro- 
negativities. Their effects, which do not seem as 
pronounced as that of acetylenic C, are neglected. 


Effects of Electric Charges 


The x values listed in Table VIII do not hold 
for charged atoms. However, Pauling”* has given 
simple rules, based on changes in nuclear screen- 
ing, for correcting these values for the effects of 
electric charges. If x* represents the electro- 
negativity of an element with a +1 charge, then 


xt=x+2(x,—x) 


approximately, where x is the value for the 
neutral atom and x, is the value for the atom 
immediately to the right of it in the atomic 
table. Similarly for a negative unit charge, 


x= x—3(v—m), 


where x; is the electronegativity of the element 
immediately to the left in the atomic table. 
These rules are used to correct for the effects of 
formal charges in ionic groups treated in this 
paper. 

For most applications of relation (1), the 
effects of partial ionic character of a bond on the 
x values may be neglected as the electronega- 
tivity product tends to remain constant when 
the x value of one of the bonded atoms is in- 
creased at the expense of the other. 


SPECIAL APPLICATIONS 
(MH)* 
Assuming a normal bond order of unity for 
these ions, higher force constants would be pre- 


*4 J. B. Conant and G. W. Wheland, J. Am. Chem. Soc. 
54, 1212 (1932). 

*%S. C. Stanford and W. Gordy, J. Am. Chem. Soc. 63, 
1094 (1941). 

*6L. Pauling, reference 13, p. 65. 


dicted for them, because of their greater electro- 
negativity product, than for the corresponding 
neutral group. When Mt has a large x value 
compared to H, however, contributions of the 
form, M Ht, with zero bond order, could reduce 
the bond order enough to more than counter- 
balance the effects of the greater electronega- 
tivity product. Thus, understandably, force con- 
stants for (HCl)*+ and (OH)* are lower than 
those of the corresponding neutral group, whereas 
those for (BH)* and (BeH)* are higher. 

From relation (1), the bond order of M+—H 
is approximately 


N’ = (k'/k) (xmd"?/xm’d?)*, 


where the primed constants are for M+—H. 
This assumes that a and } have the same value 
for M—H and M+-—H and that 3 is small com- 
pared to k or k’. The former assumption is 
certainly questionable in some instances, and 
the computed N’s should be regarded as only 
approximate predictions. The bond orders for 
(OH)+, (HCl)+, (BH)+, and (BeH)*+ computed 
in this way, using Herzberg’s® constants, are 
0.6, 0.8, 0.9, and 1.0, respectively. 


(NH,)* and (H;0)* 


The principal structures contributing to the 
ground state of (NH,4)*+ would be I and the four 
equivalent forms of I]. 


H H 


| | 
1. H—N+—H ll. H—N H+ 


| | 
H H 


Pauling?” has predicted 20 percent ionic char- 
acter for NH. I shall assume equal resonance 
between I and each of the forms of II. This 
would mean that a particular NH is bonded only 
80 percent of the time—since no effective bond 
exists in N Ht+—and that N=#. The x value 
of N is increased by about .1 since it has a + 
charge } of the time that a given NH is bonded. 
The NH distance will not be greatly different 
from that in NHs3, and | shall take it to be 1.02A. 
With these assumptions relation (IA) gives 
kun = 5.64. If this value is decreased 8 percent to 
correct for anharmonicity and the simple har- 


27. Pauling, reference 13, p. 72. 
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monic formula applied, the predicted vibration 
frequency is 3070 cm-'. This is in good agreement 
with the absorption of the ammonium halides*® 
in the region of 3100, which presumably arises 
from the valence vibration of (NH,)*. 

A similar treatment can be given (H;0)?*. 
Assuming equal resonance between I and each 
of the three equivalent forms of type II, 


H H 


I. HOH II. a, H+, 
the bond order is effectively $, and the average 
electronegativity of O is increased about 0.1 by 
the + charge. With dou =0.97, the estimated 
kon is 5.96. The ion should absorb in the region 
of 3120 cm~ or 3.2u. This, as it should be, is a 
little beyond the position to which the HO 
band is shifted by the strong hydrogen bonds 
existing in pure water. These ions would exist in 
a solution of a strong acid in water, but it would 
hardly be possible to resolve their absorption 
from that of liquid water and ions of the type 
(H,O)++. A strong acid should, however, tend 
to broaden the 3u water band on the long wave- 
length side. 


CH-NH., CH;OH, C,;H;OH, and C,H;SH 


Linnett?® regards the NH and OH stretching 
force constants which he computes for these 
compounds as lower than normal. If his values 
for methyl amine and methyl alcohol, 6.30 and 
7.65, respectively, are increased 8 percent to 
correct for probable effects of anharmonicity, 
the resulting values, 6.80 and 8.26, are in almost 
exact agreement with the values 6.81 and 8.24 
computed from relation (IA) assuming normal 
single bonds. 

The OH stretching constant obtained by 
Linnett for phenol does, however, appear to be 
significantly lower than that for water or for 
methyl alcohol. Phenol® is known to be appreci- 
ably more acidic than the aliphatic alcohols. 
Pauling* attributes this to structures of Type II 

28 J. H. Hibben, reference 11, p. 466. 

29 J. W. Linnett, Trans. Faraday Soc. 91, 223 (1945). 
The value given is increased by 8 percent to correct for 
anharmonicity. 

80H. J. Lucas, Organic Chemistry (American Book Com- 


pany, New York, 1935), p. 131. 
LL. Pauling, reference 13, p. 189. 





in resonance with those of Type I. 
I. CsH;—O—H, II. C-sH,=O+t—H, 
III. C-sHs==O Ht. 


~The OH bond order predicted by relation (1A) 


is 5 percent lower than unity. Hence, it appears 
that structures of Type III are also of im- 
portance. 

Linnett’s value for the SH bond-stretching 
constant of thiophenol 3.8 (or for CH;SH) after 
the estimated correction for anharmonicity is 
made becomes 4.1, the value predicted for a 
normal SH bond by relation (IA). 


C-=N 


The force constant of C-=N computed from 
the CN vibrational frequency of NaCN® is 
16.5, a value somewhat lower than that for the 
CN group of HCN. The lower force constants for 
the ion can be accounted for without assuming 
any decrease in the bond order, by decreasing the 
x of C by 0.3 and increasing the bond distance 
by 0.05A to correct for the negative charge® 
on C. The value predicted in this way from 
relation (IA), assuming a triple bond, is 16.2, in 
excellent agreement with the observed value. 





Guanidonium Ion 


According to Pauling** the guanidonium ion 
resonates equally between the equivalent struc- 
tures of the type, 


+NH» 
| 
C 
¥. 
H.N NH; 


giving the CN bond 33 percent double bond 
character. Assuming N=1.33, and using the 
kcon=7.1 computed by. Kellner,*® the CN dis- 
tance calculated from (IA) is 1.34A. This seems 
to be a reasonable value. Pauling’s relation*®* be- 
tween bond distance and double bond character, 
assuming 33 percent double bond character, gives 
don = 1.36A. 

#W. Gordy and D. Williams, J. Chem. Phys. 3, 664 
(1935). 

% LL. Pauling, reference 13, p. 159. 

4 LL. Pauling, reference 13, p. 198. 


% L. Kellner, Proc. Roy. Soc. 177, 456 (1941). 
%°L. Pauling, reference 13, p. 164. 
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Urea 


The principal structures contributing to the 
ground state of urea are: 


pu N*H, NH, 


O=C 2 Si, « OE , o-c 
‘s \ \ 
NH: NH, Nt+H: 
Kellner** has computed from spectra kco =9.7 
and kcn = 6.6. Using these values and the inter- 
atomic distance values dco=1.25A and don 
= 1.37A, measured by Wyckoff and Corey,*” the 
bond orders obtained from (JA) are 1.5 and 1.3 
for CO and CN, respectively. Pauling* predicts 
20 percent double bond character for CN and 
60 percent double bond character for CO. 


BoH; 


Because of its unusual and still somewhat 
uncertain structure this molecule is of great 
interest. Lewis** has suggested that the total of 
twelve valence electrons available combines into 
six electron pair bonds which resonate among the 
seven bonds to give 6/7 of a normal electron 
pair bond for each. Pauling*® assumes that in 
resonance with these are several structures of 
the type, 


HH HH 
H-B:B-H, H:B-B-H, 
HH HH 


containing one-electron bonds. Since a one- 
electron bond is approximately equivalent to 
one-half an electron-pair bond, each of the bonds 
would still average about 6/7 of a normal single 
bond. Taking N=6/7 for BH and BB and the 
bond distances measured by Bauer, dgy = 1.27A, 
dgp=1.86A, the force constants computed from 
relation (IA) are kpu=3.1 and kgg=1.8. This 
value of kpn does not agree satisfactorily with that 
of 3.9 computed from the observed spectra.?* 





7 R. W. G. Wyckoff and R. B. Corey, Zeits. f. Krist. 89, 
462 (1932). 

88 G. N. Lewis, J. Chem. Phys. 1, 17 (1933). 

*° L. Pauling, reference 13, p. 241. 

“0S. H. Bauer, J. Am. Chem. Soc. 59, 1096 (1937). 

“ Anderson and Burg (reference d of Table VII) ten- 
tatively assigned the widely separated frequencies 2102 


The B—B vibrational band for kgg = 1.8 should 
occur at about 660 cm. No very strong band is 
observed near this frequency in either the Raman 
or infra-red spectra.” 

Recently Longuet-Higgens and Bell* have 
proposed the bridged structure 


H. H H 


held together with hydrogen bridges, with no 
direct bonding between the borons and with 
normal single bonds between the terminal BH 
groups. Relation (IA) supports their interpreta- 
tion. Using their value, 1.18A for the B—H 
distance, and assuming a normal single bond, 
N=1, the B—H force constant computed from 
(IA) agrees closely with that computed by 
Linnett”® from the data of Stitt. See Table I. 


CO,, SCO, and CS, 


Normal electronic structures with double 
bonded links can be assigned to carbon dioxide. 
Nevertheless, the CO distance, 1.16A, is some- 
what shorter than that for a double bond, and 
Pauling“ has shown that the three structures, 


I. O=C=O, II. O-—C=0', 
II]. Ot=C—O-, 


probably contribute about equally to the ground 
state. Taking account of the effects of formal 
charges on the x value of O and of the triple 
bond on the x of C for the second two structures 
in proportion to their contribution and with N 
equal to the average bond character (2.00), the 
bond-stretching constant from (IA) is 14.4. This 
is in fair agreement with the observed value,*® 


cm and 2523 cm both to a BH, stretching vibration, 
offering resonance splitting as a possible cause of the 
separation. From this questionable assignment they ob- 
tained a value of 3.0 for kpy. They point out, however, 
that the separation of 421 cm™ is abnormally large for 
resonance splitting. 

“F, Stitt, J. Chem. Phys. 9, 780 (1941). See also 
reference d of Table VII. 

48H. C. Longuet-Higgens and R. P. Bell, J. Chem. Soc. 
250 (1943). 

“L. Pauling, reference 13, p. 179. 
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15.2. A similar interpretation can be given 
for SCO. 

The resonance energy of CSe, only 11 kcal./ 
mole** as compared with 33 kcal./mole for COs, 
indicates that structures of Types II and III are 
much less significant in this molecule. I have 
consequently regarded it as having normal 
double bonds. See Table I. 


O,, so, and S. 


One might expect, from their covalencies, two 
normal electron-pair bonds in these molecules. 
Their *Y ground states,’ nevertheless, indicate 
that they have unpaired valence electrons. 
Pauling** has postulated for them a normal single 
bond plus two three-electron bonds. Though the 
two three-electron bonds in these molecules 
should be about equal in energy to a single 
bond,** it is not possible to regard their -bond 
order as 2.00 in the present relation. The bond 
orders required to make the relation fit are 1.4, 
1.7, and 1.9 for Os, SO, and Se, respectively. 
Ses and Tee behave similarly. 


The force constant calculated from (IA), 
assuming triple bonds for these molecules, is too 
large by 8 percent for PN and by 10 percent 
for Ps. It is not probable that these deviations 
are caused by errors in the observed k’s or d’s, 
which are those listed by Herzberg.’ They could 
be made to fit the relation adequately by an 
8 percent reduction in the x value for P. However, 
the x used for P seems to be satisfactory for PHs. 
There is no apparent reason why these molecules 
should be exceptions to the relation since the 
similar molecule Ne is not an exception. PN and 
P, are usually regarded as having pure triple 
bonds. The tendency to form triple bonds, 
however, decreases rapidly with the increase in 
size of the atoms bonded, and tentatively I am 
attributing the misfit of the relation to this 
tendency. 

The necessary reduction in bond order may 
be achieved in these molecules by a small con- 
tribution from structure II, containing a single 


4 L. Pauling, reference 13, p. 131. 
46 L. Pauling, reference 13, pp. 253-4. 
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bond and incompleted octets on both atoms. 


BR Isis: . sa. “tri: . 

The incompletely-filled valence shells should 
allow the second structure to exist with the same 
interatomic distance as that required by |, 
though because of its much higher energy the 
contribution from I] may not be very significant. 


H:0, and H,N, 


In contrast to the satisfactory agreement of 
the similar molecules HeS2 and H¢Ce, both 
H:O2 and H,Ne appear to be exceptions to 
relation (JA). This behavior will be treated in a 
later communication. 


Pyramidal AB; Molecules 


When A is not large compared to B, these 
molecules are complicated by the possibility of 
interaction of the non-bonded atoms. The radius 
ratio 74/rg should give a relative indication of 
this interaction. Consequently, I have listed the 
molecules in Table X in the order of increasing 
values of r4/rg. From an examination of the table 
it is seen that the bond-stretching force constants 
calculated from (IA), assuming normal single 
bonds, agree satisfactorily with the observed 
values for all molecules with r4/rg greater than 
1.5. For those with 74/rg between 1.2 and 1.5 
the agreement is fair except for BiCl3. The 
exceptional disagreement for BiCl; makes the 
observed k appear questionable. 

No observed k’s are available for the iodides. 
Because of the large radius of iodine the radius 
ratio would be unfavorable to a good fit of (IA) 
for them. 

Bond-stretching constants predicted from (IA) 
for SbBr; and BiBrs (radius ratios 1.24 and 1.28) 
are 1.71 and 1.59, respectively. 


Tetrahedral AB, Molecules 


Except when B is the small hydrogen atom, 
there is considerable interaction between the 
non-bonded corner atoms of AB, tetrahedral 
molecules. The valence force function is in 
general unsatisfactory for this type of molecule. 
For some of them Urey and Bradley*’ and 


47H. C. Urey and C. A. Bradley, Phys. Rev. 38, 1969 
(1931). 
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TABLE X. Calculated and observed bond-stretching force 
constants for pyramidal AB; type molecules. 








Bond Bond-stretching kcalc.b & calc. 





length force constantsin by from Radius 

in A dynes/cmX10-5 Eq. k obs. ratio 

Molecule d obs.* k obs. (IA) % dev. vA/’B 
AsH3 1.56 2.81¢ 2.84 1.1 3.32 
PH; 1.45 3.33¢ 3.24 24 2.98 
NH; 1.014 6.86° 6.81 er 2.00 
Bi(CHs)s 2.23 1.754 1.87 6.9 1.90 
Sb(CH3)3 = (2.13) 2.09 1.99 4.8 1.83 
AsF3 Lae 3.92° 3.79 3.4 1.75 
As(CHs)s3 1.98 2.444 2.34 4.1 1.63 
PF; 1.52 4.59¢ 4.66 1.5 1.56 
BiCl; 2.48 1.19¢ 1.80 (51.0) 1.48 
SbCl; 2.37 1.78¢ 1.91 7.3 1.43 
P(CHs)s 1.87 2.144 2.60 21.5 1.43 
AsCl3 2.16 2.03¢ 2.30 14.8 1.27 
PCI; 2.00 2.12¢ 2.63 24.0 1.11 
N(CHs)3 1.47 3.974 4.58 15.4 .96 
PBr3; 2.20 1.63¢ PA EY sak .96 








® The hydrides are repeated from Table I. Bond distances for other 
molecules are from reference c of Table V except those in parentheses, 
which are the added radii. 

b Normal single bonds are assumed. 

¢ See reference e of Table I. 

4d J. E. Rosembaum, D. J. Ruben, and C. R. Sandberg, J. Chem. Phys. 
8, 366 (1940). 

e G. Herzberg, see reference w of Table VI, p. 456. 


others** have successfully used modified valence 
potential functions with terms included to account 
for the interaction of the corner atoms. 

Herzberg*® and Silver® have computed bond- 
stretching constants for several tetrahedral type 
molecules assuming the simple valence force 
system. These values do not agree with those 
computed from (JA) but agree, except when A is 
carbon, with values computed from relation I in 
the form, 





9 


Xaxp\? 
= 3.29( *) —0.40. (IF) 


The results of this comparison are shown in 
Table XI. They imply that the interaction of the 
non-bonded atoms varies in a systematic manner 
for the different types listed. An equally satis- 
factory fit is obtained for the series CFy, CCla, 
CBr4, and C(CHs3), if the constants a and 3} of 
relation (1) are given the values +2.94 and 
— 1.40. 

A force stretching the bond AB in these 
molecules would be opposed not only by the 
chemical binding of AB but also by components 
of the interaction forces of the corner atoms. 


48 N.S. N. Nath, Ind. J. Phys. 8, 581 (1934). See also 
K. W. F. Kohlrausch, reference r of Table I, p. 157. 

“9G. Herzberg, reference w of Table VI, p. 182. 

80S. Silver, J. Chem. Phys. 8, 919 (1940). 


Thus quite understandably the AB stretching 
constants computed from (IA), assuming only 
the forces of a normal single bond, are lower 
than those obtained from spectra using the 
valence force field system. An additional factor 
tending to make the values from (IA) (with 
N=1.00) lower than the spectral constants is 
the probable existence®! of some double bond 
character in the bond AB. 


Boron Halides 


The boron halides are symmetrical planar mole- 
cules with bond distances considerably shorter 
than those expected for normal single bonds. 
Pauling® has predicted the three equivalent 
structures of type, 


xX Xt xX 
| | 

ma” = — os 
, a FF i 


4 


/N 
+X X X X X Xt 


to be the principal ones contributing to their 
ground states. Neglecting all structures except 
the above type and counting in the ionic bond, 
the average bond character is 5/3. Applying 
relation (IA) with N=5/3 gives k=4.06, 4.77, 
and 8.80 for BBrs, BCl3, and BF; in this order. 


. Assuming the valence force system, which holds 


reasonably well for these molecules, Herzberg** 
computes from their spectra 3.66, 4.63, and 
8.83 for kpsr, kaci, and kpr, respectively. 


SOME GENERAL DEDUCTIONS 
1 


If the bond order remains unchanged, a formal 
positive charge on one of the atoms of a bond 
should increase its stretching force constant, and 
a negative formal charge should decrease it. 
This follows from the fact that an increased 
positive charge on either atom would increase 
the electronegativity product. and would also— 
unless the bond order is lowered—tend to make 
the bond length shorter. 

The implication of the above statement is 
that when, as frequently happens,”® a positive 





5tL. Pauling, reference 13, Ch. VII. 
8]. Pauling, reference 13, p. 219. 
83 G. Gerzberg, reference w of Table VI, p. 178. 
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charge actually lowers the force constant, it also 
effectively reduces the bond order. That a de- 
crease in the electron density of one of the bonded 
atoms would tend to reduce the bond order is 
not surprising. For illustrations, see the sections 


on (MH)* and (NH,)t. 


2 


For bonds of different orders between the 
same pair of atoms the force constant k,, for any 
given bond order N,, is related to that k, for 
some other order NV, by 


km—0.30 Nn 
kn —0.30 — 0.30 “(5 *) 
This follows obviously from relation (IA). More 
approximately, 
km Nmfdn\} 
eae 


In the above equations d,, and d, represent the 
interatomic distances for the respective bond 
orders, V,, and N,. 

Frequently it is desirable to know the approxi- 
mate ratios of the force constants of single, 
double, and triple bonds between a given pair of 
atoms. If the interatomic distances have not 
been measured for the different bonds, they can 
be determined approximately by adding the re- 
spective single, double, or triple covalent radii 
given by Pauling.** Whenever k is known for 
one of the orders, the first and more accurate of 
the above forms can be applied. Thus the bond- 
stretching force constants for single, double, and 
triple CN bonds should have the approximate 
ratios, 1.00:2.37:4.10. Probably a more useful 





LL. Pauling, reference 13, p. 154. 
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TABLE XI. Calculated and observed bond-stretching force 
constants for tetrahedral AB, type molecules. 








Bond-stretching 
force constant in 





Bond dynes/cm X1075 k calc. 
length k calc. from 

in A from k obs. 

Molecule d obs.* k obs.> Eq. (IF) % dev. 
GeBr, 2.32 2.58 2.56 0.8 
GeCl, 2.08 3.27 3.28 oil 
SiBr, 2.19 2.92 2.94 8 
SiCl, 2.02 3.75 3.63 Sut 
Si(CHs),4 1.93 3.31 3.45 4.2 
iF, 1.54 7.16 7.10 8 
SnBr, 2.45 2.28 § 1.8 
SnCl, 2.30 2.80 2.77 1.1 
TiCl, 2.21 $.11 2.82 9.3 








* SiBr distance in SiBrs is assumed to be equal to that in SiBrCls; 
SnBr distance is determined from the Schomaker-Stevenson rule, 
reference 4. Other distances are from Maxwell, reference c of Table V. 

b k for Si(CHs)3 is from Silver, reference 50. k's for the other molecules 
are from G. Herzberg, reference 49. 


application is the prediction of bond orders from 
force constant ratios. 


3 


The concept of electronegativity, which in the 
past has been used principally in predicting 
properties associated with the partial ionic char- 
acter of bonds, is of real importance in the 
interpretation of molecular spectra. Previous 
applications have been primarily concerned with 
the quantity, x,4—xs, a fact which apparently 
led Wheland®® to comment that it is only the 
differences between electronegativities which are 
significant. In the present work, wide application 
has been made of the electronegativity product, 
XAXB. 

I wish to thank Vida Miller Gordy for assisting 
with the calculations and the preparation of the 
manuscript. 


55 G. W. Wheland, The Theory of Resonance (John Wiley 
& Sons, Inc., New York, 1944), p. 84. 
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In the first paper with this title, a relation was derived from the available experimental 
data giving the Debye characteristic temperature of solid argon as a function of a (the inter- 
atomic distance) or V (the molal volume of the solid). This relation was not consistent with 
the usual Griineisen assumption, which for a small range of volumes may be written approxi- 
mately —@-'d@/dV=constant. In this paper it is shown that the Griineisen assumption, with 
the constant set equal to 0.18 mole cc™ is reasonably consistent with the experimental data. 
On this basis a new relation between the potential energy E, of solid argon and a is derived. 
Most of the earlier conclusions are unaltered, but the present basis seems more reasonable, 
and, as will be shown in the following paper, avoids difficulties when application is made to 


liquid argon. 





N a paper of the same title! (referred to herein- 
after as 1) we have used the available thermal 
and equation-of-state data on solid argon to 
determine the potential energy, EZ,, and the 
Debye characteristic temperature, 9, of the solid 
as functions of a, the equilibrium distance be- 
tween nearest neighbors in the lattice. The rela- 
tion thus obtained between © and a has the 
peculiarity that d@/da varies extremely rapidly 
with a. In fact, it gives dO/da=—211 degree 
per angstrom at a=3.86A (with © itself, of the 
order of 80°), and d0/da=0 at a=3.924A. This 
rapid variation seems extremely improbable, and 
when an attempt is made to apply such a func- 
tion to a study of liquid argon it leads to diffi- 
culties. To be sure, the experimental data on 
which our expression for © was based extend 
only over a range of a from about 3.82 to 3.86A, 
so that anything outside this range may be said 
to be an extrapolation; however, even in this 
range, —d@/da varied from about 273 to 211. It 
appears quite possible that this apparent, rather 
violent variation of dQ/da may be due to experi- 
mental error. In fact, as explained in the follow- 
ing article,? an examination of Bridgman’s data* 
indicates that they do not conform exactly to 
the values given by him for the volume of the 
liquid in equilibrium with the solid at various 
temperatures, and a new set of values has been 
estimated. These give rise to new values for the 
10. K. Rice, J. Chem. Phys. 12, 289 (1944). Erratum: 
Sign inside bracket on right-hand side of Eq. (4) should be 
+ instead of —. 


2Q. K. Rice, J. Chem. Phys. 14, 324 (1946). 
* P. W. Bridgman, Proc. Am. Acad. 70, 1 (1935). 


volume of the solid in equilibrium with the liquid, 
which differ appreciably from those given by 
Bridgman. 

Since our work was based on these volumes, 
some alteration in the results is obviously indi- 
cated, and it also appears that the margin of 
error may be somewhat greater than previously 
assumed. It, therefore, now appears reasonable 
to try the simple Griineisen assumption, 


d1n® /d1|1n V=constant, 


where V is the molal volume; actually, since our 
range of V is small, and it is just as good theo- 
retically and much more convenient, we set 


—O-d0/dV=T, (1) 


where I is a constant. This is also nearly equiva- 
lent to setting dQ/da constant, in view of the 
small range of a and 9 involved, and so in this 
respect the behavior of © is quite different from 
that supposed previously. The definition of the 
functional dependence of © is completed by 
setting 


0 =70.56° at V=24.493 cc/mole 
(a=3.8600A). (2) 


These equations replace Eq. (6) of I. Such a 
relationship must have an empirical basis, and 
can be considered as satisfactory only if it is 
consistent with the experimental data. To test 
this consistency we make use of Eq. (5) of I, 
which may be rewritten in the form 


—P=0"(d0/dV)((E.+ E,.) —(E:+ E.)p-0]. (3) 
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TABLE I. Test of relation between © and V. 
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V(obs.), cc/mole @(calc.) V(calc.) @(calc.) V(calc.) 

TT: P cal./cc Bridgman Revi r=0.18 T =0.18 T =0.195 Ir =0.195 
83.78 0 24.507 24.50t 70.62* 24.489 70.62* 24.489 
126.3 46.9 24.49 24.43 71.31 24.434 70.24 24.516 
161.9 93.7 24.37 24.17 73.86 24.239 71.88 24.398 
192.9 140.6 24.25 24.01 ry Bi 23.999 74.66 24.203 








* Calculated from Table I of I. _ 
+ From Baly and Donnan; see reference 4. 


In this equation P is the pressure at any given 
value of V and of the temperature, 7, and E, 
and £, in the first parenthesis are the thermal 
and zero-point energies at this V and this tem- 
perature, while in the second parenthesis they 
are evaluated at the same value of V but at 
zero pressure and whatever temperature corre- 
sponds thereto. 

To make use of Eq. (3) we note, first, that 
- Table I of I gives an empirical relation between 
(E.+£,)p-o and ©, both of which are determined 
there directly from the experimental data. On the 
other hand, the value of £,+£, can be obtained 
from the Debye theory as a function of 0/T. 
These two relations, together with Eq. (1), 
enable us to express the pressure P at any given 
T as a function of ©, and, hence, through Eq. 
(1), as a function of V (or its equivalent a). 
The result of this procedure is thus to give a 
relationship between P, V, and 7, which can 
be used to calculate V for the solid at values of 
P and T corresponding to various points along 
the melting curve. These calculated V’s can 
then be compared with those given by Bridgman 
or with our revised values. The results are given‘ 
in Table I. It will be seen that calculations of V 
based on '=0.195 agree well with Bridgman’s 
original estimates of V, while calculations based 
on '=0.18 agree well with my revised values. 


4 The values of a given in Table II of I were slightly 
altered from those calculated directly from the values of 
V given by Bridgman, because the estimated value 
V=0.616 cc/g (see end p. 291 and top p. 292 of I) does 
not correspond exactly to the value of a at the melting 
point obtained from Table I of I. In Table I of the present 
paper the values of V(obs.), Bridgman, at the higher tem- 
peratures correspond exactly to those given by Bridgman, 
but the value at 83.78° is obtained from the measurements 
of the density of the liquid by E. C. C. Baly and F. G. 
Donnan (J. Chem. Soc. London 81, 907 (1902)), and the 
measurement of A V of fusion by K. Clusius and K. Weigand 
(Zeits. f. physik. Chemie B46, 1 (1940)). Bridgman’s value 
of AV of fusion at 83.78 is probably in error, but his value 
of the density of the liquid is close to that of Baly and 
Donnan. 





Accepting the latter, it is clear that the results 
obtained by the use of ' =0.18 are certainly well 
within the limits of experimental error. It may 
seem surprising that, by fitting Bridgman’s V 
exactly as we did in I, we obtained results so 
different with respect to the behavior of dO/dV, 
but this happens because all the other properties 
are very sensitive to small changes in the volume. 
On the other hand, experimental measurements 
of the volume are not correspondingly highly 
accurate, so that the representation given by 
Eq. (1) appears to be quite adequate. 

The agreement with the measurement of 
(0P/dT)y at 73°K and zero pressure, made by 
Simon and Kippert® is scarcely altered as com- 
pared with that in I. Calculating (@P/dT)y 
from the relation 


(0P/dT)y= —O-'(dO/dV)C,=1TC,, (4) 


(see Eq. (7) of I), using [T'=0.18, we find 
(0P/dT)y=42 atmos./deg.; in I it was evalu- 
ated as 43 atmos./deg. The experimental value 
is 38 atmos. /deg. 

In spite of the differences between Eq. (1) of 
the present article and Eq. (6) of I, the con- 
clusions drawn in the earlier work are not 
substantially modified. In the experimental range 
of volumes involved in Bridgman’s measure- 
ments of the solid, —@-—'dO/dV as calculated 
from I, ranges from 0.157 to 0.193, of the same 
order of magnitude as the 0.18 found here.® 
On the other hand, the calculations of my earlier 
article’ in which an attempt was made to 
establish a theoretical connection between 0 
and E,, made it appear that —9-'d@/dV could 


5 F. Simon and F. Kippert, Zeits. f. physik. Chemie 135, 
113 (1928). 

6 This means that the total expansion of the solid 
between O°K and the melting point is scarcely changed 
from that shown in Table I of I. 

70. K. Rice, J. Am. Chem. Soc. 63, 3 (1941). 
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not be greater than about 0.10 in the experi- 
mental range, which is a very significant differ- 
ence. It was already pointed out in this earlier 
article that these results were not consistent 
with the work of Simon and Kippert, nor are 
they consistent with Bridgman’s data. i 

In I it was concluded that the previously de- 
rived’ theoretical relation between © and E, 
probably was not valid because the vibration of 
an atom in the fields of its neighbors would be 
markedly anharmonic, even when only zero- 
point vibrations were involved. This situation 
arose because the curvature of the E, vs. a 
curve changed rapidly. This conclusion is not 
invalidated by the results of the present paper, 
when they are combined with Table I of I, to 
give a new relation between £, and a. 

The direct method of getting this relationship 
would be to use Eqs. (1) and (2) above to get V, 
and hence a, as a function of 0. In this way we 
could assign new values for a for the temperatures 
in Table I of I, and so obtain a new relation 
between E, and a. However, in the course of 
some other calculations being made in collabora- 
tion with Mr. G. W. Murphy, it appeared that a 
much more accurate relation could be obtained 
by an indirect method, using the Griineisen 
equation for the pressure [Eq. (3) of I]. For the 
special case of zero pressure, this equation 
becomes 

0O=l(E.+£.) —dE,/dV. (5) 


From Table I of I, we can find £,+£, as a func- 
tion of ©. The following empirical equation 
expresses this relationship quite well: 


E, +E, = 296.7 +29.5(80—@) 
—0.643(80—@)2. (6) 


Eq. (6) comes directly from the experimental 
data and the Debye theory, and is not affected 
by the dependence of © on a. With any given 
value of I’, however, it can be changed into a 
relation between E,+£; and a, and inserted into 
Eq. (5); at the same time dV may be expressed 
in terms of da. We then obtain dE,/da as a 
function of a. It is convenient to set up a new 
empirical expression for dE,/da. This has been 
done, and it has been integrated to give 


— E, = 2007.6 — 1410.0x — 10,710x2 
+53,947x? —76,500x', (7) 


where x =a—3.84. The greater accuracy of this 
method of obtaining the expression for E, arises 
because the absolute error of E,+£; in Table I 
of I is about the same as that of Z,. Small but 
changing errors in E, can make appreciable 
differences in dE,/da; on the other hand, such 
small errors are a very small percentage of 
E,+£; which can, therefore, be used in Eq. (5) 
to get a much more accurate value of dE,/da. 
This, when integrated, using an integration con- 
stant to match E, in Table | of I at one value 
of © (correlated with x from Eqs. (1) and (2)), 
should give a better expression for E, than that 
obtained directly. 

We may now compare Eq. (7) above with 
Eq. (8) of I. The coefficients of x* and x* in 
Eq. (7) are much smaller than the corresponding 
ones of Eq. (8) of I, while the coefficient of x’, 
with which the other coefficients must be com- 
pared, is larger. Nevertheless, the rate of change 
of curvature of Eq. (7) is still quite substantial. 
Thus the value of d?E,/da? obtained from Eq. (7) 
changes by about 150 percent in a range of a of 
about 0.1A. However, this estimate is based on 
the third-degree term in the expression for E>, 
and when we consider the motion of an atom 
surrounded by a highly symmetrical distribution 
of neighbors, the effects of any term of odd order, 
in causing anharmonicities in the effective field 
in which the atom vibrates, tend to cancel out in 
first approximation. It is, therefore, desirable to 
consider the fourth-degree term as well. This is 
also considerably less than the fourth-degree 
term in Eq. (8) of I, but by no means negligible. 
In fact, the fourth-degree term is itself sufficient 
to cause a change of more than 40 percent in 
@E,/da* over a range of a of 0.1A. It is, of course, 
true that the fourth-degree term can be de- 
termined only roughly. But, since greater dis- 
placements than 0.1A may be expected, even for 
the zero-point vibrations, there seems to be no 
reason to doubt that the reasons given in | for 
the breakdown of the previously derived’ relation 
between © and E, are substantially correct as 
there given. 

Over the range of a from 3.82A to 3.86A, 
which about covers the expansion of the solid 
from 0°K to the melting point, the appearance of 
the E, vs. a curve is not greatly changed from 
that given by Eq. (8) of I. (See Fig. 1 of I.) 

































Equation (7) above gives a value of —E, about 
2.2 calories per mole lower at 3.82A, about 0.7 
calories lower at 3.83A, about 0.8 higher at 
3.84A, about 1.4 calories higher at 3.85A, and 
about 0.2 calories higher at 3.86A. Outside this 
range, where the E, vs. a@ curves represent 
extrapolations, the differences are naturally 
greater. In particular, the minimum of £, (and 
correspondingly the minimum of the interatomic 
potential curve) given by Eq. (7) is displaced by 
approximately 0.01A toward smaller values of a 
as compared to that given by Eq. (8) of I. 
However, these differences can hardly affect the 
qualitative relationship between the curves shown 
in Fig. 2 of I. It should be emphasized that the 
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most important difference between the two 
curves shown in Fig. 2 of I is their difference in 
slope in the experimental range of a. It will be 
clear that this difference is scarcely altered by 
the present work. We, therefore, conclude here, 
as in I, that the slope of the E, vs. a curve in 
the experimental range, obtained by direct appli- 
cation of the Debye theory to the experimental 
data, is much greater than that expected from 
the attempt’ to derive a relation between £, 
and 0; and, as stated above, that the attempted 
derivation fails because of the extreme anhar- 
monicity of the interatomic potential curve, 
which means that the theory of small vibrations 
cannot be applied. 
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The thermodynamic properties and equation-of-state 
data of liquid argon are interpreted in terms of a partition 
function previously derived. This partition function has a 
part due to vibratory motion of the atoms in the liquid 
and a part due to translational motion. The vibratory part 
contains the Debye characteristic temperature © for the 
solid, evaluated for a certain interatomic distance do 
which appears as a parameter in the partition function. 
ao is interpreted as the average interatomic distance 
between nearest neighbors in the liquid, and is a function 
of a, the latter being defined as (24V/N)!, where V is the 
molal volume and N Avogadro’s number. With V set 
equal to the molal volume of the liquid at its normal 
melting point, do is assigned a value which appears reason- 
able from the atomic distribution curves obtained by x-ray 
studies of the liquid. With this value of ao, the partition 
function for this fixed value of V gives reasonably good 
values of the entropy over a large range of temperatures. 
A small, entirely reasonable change in the partition func- 
tion improves the agreement with the experimental data. 
The parameters in the partition function have been 
evaluated for other values of V by using equation-of-state 
data. This leads to a calculation of the specific heats of 
saturated liquid argon over a range of temperatures. The 
agreement with experimental data is not exceedingly close, 
but is probably within experimental error. A calculation 


1. THE GENERAL THEORETICAL EQUATIONS. 


N a paper entitled “On the Statistical 
Mechanics of Liquids, and the Gas of Hard 





is also made of the free energies of the solid and liquid 
phases along the melting curve. They are found to be 
equal to within reasonably close limits. The application to 
various values of V leads to the evaluation of ado as a func- 
tion of a. The experimental curve for (0S/8V)r plotted 
against V appears to show a change of slope near the value 
of V characteristic of the liquid at its normal melting point. 
There is also a change in slope of the curve for dao/da vs. a, 
at the corresponding value of a. These facts are interpreted 
in terms of the atomic distribution curves, which appear to 
indicate a change in structure of the liquid which starts 
near the same value of V or a. The liquid appears to 
change from a distorted face-centered structure to a 
distorted body-centered structure as the volume increases. 
Information about the potential energy of liquid argon is 
obtained from the experimental data by application of the 
partition function. It is found that the potential energy of 
the liquid can be rather effectively correlated with the 
potential energy of the solid. Finally, certain peculiarities 
of the equation-of-state data are believed to indicate a pre- 
freezing phenomenon, in which particles containing appre- 
ciable numbers of atoms and having the structure of the 
solid, appear in the liquid. A state of this sort can be 
stabilized by the entropy of mixing of the solid particles 
with the liquid. 


Elastic Spheres’”’ ! (referred to hereafter as SML) 
| proposed the following partition function for 


10. K. Rice, J. Chem. Phys. 12, 1 (1944). 
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LIQUID 


the thermal part of the energy of a monatomic 
liquid in a not too expanded state, and at a 
temperature such that classical considerations 
apply : 


Q=Qo° 
=[(2e)!T'/O+ (2emkT/h?)}2—’e(a—ay) ]®. (1) 


In this equation Q is for one atom; Qo, which is 
introduced here for convenience, is for one 
degree of freedom. The symbols used have the 
following significance : 

T =absolute temperature, 

@=Debye characteristic temperature of the solid, evalu- 
ated for such a degree of compression of the solid 
that the distance between nearest neighbors is do, 

m=miass of an atom, 

k=Boltzmann’s constant, 

h=Planck’s constant, 

a=the distance between neighbors which would be found 
if the liquid were a close-packed cubic lattice having 
the same density. 


We have 
a'=2!V/N, (2) 


where V is the molal volume and N is Avogadro’s 
number. The quantity dp is a parameter; do is 
always less than a. For the present purposes it is 
to be considered as a function of a (though when 
a more exact theory becomes possible, it may be 
found better to take it as a variable parameter 
to be found by minimizing the free energy). 
The interpretation of a) and the justification 
for the use of Eq. (1) are explained in SML. It 
there appears that if we had a gas or liquid of 
hard elastic spheres of radius ¢ we could write 


Qo = (2amkT /h?)32-te(a—0). (3) 


In the case of the real liquid with a force field 
between the atoms, part of this quantity is 
replaced by the term (2e)!7/@. ay may, there- 
fore, be described as a parameter. which divides 
the space available to any degree of freedom into 
a part where the motion is translational and a 
part where it is vibrational. As pointed out in 
Section 4 of SML, we can consider two types of 
motion of atoms, one a long range motion in 
which the atoms move relatively great distances 
and exchange places, and one a short range 
motion. The first type is translational and the 
second type vibrational; any one degree of 
freedom, however, rapidly and frequently shifts 
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from one type to the other. Both types of motion 
are possible even in the case of hard elastic 
spheres, though in this case the “‘vibrational’”’ 
type of motion is actually force-free. The vibra- 
tional motion is the short range type, one atom 
between its neighbors, not involving the longer 
motions involved in exchanges of position. It, 
therefore, seems reasonable to identify do, at 
least in simple cases, and at least approximately, 
with the actual average distance between centers 
of molecules which can be described as nearest 
neighbors in the liquid. Thus we understand that 
© is to be considered as a function of ao, and we 
may expect the potential energy of the liquid to 
depend on dy. This, of course, is not an exact 
statement; it may, however, be expected to con- 
tain a considerable element of truth, and does in 
fact lead to interesting correlations of the 
properties of the liquid with those of the solid, 
as we shall see. 

We shall find it convenient to rewrite Eq. (1) 
in the form 


Qo=LT/0+ MT (a —A), (4) 


where the meanings of Z and M are defined by 
the equations. We can then concisely define 
certain useful quantities: 


B=T(0Q)/dT)y =LT/O+4MT'(a—a), (5a) 
D=T(0B/0T)y =LT/O+1MT'(a—a). (5b) 
Also we have 


(0Qo/da) r= — (LT /0?) (dO/day)day/da 
+MT'(1—dao/da). (6) 


As noted above, the partition function takes 
into account only the thermal energy. The 
thermal energy, by the normal process of loga- 
rithmic differentiation of the partition function, 
is found to be, per mole (i.e., per 3N degrees of 
freedom), 


E.=3RTB/Qo. (7) 
The molal work function, A, is given by 
A=-—3RT In Qo+Ep,, (8) 


where E, is the potential energy per mole and 
is a function of a only. E, is the energy the 
system would have if, without change of volume 
and without crystallization, it were reduced to 
absolute zero, and if there were no zero-point 
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TABLE I. Properties of liquid argon at its normal melting point. 











Property Source Value 
a Baly and Donnan*®» 4.038A 
a Mathias et al.°4 4.048A 
a Bridgman? 4.038A 
Used in this paper‘ 4.037A 


a 
AV, fusion 
AV, fusion 


Clusius and Weigand« 
Bridgman °¢ 


3.53 cc mole™ 
3.18 cc mole 


al = V1(0V/dT)p | Baly and Donnan*» 0.0045 deg.“ 
a Mathias et al.*4 0.0042 deg. 
C, Eucken and Hauck>:* 5.5 cal. mole deg. 
Ce Eucken and Hauck*:® 10.5 cal. mole deg. 
Cp» Clusius® } 10.0 cal. mole deg. 
(@P/dT )v ’ Bridgman® 0.556 cal. deg.-! cc“ 
(@P/dT )v Simon and Kipperti 0.60 cal. deg.~! cc 
AS, fusion Clusius and Weigand 3.35 cal. mole deg.~! 
AS, fusion Bridgman*® : 3.12 cal. mole deg. 








* Involves slight or quite certain extrapolation. 

b See reference 8. 

¢ Involves appreciable extrapolation. 

4 See reference 7. 

© See reference 6. 

f Consistent with V(calc.) of Table I of reference 9. 


energy. The zero-point energy does not need to 
be considered at all, in our case, since we deal 
only with temperatures sufficiently high so that 
the system may be thought of as classical. The 
energy E; does actually include some potential 
energy in virtue of the vibrations of the atoms 
about their positions of equilibrium in the field 
of their neighbors. E, might perhaps be more 
exactly defined as the lowest possible potential 
energy for the liquid at the molal volume V. 
The various thermodynamic quantities can 
now be obtained in terms of the quantities 
already defined. We have for the molal entropy, 


S=—(0A/dT)y=3R In Qo+3RB/Qo, (9) 
and for the pressure 


P=-—(0A/dV)r 
=(3RT/Qo)(0Q0/8V)r—dE,/dV. (10) 


Also 


(0S/dV)r=(0P/dT)v 
= (3R/Qo)L(8Q0/8V) r+ (0B/dV)r 
— (B/Qo)(8Qo0/AV)r J, 


and for the specific heat at constant volume 
C.=3R[(B/Qo) + (D/Qo) — (B?/Qo’) J. 


Derivatives with respect to V may be replaced 
by derivatives with respect to a by means of Eq. 
(2), and Eq. (6) and a similar equation for 
(0B/da)r may be used to get explicit expressions. 
For future reference we write down the expres- 


(11) 


(12) 


& See reference 11. 

bh A, Eucken and F. Hauck, Zeits. f. physik. Chemie 134, 161 (1928). 
i See reference 10. 

i F. Simon and F. Kippert, Zeits. f. physik. Chemie 135, 113 (1928). 
For numbered references see citations in text, 


sion for P: 














RTT LT dO day dao 
=— |- eur 1- 
a’*NQoL ©? dao da da 
2: dE 
_ ?. (13) 
3a°N da 


It is now proposed to use the equations just 
developed as a guide in the study of the proper- 
ties of liquid? argon. In doing this we shall adhere 
to the philosophy followed in a number of earlier 
papers’ in which the approach was frankly semi- 
empirical. The equations do not suffice in them- 
selves for a complete description of the liquid 
without reference to the experimental data, since 
we do not have sufficient information concerning 
some of the quantities involved in them. In par- 
ticular, we do not know dp as a function of a, 
since we do not have a sufficiently detailed 
knowledge of the complicated geometry of a 
more or less random distribution of atoms. It 
does appear, nevertheless, that the equations 
offer a reasonably good guide for the handling 
of the experimental data, and for finding rela- 
tions which might not otherwise be discovered. 


It may be well to note that our method of handling the 
parameter @ is somewhat different from the procedure 


__ ? We use the term liquid to describe the amorphous phase 
if its density is greater than the critical density, regardless 
of the temperature. This is a consistent usage from our 
point of view since the properties are continuous at the 
critical temperature. 

30. K. Rice, J. Chem. Phys. 7, 136, 883 (1939). 
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used in the earlier work.* In this earlier study, instead of do, 
we used the concept of the average coordination number of 
the atoms in the liquid. This quantity was treated as 
though it were an extra variable of state, determined under 
any fixed set of conditions by the fact that the free energy 
of the system had to be a minimum. An attempt®* was 
made to estimate the effect of this variable on the energy 
and entropy of the system, in order that the result of 
minimizing the free energy might be determined, and these 
attempts have been continued since. However, they were 
not entirely successful, and it became evident that the 
geometry of the system was too complicated to permit 
this program to be carried to completion. The later work 
indicated’ that the approach outlined above might be 
more practicable and successful. 


2. EXPERIMENTAL PROPERTIES OF ARGON 


It is convenient to start with a summary of 
the properties of argon at its normal melting 
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Fic. 1. Equation-of-state data for argon. The curves 
connect points for constant volume; the number at the 
right-hand end of each curve gives the volume in cc/g. 
he number just below the volume gives the value of a. 
Open circles are from Bridgman; black circles are from 
Mathias, Onnes, and Crommelin. The relation between the 
melting curve and the constant-volume curves is such as 
to indicate somewhat different values of the volume of the 
liquid along the melting curve than those given by 
Bridgman. At 2000, 4000, and 6000 kg/cm? we find, 
respectively, 0.654, 0.633, and 0.622 cc/g. 


. 0. K. Rice, J. Chem. Phys. 10, 654 (1942) (in this note 
the number of holes” was used as a variable instead of 
the coordination number). 

° See reference 1, page 17. 
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Fic. 2. (0S/dV)r curves. The open circles are taken or 
interpolated directly from Table V (this table was also 
used to assist in construction of the curves of Fig. 1, as 
noted in text). The black circles are taken from measure- 
ments of slopes of curves in Fig. 1. Each circle in every case 
belongs to the curve nearest it. 


point, 83.78°K, and atmospheric pressure; this 
is given in Table I where sources are indicated. 

The other data which are available consist 
principally of equation-of-state data due to 
Bridgman* and to Mathias, Onnes, and Crom- 
melin.? These are summarized in Fig. 1. The 
open circles are interpolated (in some cases there 
are slight extrapolations) from Bridgman’s data, 
and give the variation of pressure with tem- 
perature for a number of constant volumes, as 
indicated. The points for a given volume are 
connected by curves, whose shapes are in part 
determined by theoretical considerations, as will 
be described below. The black circles are from 
Mathias, Onnes, and Crommelin’s data. Their 
volumes have been corrected slightly so that 
they will be consistent with the specific volumes 
of Baly and Donnan and of Bridgman, which 
agree well with each other (see note to Table V, 
below). 


6 P. W. Bridgman, Proc. Am. Acad. 70, 1 (1935). 

7E. Mathias, H. K. Onnes, and C. A. Crommelin, Proc. 
Roy. Acad. Sci. Amsterdam, 15, 667, 960 (1912); C. A. 
Crommelin, ibid. 13, 54 (1910); 16, 477 (1913). 

8E. C. C. Baly and F. G. Donnan, J. Chem. Soc. 
London 81, 907 (1902). 
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TABLE II. Entropies of solid and liquid argon at various 
temperatures in calories per mole per degree. 











. Stiquid 
T(°K) Ssolid AStus ASexpan (a =4.037) 
83.78 9.07* 3.30 0.00 12.42 
126.3 11.39 2.07 1.28 14.74 
161.9 12.65 1.59 1.78 16.02 
192.9 13.43 1.34 2.01 16.78 








oleae oy ag boy pom rank 'Zelte f \ Bicktrochenie 49. 
308 (1943), which has recently come to my attention; the difference is 
not important. 

Figure 2 gives the value of (@P/dT)y or 
(0S/8V)r, converted to cal. cc! deg.—!, for 
several selected temperatures, as obtained from 
the curves of Fig. 1. Figure 2 will be discussed 
in detail in a later section, but will be found 
useful here in enabling us to determine the 
entropy of the liquid at its normal volume, 
where a=4.037, at various temperatures. For 
this purpose it is necessary to find first the 
entropy of the solid in equilibrium with the 
liquid at various temperatures. This is obtained 
from the values of © found for '=0.18 in Table 
I of the preceding article.’ We then obtain the 
change of entropy on fusion (A.Stus) from the heat 
of fusion as given by Clusius’" at the normal 
melting point of 83.78°K and by Bridgman at 
the higher temperatures. Finally, for the higher 
temperatures, we obtain the entropy of expansion 
(ASexpan) Of the liquid from its volume when in 
equilibrium with the solid to its ‘normal’ 
volume, for which a=4.037, by integrating the 
(0S/8V)r values. The limits of integration are 
shown on the curves in Fig. 2 by short vertical 
lines. For consistency these volumes are adjusted 


RICE 


to conform to the calculated V’s for the solid 
given for [=0.18 in Table I of the preceding 
article, using Bridgman’s values of A V of fusion. 
The results are given in Table I]. 


3. APPLICATION OF THE EQUATIONS TO 
ARGON AT ITS “NORMAL” VOLUME 


In attempting to apply the equations, we come 
at once to a difficulty, since, as noted above, we 
do not have any theoretical way to find dp» as 
a function of a. However, the considerations of 
Section 1 do indicate that a) may be approxi- 
mated by the average distance between nearest 
neighbors. The atomic distribution in argon has 
been studied by Lark-Horovitz and Miller’ and 
by Eisenstein and Gingrich,'* who agree that 
near the freezing point the first maximum in the 
distribution curve of liquid argon occurs at 3.79 
or 3.80A, within about 1 percent. The average 
distance of nearest neighbors is slightly displaced 
toward greater distances, though it is not too 
exactly defined because of the broadness of the 
maximum. We have, therefore, carried out com- 
putations of a number of the thermodynamic 
quantities for a=4.037A and a9=3.82A, 3.85A, 
and 3.88A. These are given in Table II]. In this 
table S and C, are calculated directly from Eq. 
(9) and Eq. (12) from the value of a» given, 
determining 0 from da» by use of the preceding 
article. L and M are evaluated as indicated by 
comparison of Eqs. (4) and (1) using Birge’s 
values of the constants.'* We have evaluated 
da,/da and dE,/da empirically by using two 
experimental values of the pressure (given in 


TABLE III. Calculated thermodynamic properties of liquid argon at a=4.037A, L =1.7584, 
M=0.8769A™ deg.-3. 

















ao 3.82 3.85 3.88 
dao/da 0.487 0.441 0.411 

dE,/da 993 1067 1128 
T 83.78 192.9 320 83.78 126.3 192.9 320 83.78 192.9 320 
S 12.08 16.27 18.92 12.19 14.29 16.51 19.23 12.35 16.80 19.59 
C, 4.88 5.17 5.31 5.05 5.19 5.31 5.43 5.23 5.45 5.54 
sg (0) 58.1 (123.9) (0) 23.3 58.7 (123.9) (0) 59.2 (123.9) 

(0S/AV)r 6.544 0.524 0.514 0.562 0.540 0.523 0.505 0.571 0.523 0.499 








Units: ao and a in A; Ep in cal. mole“; T in °K; S and Coin cal. mole~! deg.~!; P incal. cm~3; V in cm? mole. To convert P to kg/cm? divide by 
0.02343. 


90. K. Rice, J. Chem. Phys. 14, 321 (1946). 


10K. Clusius, Zeits. f. physik. Chemie B31, 459 


(1936). 


4K. Clusius and K. Weigand, Zeits. f. physik. Chemie 


B46, 1 (1940). 


12K. Lark-Horovitz and E. P. Miller, Nature 146, 459 
(1940). 

13 A, Eisenstein and N. S. Gingrich, Phys. Rev. 62, 261 
(1942). 

14 R. T. Birge, Rev. Mod. Phys. 13, 233 (1941). 
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TABLE IV. Calculated thermodynamic properties of liquid argon at a=4.037A, L=1.9, M=0.8769A~ deg.—. 
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ao 3.82 3.85 3.88 

dao/da 0.494 0.453 0.425 

dE,/da 984 1050 1109 
: 83.78 192.9 320 83.78 126.3 192.9 320 83.78 192.9 320 
5 12.38 16.61 19.29 12.51 14.63 16.87 19.61 12.70 17.19 19.99 
A (0) — 580.0 — 1634.2 
C, 4.93 5.20 5.35 5.10 5.24 5.35 5.47 5.28 5.47 5.58 
P (0) 58.0 (123.9) (0) 23.3 58.6 (123.9) (0) 58.9 (123.9) 

0.538 0.524 0.507 0.565 0.522 0.502 


(A.S/dV)r 0.542 0.526 0.514 0.556 





Units: A in cal. mole; for the others see Table III. 


parentheses) in Eq. (13), thus setting up two 
simultaneous equations of the form of Eq. (13). 
Since © and d@/dap are known from the pre- 
ceding article, and all other quantities are also 
known except dao/da and dE,/da (which are the 
same in the two simultaneous equations, since 
they depend only on a, not on 7), we can solve 
for the latter two quantities. Having done this, 
the values of P not enclosed in parentheses are 
obtained from Eq. (13), and the values of 
(0S/0V)r from Eq. (11). 

It will be observed that none of the thermo- 
dynamic quantities depend strongly on the value 
chosen for do. However, the values of S are a 
little low as compared to those given in Table II. 
Since S does show a slight increase with do, this 
discrepancy could be lessened by going to larger 
values of do, but it is felt that a value of a» too 
much larger than that for the maximum in the 
atomic distribution curve should not be used. 

On page 16 of SML the value of LZ was deter- 
mined by certain considerations concerning the 
relaxation of the effect of ‘‘side-atoms”’ on going 
from the solid to the liquid state. If we imagine 
an atom vibrating in a crystal lattice directly 
toward one nearest neighbor and directly away 
from the opposite one, the field of force in which 
it moves is also influenced by the atoms off at 
an angle, these being fixed in their positions in a 
lattice. This means that 0 for the solid is roughly 
2! times as great’ as it would be if these side- 
atoms did not exist. In the liquid it is assumed 
that these side-atoms, not being fixed in position, 
but being free to move, on the average adjust 
themselves about an atom which happens to be 
displaced toward one of its neighbors and away 
from the opposite one. Thus, since we use 9 for 





% Exactly, if only nearest neighbors are effective, and 
the lattice is cubic close-packed. 








the solid in our partition function, we should 
divide 9 by 2?. Actually, however, the side-atoms 
do have space-filling properties, which reduces 
the free space available to the vibrating atoms, 
and the net result is that @ should be divided 
by 2/2? or 2!. Now since there are fewer nearest 
neighbors in the liquid than in the cubic close- 
packed solid argon, the space-filling properties 
of the neighbors may not be quite so great as for 
the solid; furthermore, if an atom in a liquid is 
displaced toward one neighbor, there will not in 
general be another neighbor exactly opposite. 
This suggests that there may be a little greater 
relaxation in the liquid than is given by dividing 
© by 2!, and that, therefore, we should take a 
slightly larger value for Z than that given" in 
Table III. We have, therefore, tried L = 1.9, and 
have obtained the results shown in Table IV. To 
this table we have, for reasons which will appear 
later, added, in one case, values of the molal work 
function A, taken arbitrarily zero at the normal 
melting point. 

It is seen that with d)=3.85 we have achieved 
quite satisfactory agreement with the entropy 
values given in Table II. It may, of course, be 


16 The considerations of SML dealing with the gas of 
hard elastic spheres (before the effect of interatomic forces 
was taken into account at all) indicate the possibility that 
the partition function in Eq. (3) above may be an under- 
estimate from another cause (see p. 13, especially Fig. 5, 
of SML). This would mean that Qo in Eq. (1) would also 
be an underestimate. This would cause an increase in the 
whole expression in brackets in Eq. (1), rather than just 
in the first term, but would have about the same effect. 
On the other hand, the communal entropy was handled 
in such a way as to make the equations give an over- 
estimate for Qo (see Section 1 of SML), so these effects 
counterbalance each other to a greater or less degree. It 
may be remarked, however, that the relaxation of the effect 
of the side-atoms could possibly cause some increase in the 
communal entropy carried by the vibrational term in Qo 
over that characteristic of a Debye solid, especially since 
in the latter the communal entropy is not fully excited (see 
page 7 of SLM). 





330 


O. K. RICE 


TaBLE V. Calculated thermodynamic properties of argon at various volumes, L = 1.9, M=0.8769. 








a 4.15 4.10 4,06 4.037 4.015 3.97 3.94 
ao 3.8912 3.8753 3.8600 3.85 3.8400 3.8195 3.8058 
day/da 0.287 0.351 0.415 0.453 0.456 0.456 0.456 
dE,/da 1057 1070 1065 1050 962 740 547 
Toas* 101.58 94.24 87.75 83.78 100 
P at Tot (0.1) (0.0) (0.0) (0.0) 15.1 
P at 192.9° 36.5 44.9 53.2 58.6 66.2 
P at 320° (83.6) (99.3) (114.4) (123.9) (133.2) (155.4) (173.8) 
(@P/aT vy at Trat 0.416 0.474 0.525 0.556 0.565 
(aP/dT )y at 192.9° 0.384 0.440 0.489 0.524 0.538 
(aP/aT)y at 320° 0.361 0.419 0.473 0.507 0.520 0.544 0.563 
(10.5/12)(dao/da)(dEy,,/dav)t 536 610 645 639 565 360 190 

















* Tsat is the temperature at which the liquid having the given value of a is in equilibrium with its vapor (except for a =4.015, where the lower 
temperature was taken arbitrarily). These temperatures were obtained from the data of Mathias, Onnes, and Crommelin (see reference 7), their 
values being slightly corrected so as to conform to Baly and Donnan’s value (see reference 8) at the melting point, but so as to be consistent with 
their own values of a, as given in Table VI. 





t See Section 7, subsection, “Potential Energy of Liquid Argon.” 


said that this is obtained only by the use of two 
adjustable constants, ad) and L. However, the 
adjustment of L is really quite small, it being 
increased over the originally estimated value by 
only 8 percent, while the value of do used is 
entirely reasonable, and the results are not 
extremely sensitive to it. In any event, we have 
matched three values of S reasonably well with 
the value of a used, and the only slightly ad- 
justed value of L. Furthermore, the curve for P 
against J comes out to be almost a straight line, 
which is obviously correct as indicated in Fig. 1, 
though too much stress on this point must not 
be laid, since presumably any reasonably sensible 
theory would give this result. 


TABLE VI. Specific heats of saturated liquid 
(calories per mole per degree). 








T(°K) 83.78 87.75 94.24 101.58 
a(A) 4.037 4.06 4.10 4.15 
a(deg.—) x 10°* 4.20 4.38 4.74 5.23 
C,(obs. )** 5.50 5.50 5.50 
C,(obs. )** 10.5 10.72 11.04 
C,(obs. )f 10.06 
C,. (calc. } 5.10 5.10 5.08 5.04 
Cp —C,(calc.) 5.49 by gh 6.20 6.73 
C,(calc.) 10.59 10.85 11.28 ER ry 








* Mathias, Onnes, and Crommelin, reference 7 (values at 83.8 and 
87.75 obtained by extrapolation). 

** Eucken and Hauck, reference h, Table I. (Actually the specific 
heat at saturation was measured, but over the range involved this is 
essentially Cp.) 

+ Clusius, reference 10. 


TABLE VII. Entropies and specific heats at a=4.037A 
(calories per mole per degree). 








aU) 83.78 126.3 161.9 192.9 
S (Table II) 12.42 14.74 16.02 16.78 
S (Table IV) | 14.63 15.94 16.87 
Cv, w. (Table IT) 5.63 5.1 4.3 

Cv, ww. (Table IV) 5.17 5.3 5.3 








We shall return to a discussion of the specific 
heat values in Section 5. 


4. APPLICATION OF THE EQUATIONS TO 
ARGON AT VARIOUS VOLUMES 


By combining the results of Bridgman and 
Mathias, Onnes, and Crommelin, as presented in 
Fig. 1, we have a series of data at high temper- 
atures and at one much lower temperature for 
each of a number, of volumes greater than the 
“normal” volume of melting argon; i.e., values 
of a greater than 4.037A. These data can be used 
in the same way as the data at a=4.037A to 
obtain values of dao/da and dE,/da, as well as 
the various thermodynamic quantities. Starting 
with a@=3.85A at a=4.037A, ao at any other 
value of a will be determined by integration of 
day/da; on the other hand, dao/da at any given 
a itself depends upon the corresponding value of 
dy. It is, therefore, obvious that it is necessary 
to make successive adjustments to assure con- 
sistency in this matter. This, however, is rela- 
tively easy, because the dependence of dao/da 
on @» at a given a is slight. 

For values of a less than 4.015, we see from 
Fig. 1 that the slopes of the P vs. T curves change 
sharply at the lower temperatures. On the basis 
of Eq. (13) this is a distinctly ‘‘abnormal’’ be- 
havior. This “abnormal’’ behavior will be dis- 
cussed later ; we believe that the high temperature 
ends of these curves are ‘‘normal.’’ Unfortu- 
nately, there is in each of these cases only one 
experimental point at a high temperature, so the 
slope in this region can only be inferred. The 
curves for a<4.015 shown in Fig. 1 are drawn 


,*, — ss 





and 
d in 
per- 
- for 
the 
lues 
used 
\ to 
ll as 
‘ting 
ther 
n of 
iven 
1e of 
sary 
con- 
rela- 
1o/da 


from 
ange 
basis 
* be- 
. dis- 
ature 
ortu- 
7 one 
o the 
The 


rawn 





with slopes (for the high temperature ends of 
the curves) calculated from Eq. (11) on the as- 
sumption that as a goes below 4.037 the value 
of dao/da reaches a limit of 0.456. It is obvious 
that the slopes shown in Fig. 1 cannot be greatly 
in error, but some more experimental work in 
this region would be most welcome. 

The calculations mentioned in the above two 
paragraphs are summarized in Table V. Again 
pressures which are obtained directly or by 
interpolation from the experimental data are 
given in parentheses. It was these which were 
used for a >4.037 to calculate day/da and dE,/da. 
(For a=4.015 it was assumed, as noted above, 
that day/da=0.456 and P at 100° was then cal- 
culated). In all cases the values of 0 necessary 
for the calculations were obtained from do, using 
Eqs. (1) and (2) of the preceding article, taking 
r=0.18. 

The curves given in Fig. 1 were all drawn with 
the aid of the values of P and (dP/dT)y given 
in Table V (except for the low temperature ends 
of those with a <4.015). Since the values of a in 
Table V and Fig. 1 do not all coincide, inter- 
polation, or, in some cases, extrapolation, was 
necessary in order to construct the curves. 


5. THE SPECIFIC HEATS 


The data on the specific heat of the liquid are 
summarized in Tables VI and VII. The data in 
Table VI refer to liquid in equilibrium with its 
vapor. C,(calc.) comes directly from Eqs. (5) and 
(12) using Table V. C,—C,(calc.) is obtained 
from the formula 


C,—C,=aVT(dP/dT)y, (14) 


a is obtained from the density measurements of 
Mathias, Onnes, and Crommelin? (these all being 
made at what amounts essentially to zero pres- 
sure). (0P/dT)y =(8S/dV)r is taken from Table 
V, or, what is the same thing, from Fig. 1. 

In Table VII the average C,’s were calculated 
from the changes in entropy with temperature; 
the values from Table II are essentially the 
“experimental” values, though some calculation 
is involved in them, also, since they depend on 
the entropy of the solid, which has not been 
measured directly. To the values of Table IV 
one at 161.9° has been added. 

It will be seen from Table VI that our cal- 
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culated values of C, are somewhat lower than the 
experimental values of Eucken and Hauck. This 
discrepancy seems to be confirmed by Table 
VII in view of the difference between the experi- 
mental and calculated C,’s for the range between 
83.78° and 126.3°. On the other hand if the cal- 
culated value of C, at 83.8°, given in Table VI, 
is too small, then the calculated C,—C, must be 
too large; otherwise, C, would considerably 
exceed the experimental value. However, C,—C, 
is actually pretty well tied to the experiments. 
The value of (@P/dT)y used in Eq. (14), from 
which C,—C, is calculated, is one which fits 
Bridgman’s data as shown in Fig. 1; in fact, the 
only role of the theory here is to indicate the 
variation of (0P/dT)y with the temperature. The 
data do not definitely exclude a slightly different 
curvature of the P vs. T curve. (@P/dT)y might 
have a lower value than that given in Table V, 
at the lower end of the curve, and a higher value 
at the upper end; the theory may be wrong in 
this respect, though it seems unlikely.'? Further- 
more, the independent direct measurement by 
Simon and Kippert of (¢P/dT)y for the normal 
liquid gives a value greater than that inferred 
from Bridgman’s data (see Table 1). Also, in 
applying Eq. (14) we have used the smaller of 
the two values of a, a value which was in fact 
obtained by extrapolation. In other words, our 
estimate of C,—C, appears to be about the least 
possible. It, therefore, seems justifiable to say 
that there are inconsistencies between various 
sets of experimental data as great as any errors 
in the theory. We may note further that, while 
we have no direct check on the C, values of 
Eucken and Hauck, disagreements appear in 
other thermal measurements. Thus the C, values 
of these authors do not agree exactly with 


17 A slight decrease in (@P/dT)y, that is (0S/8V)r, with 
increasing temperature seems normal from general con- 
siderations. In contrast to a liquid composed of hard 
elastic spheres, which would, at any given total volume, 
have a free volume independent of temperature, the free 
volume of a real liquid increases with temperature. Hence 
a given increase in total volume should have a smaller 
relative effect on the free volume at a higher temperature, 
so that (0.8/8V)r should be smaller. 

We may look at this matter in another way. We have 


[a(aS/dV)r/dT lv = [0(0S/0T )y /dV Jr =T-(8C,/dV)r. 


But as the liquid expands and becomes more gas-like, C, 
should decrease, and (0C,/dV)r should be negative, which 
means, then, that (0S/0V)7r decreases with increasing tem- 
perature. 
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Clusius. Clusius has noted that his C,’s for solid 
argon are 1 to 6 percent lower than earlier values 
of Eucken, while his value of the heat of fusion is 
considerably higher than that of Eucken. 

Turning to Table VII, we see that the cal- 
culated trend of C, values is upward with in- 
creasing temperature, while the experimental 
values appear to decrease. Here, again, the theory 
reflects the normal expectation. A slight rise of 
C, with temperature is to be expected because 
the vibratory motion, as compared to the trans- 
lational, should become relatively more promi- 
nent at the higher temperatures where its 
amplitude increases. 

We note that if the entropies at 126.3° and 
161.9° from Table 11 were lower by about 0.2 
entropy units, the discrepancy in the C, values 
would be removed. An error of this magnitude 
is perhaps not beyond the limits of experimental 
_ inaccuracy, since it would be only 6 percent of 
the entropy change of about 3.35 units which 
occurs in going from melting solid to liquid with 
a=4.037A. 


Since the relation (0S/8V)7=(0P/dT)y holds, whether 
the system is undergoing a change of phase or not, the 
discrepancy would be removed if Bridgman’s dP/dT for 
the melting curve, and the average ordinates of the curves 
for 126.3° and 161.9° in Fig. 2 above, were all too high by 
6 percent. The curves of Fig. 2 certainly may well be in 
error by this amount. And there is some suggestion from 
the more recent though less extended melting curve of 
Clusius and Weigand" that Bridgman’s melting curve may 
be this much in error. At 0° and the normal melting point, 
Clusius and Weigand found dP/dT =40.5 kg/cm™ deg.“, 
while Bridgman’s value is 42.0. (The discrepancy is not so 
great as indicated by Clusius and Weigand, who must 
have made an incorrect conversion of units.) However, 
the curvature of Clusius and Weigand’s P vs. T curve is 
very much greater than that of Bridgman’s, so that already 
at about P = 200 kg/cm?, the limit of Clusius and Weigand’s 
pressure range, the two slopes become equal. Assuming 
the careful measurements of Clusius and Weigand to be 
correct over their limited range of pressures, the question 
arises as to whether, with rising pressure, the slope is 
correctly given by Bridgman’s data, or whether it con- 
tinues to increase. If the latter occurs, the true pressures, 
which (assuming them to be given by Clusius and 
Weigand’s data) at temperatures very close to the normal 
melting point lie below Bridgman’s curve, would soon 
cross it. If the true pressures again crossed Bridgman’s 
curve around 135°K, there would be a region where the 
slope was a little less than that given by Bridgman. This 
could happen without any of Bridgman’s temperatures 
being in error by more than about +1 to 1.5°, but there 
would be a slight inflection in the P vs. T curve. Such an 
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inflection would not be surprising if it were associated with 
the prefreezing phenomenon mentioned later in the paper. 
However, as noted in Section 6, there is a small discrepancy 
between the free energies of the liquid and solid along the 
melting curve, which would be somewhat increased if new 
experiments should confirm changes such as are suggested. 
Furthermore, it is difficult to make a complete analysis of 
the situation because of the way various experimental 
errors may interact with each other.'* So the above analysis 
is not to be taken too literally, but it does suggest that 
the calculated values in Table VII are correct to within 
experimental error. Complete consistency between all 
possible data analyzed on the basis of our theory, is 
perhaps too much to hope for, in view of our rather bold 
use of the Debye theory of the solid, when there are strong 
indications” that it may need some modification on account 
of anharmonicities in the vibrations of the atoms. 


The calculated value of C,—C,, as presented 
in Table VI, rises much more rapidly with tem- 
perature than the experimental values of Eucken 
and Hauck. We again recall that the calculated 
value of C,—C, is closely tied to experimental 
data, but it is again true that the data of Fig. 1 
do not exclude somewhat smaller values of 
(0P/dT)y at the lower end of the curves. Any 
change in curvature of the P vs. T curves suff- 
cient to make (0P/dT)y increase with increasing 
temperature has the same objections as those 
previously noted.'? Even if we assume that 
(@P/dT)y is independent of temperature at 
a=4.15, say, C,(calc.) for 101.58° is reduced 
only to 11.20 cal. mole deg.—', and this value 
is based on the calculated C, of 5.04, which is 
much lower than Eucken and Hauck’s experi- 
mental value. I believe, however, that while the 
decrease in (0P/dT)y may not be as great as 
given in Table V, (0P/dT)y certainly will de- 
crease with temperature to some extent. This, 
therefore, is further indication of error in some, 
at least, of the experimental results. The work of 
three sets of investigators enters into the com- 
parison of values of C, in Table VI, namely, that 
of Bridgman, that of Mathias, Onnes, and 
Crommelin, and that of Eucken and Hauck. 
While it seems unlikely that any one of these 

18An error in the melting curve would not itself be 
without other complicating effects. A change would have 
to be made in the calculations of Table I of the preceding 
paper. These changes would be rather small, and the effect 
on the estimated entropy of the liquid (aside from that 
caused directly by changes in the slopes of the curves) 
could be neglected; but other changes in the experimental 
results might well result in unforeseen complications. 


190. K. Rice, J. Chem. Phys. 12, 289 (1944); and see 
the preceding paper (reference 9). 
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pieces of work is sufficiently in error to account 
for the discrepancies involved, it may well be 
that the combined errors in all of them are great 
enough. 

That this is a possibility is also suggested by 
some experiments on the velocity of sound in 
liquid argon made by Liepmann.” A slight extra- 
polation of his results gives a velocity of 8.30 X 104 
cm/sec. at 87.75°K. The velocity of sound is 
given in terms of thermodynamic quantities by 
(C,—C»)*#(Cp/C.a2T)' (see Liepmann’s Eq. (5)). 
Using the values of a from Table VI and our 
calculated values of the specific heats (converted 
to ergs per gram) we obtain 8.73 X10‘ cm/sec., 
about 5 percent too high. If we use, instead, the 
specific heats of Eucken and Hauck, the velocity 
is calculated at 7.73104 cm/sec., about 7 per- 
cent too low. 

In whatever way the questions raised in this 
section are ultimately settled, they do indicate 
one thing which is most important for the 
present; namely, there is a great need of new, 
highly accurate experimental work both on the 
thermal properties and the equation of state. 
The data now at hand represent excellent and 
beautiful pioneer work; but it is pioneer work 
and does not give a sufficiently detailed account 
of liquid argon to furnish a satisfactory basis for 
the theory. New measurements, using the best 
modern techniques, are urgently required. 

In the meantime, since the theory seems to 
represent the various sets of experimental results 
within approximately the limits indicated by 
their mutual inconsistencies, it seems worth 


TABLE VIII. Free energy changes along the melting curve 
(calories per mole.) 








T(°K) AA ts PmVs AEp AF, AAy, l PyV I AFi 


126.3 —432 1145 —5 708 —580 652 639 711 
192.9 —1190 3374 -40 2144  -—1634 1642 2145 2153 











“°H. W. Liepmann, Helv. Phys. Acta 12, 421 (1939). 
Liepmann has noted the inconsistencies between his values 
and the thermal data. The inconsistency between his 
results and Eucken and Hauck’s specific heat data would 
not appear to be quite so great had he used the values of 
a given in Table VI instead of the values given by Donnan 
and Baly; however they would still exist. It seems to me 
unlikely that the possible explanation offered by Liep- 
mann, namely, that the sound wave is too fast to allow 
internal equilibrium to be established in the liquid, can 

correct, since @=80 corresponds to a frequency of 
about 1.710", far higher than the frequency of the 
sound wave used by him. 





5) 






































| 
|| 

40 4 

a (Angstrom units) 


Fic. 3. 





while to use it to make some calculations and 
deductions. 


6. THE FREE ENERGY ALONG THE 
MELTING CURVE 


Let us write AF, for the free energy of the 
solid at any temperature and pressure, minus its 
free energy at zero pressure and its normal 
melting point at zero pressure; and let us write 
AF, for the similar quantity for the liquid. Then 
at any point along the melting curve we should 
have 

AF,=AF). (15) 
We write 
AF,=AA;, s+ PmnVst+AE>. (16) 


Here AA;,, is the change in A,,, the thermal 
part of the work function of the solid, along the 
melting curve; it is readily computed, using the 
Debye theory, from the 0’s given in the pre- 
ceding paper. P,, is the pressure along the melting 
curve, at the temperature of interest, and is 
known from Bridgman’s work. V, is the volume 
of the solid at the same temperature obtained 
from the V(calc.) of Table I of the preceding 
paper. The change in potential energy AE, can 
be readily calculated, since E, is known as a 
function of @ from Eqs. (7), (1), and (2) of the 
preceding paper. As always, we use I =0.18. 
For AF; we write 
Pn 


AF,=AAy,i+PyV+ VdP. (17) 


Py 


Here AAy,; is the change in work function at 
constant volume (a =4.037) (the A of Table IV). 
Py is the pressure reached at the temperature of 
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interest at the constant volume V (with 
a=4,037) ; the values given in Table IV are used. 
The integral takes care of the change in F when 
the liquid is now compressed at constant tem- 
perature until solid begins to form, and can be 
evaluated from Bridgman’s equation-of-state 
data (Fig. 1). 

AF, and AF; and the other quantities in Eqs. 
(16) and (17) are given in Table VIII. The 
integral in AF; is designated as J. Equation (15) 
is reasonably well satisfied. 










We have noted that AAy,: is calculated from the theo- 
retical results of Table IV. Since some of the entropies in 
this table are higher than the experimental values and 
some lower, the average which goes into the equation 
(8A /dT)y=—S is about the same as if the experimental 
values had been used. However, if the experimental values 
in the middle temperature range were adjusted downward, 
as proposed in the preceding section, the discrepancy in 
Table VIII would be increased by use of the experimental 
values. This suggests that possibly the value of S at the 
normal melting point should be adjusted upward. This 
would have to represent a change in the AS of fusion, 
because if the entropy of the solid were also changed no 
net effect would result on the difference between AF, and 
AF, since solid and liquid would be affected similarly. It 
seems unlikely, however, that there is any appreciable 
error in Clusius’ value of the heat of fusion, so the normal 
AS of fusion should be well known. Although the difference 
between AF; and AF, is small, it seems more difficult to 
account for the discrepancy as experimental error than it 
does to account for the discrepancies involving the specific 
heats discussed in Section 5. The difficulty may be con- 
nected with the possible breakdown of the Debye theory 
of the solid, discussed in previous work,!* and mentioned in 
Section 5. At the same time, it may be emphasized again 
that a complete reworking of the data on argon would be 
of the utmost value. 


























7. DISCUSSION 


Behavior of da,/da 


In Fig. 3, we have plotted dao/da against a. 
It will be noticed that there is a sharp break in 
the slope of this curve between a=4.03A and 
a=4.04A. This break corresponds to the break in 
(0S/8V)r for 320°K, coming at the same place 
as shown in Fig. 2. Though we cannot exclude 
the possibility that this break in the slope of the 
curves should be rounded off somewhat, there 
seems to be no doubt that the break in (0.S/0V) rp 
is real in the sense that there is a fairly marked 
change of slope in the region shown; the inter- 
pretation of it in terms of da)/da depends upon 














RICE 





our acceptance of the behavior of © with a, 
described in the preceding article.” 

It might be possible to suppose that the break 
in the slope of the (0S/0V)r curve results from 
a change in the nature of the dependence of 
d@/day, rather than of dao/da on ao. Especially 
might we consider this to be a possibility if the 
parameters were chosen (as is not excluded) in 
such a way that the value of a@ involved was 
outside the experimental range of lattice dis- 
tances in the solid where we have direct knowl- 
edge of the behavior of 0. This would mean that 
we refer the break in slope of (0S/0V)r to a 
property of the solid rather than to a special 
property of the liquid. Aside from the inherent 
improbability of such an assumption, I believe 
that we can find some circumstantial evidence 
that the break is a distinct property of the liquid. 

This rests on an examination” of the atomic 
distribution curves found by Eisenstein and 
Gingrich,'* and by Lark-Horovitz and Miller.” 
Eisenstein and Gingrich found curves for liquid 
argon at 84.35°K and 91.8°K which differed very 
noticeably from each other. Each of these curves 
showed a maximum just below 3.8A, but there 
were 10.2 to 10.9 atoms under the first peak at 
84.35°K and only 6.8 to 7.2 at 91.8°K, and the 
maximum was much lower in the latter case. 
Furthermore, at 91.8° there was a second peak 
at 4.7A, but at 84.35° there was a rather sharp 
minimum near this distance. Lark-Horovitz and 
Miller found a distribution curve at 89.2°K 
which appeared to be much more like the dis- 
tribution curve at 84.35°K than that at 91.8°, 






which the preceding paper modifies, led to a very radical 
variation of d@/day with ao. Applied to the liquid, this 
made dao/da increase rapidly with a, so that dao/da 
became greater than 1 when the liquid was expanded 
somewhat beyond a=4.037A. This seemed very unreason- 
able because it meant that the compressibility of the 
individual atoms was greater than that of the liquid as a 
whole. At large values of a, we may, in fact, expect do to 
approach a constant value, determined by the “size” of 
a single atom; da/dao, then, should decrease asymptotically 
to zero. It was these considerations which led us to make 
the recalculation for the solid.- 

A still earlier paper (O. K. Rice, J. Am. Chem. Soc. 63, 
3 (1941)) gave d In @/dV=0.10 instead of the 0.18 used 
here. The low value of dln @/dV gives a value of dao/da 
equal to about 1.9. This is obviously ridiculous, and may 
be added to the reasons already discussed for discarding 
these calculations. 


* This examination was greatly facilitated by more- 


detailed data very kindly sent me both by Dr. Gingrich 


and Dr. Lark-Horovitz. I wish to thank these investigators 
for their assistance in this matter. 








21 The calculations of the earlier article (see reference 19), ° 
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and which had about 10.0 atoms under the first 
peak. One thus gets the impression that there is 
a sharp change between 89.2° and 91.8°, amount- 
ing almost to a transition. It seems very natural 
to suppose that this change in the atomic dis- 
tribution would be associated with a change in 
the behavior of day/da with a. But if the transi- 
tion takes place between 89.2° and 91.8° one 
would expect the break in the curve in Fig. 2 
or Fig. 3 to take place at a=4.07 or 4.08 instead 
of around a=4.04. However, a closer examina- 
tion shows that the change in the atomic dis- 
tribution curve is actually more gradual than 
appears on the surface. A close scrutiny of 
Eisenstein and Gingrich’s Fig. 10, and Lark- 
Horovitz and Miller’s Fig. 2, allows us to make 
the following estimates of the density function 
4rr’p(r) at the maximum (that is, the height of 
the maximum) at the various temperatures: 
84.35°K, 12.7; 89.2°K, 10.5; 91.8°K, 9.3. Also, 
we make the following estimate of the value of 
this function at the first minimum (first two 
temperatures), or the second maximum (last 
temperature) : 84.35°K, 2.5; 89.2°K, 3.9; 91.8°K, 
5.9. Thus, certainly within the limits of error, we 
can say that we have a gradual decrease of 
density at the first maximum, and at the next 
position, around 4.7A, a gradual increase in 
density. Just how suddenly the minimum appears 
to change into a maximum only a more closely 
spaced and very accurate series of experiments 
by a single investigator could show. It may be 
that the second peak, which at 84.35°K is around 
5.3A, slides with increasing temperature into the 
position at 4.7A. In any event, at present the 
process looks like a gradual shift from a dis- 
tribution resembling a distorted face-centered 
cubic lattice to one resembling a body-centered 
cubic lattice. This goes on all through the tem- 
perature range from near the normal freezing 
point of argon to 91.8° and beyond (from the 
work of Eisenstein and Gingrich it appears to 
continue on into the high temperature range), 
which means that the change starts when a is in 
the neighborhood of 4.03 or 4.04A. (It could 
start at a slightly higher, but presumably does 
not start at a much lower value of a, for the dis- 
tribution at 84.35° resembles a face-centered 
cubic distribution about as closely as would 
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seem to be possible without having a definite 
lattice of that form). The apparently sudden 
change in coordination number between 89.2° 
and 91.8° is a psychological rather than a 
physical phenomenon. As soon as the minimum 
turns into a maximum, one divides the space 
around the central atom in a different way and 
relegates some of what were ‘‘nearest neighbors”’ 
to the class of “‘next-nearest neighbors.” 

If this interpretation of the atomic distribution 
curves is correct, it seems natural to associate the 
change in the slope of the curve da o/da vs. a 
between a=4.03 and a=4.04 with the beginning 
of the appearance of traces of the distorted body- 
centered lattice at about the same value of a. In 
fact, off hand, one would not be surprised to see 
a discontinuity in da)/da, but, according to the 
present picture, such a discontinuity would be 
associated with a discontinuity of (0S/d0V)r or 
(0P/dT)y, which would then occur at a definite 
volume independent of temperature, and this 
would be impossible by the theory of second- 
order transitions. 


The Potential Energy of Liquid Argon 


In Table V we have tabulated calculated values 
of dE,/da. These can be correlated with the 
potential curve of the solid, as given in the pre- 
ceding article. Since we believe that do is to be 
identified fairly closely with the actual average 
interatomic distance of nearest neighbors in the 
liquid, we expect E, to be related to the E, for 
the solid at the given value of a». In fact, we 
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should be able to write, at least roughly, 
E,p=(N/12)Ep, (a0) + Ea, (18) 


where E, is a constant, E,, ,(@o) is the potential 
energy of the solid at lattice distance a), and N 
is the average coordination number in the 
liquid.** If there is a ‘‘mixture’’ of face-centered 
and body-centered lattices in the liquid, we 
should perhaps add another adjective, and say 
that N is the effective average coordination 
number, as the next nearest neighbors in a body- 
centered lattice are close enough to have quite 
an influence on the central atom. If we differen- 
tiate Eq. (18) we obtain 


dE,/da = 7?3(dN/da)E,, .(ao) 
+(N/12)(dE), ./das)(day/da). (19) 


We may suppose that the greater portion of the 
variation of dE,/da is in the last term of Eq. (19). 
Also the variation of N is probably not too large. 
Therefore, in the last term of Eq. (19) we have 
assumed JN constant at 10.5. In Fig. 4 we have 
plotted dE,/da from Table V and (10.5/12) 
X (dE,, ,/dao)(day/da), taking dE,,,/day) from 
the preceding article and day/da from Table V. 
It is seen that the two curves parallel each other 
in the region of small values of a, which may, at 
least for this region of a, be taken as confirmation 
of our point of view that information about the 
energy of the liquid may be obtained in a simple 
manner from the solid. It must be granted that 
we can only make this correlation because the 
somewhat arbitrarily chosen parameter dy lies 
within the range for which we have experimental 
data on the solid. In spite of.this, we believe that 
the result must have at least qualitative sig- 
nificance. 

The difference at small values of a between the 
two curves in Fig. 4, about 400 cal. mole! A-, 
gives the first term on the right-hand side of 
Eq. (19) and, since E,,,(@o) is near 2000 cal. 
mole, we find dN/da=2.4A~-'. This rate of 
change of N is small enough to justify our taking 


N constant*4 in the last term of Eq. (19); it: 


may seem a little smaller than expected, but 
it is insensitive to the parameters. In any 


*8 Except for the constant Eg, which will be discussed 
later, this is essentially the assumption we have made 
previously (see reference 3). 

*4 A correction would decrease the parallelism of the 
curves, but only slightly. 
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event, it can have a really objective significance 
in terms of the pictorial theory of the liquid 
only below @ equal to 4.03 or 4.04A, where 
the lattice change begins to come in. For a 
greater than about 4.03A, it is seen that dE,/da 
appears to be approximately constant, more so 
than (N/12)(dE,, ,/dao)(dao/da). While we can- 
not ignore that the determination of dE,, ./day 
involves an extrapolation in this region, the 
greater degree of constancy of dE,/da may be 
significant and explainable in the following way. 
Since the position of the first maximum of the 
distribution curve remains almost fixed,'* the 
‘percentage of body-centered structure,’ which 
must thus determine E,, would be expected to 
be a linear function of a, at least over a small 
range. It is then not surprising if E, is similar in 
this respect, and dE,/da consequently nearly 
constant. 

We now return to a discussion of the constant, 
or approximately constant, quantity Ey. The 
long range or translational type of motion, dis- 
cussed in Section 1 above, in which the atoms 
slip around and exchange places, could not take 
place without changes in the distances between 
neighboring atoms. Thus the potential energy of 
the system cannot be the same as it would be if 
there were no such translational motion, and all 
neighboring atoms were at the same distance 
from each other. Even the translational motion 
will involve exchange of potential energy of 
certain pairs of molecules with other pairs, 
though the total potential energy in the entire 
large assemblage must always remain the same 
for any given a, and independent of temperature, 
if the motion is, in fact, to be properly described 
as translational. The energy Ey is the extra 
potential energy associated with the fluctuation 
in interatomic distances necessary for the trans- 
lational motion to occur. Eq is always positive 
because the atoms are near a minimum of the 
mutual potential energy curve, i.e., in a region 
where it is concave upward. 

From the available thermochemical energy we 
can estimate the total energy E of argon liquid 
at its freezing point. (This is reckoned on the 
basis of perfect gas at absolute zero having zero 
energy.) We can estimate E; from Eg. (7). The 
difference E—E,; is E,; since we are dealing with 
the classical case we do not need to consider the 
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zero-point energy. Ea can then be estimated from 
Eq. (18). We find for argon liquid at its normal 
freezing point: 
E=-—1177 cal./mole, 
E,.=398 cal./mole, 
E,= — 1575 cal./mole, 
Ea= 167 cal./mole, 


(assuming N=10.5, a9=3.85, and L=1.9). It 
will be recalled that E, itself contains some 
potential energy. The kinetic energy is equal to 
8RT or 250 cal./mole, so the potential energy 
contributes 148 cal./mole to E;, a quantity of 
the same order of magnitude as Ey. 


Prefreezing Phenomenon 


In Section 4 we have noted that the curves 
in Fig. 1 corresponding to low molal volumes 
have an “abnormal” slope at the low tempera- 
ture end. This is reflected in Fig. 2, where it is 
seen that (0S/0V)7 is ‘‘abnormally” great at 
low temperatures and small values of a. An 
indication of the reason for this apparent ab- 
normality is found when we plot 


(0@E/8V)r=T(dS/dV)r—P 


against a, as in Fig. 5. The curve for dE,/dV is 
included for comparison, and it is seen that at 
320°K the (@E/dV)r curve parallels it reasonably 
well, but at lower temperatures (0E/0V)7 be- 
comes very great at small values of a. The values 
of (#E/dV)r are subject to considerable uncer- 
tainty, because they are calculated as a relatively 
small difference, but the qualitative trend is 
unmistakeable. 

The trend of the curves in Figs. 2 and 5 can 
be explained if we suppose that at the lower tem- 
peratures the system tends on compression to go 
to some special state of low entropy and, simul- 
taneously, low energy. This leads one to specu- 
late that we are dealing here with a prefreezing 
phenomenon, in which submicroscopic particles 
or crystalloids, having essentially the crystal 
structure of the solid, are formed in the liquid. 

Such a phenomenon can, of course, occur only 
if the free energy of the entire system is thereby 
lowered. One might wonder why, if formation of 
what are essentially small particles of solid lowers 
the free energy, the system does not go the whole 
way, and form a large crystal. The answer is, of 
course, the entropy of mixing of the solid par- 
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ticles with the molecules of liquid. When the 
liquid (without crystalloids) and the solid have 
very closely similar free energies, this may be 
sufficient to turn the balance and stabilize the 
mixture. One would not expect such a system to 
be stable over a very extended range. In Figs. 1 
and 2, curves are shown® which indicate the 
points at which solid and liquid are in equi- 
librium. It will be seen that the abnormal values 
of (0S/dV)r actually occur to a considerable 
extent in the metastable region, where the solid 
has a lower free energy than the liquid. Only 
over a rather small range is our hypothetical 
mixture important and stable. It is hoped that 
it may be possible in subsequent work to give a 
more quantitative discussion of these phe- 
nomena. In our earlier discussion of the gas of 
hard elastic spheres we noted the possibility that 
such a gas on compression would form regions 
in which the “‘atoms’”’ are in a cubic-close-packed 
structure (see SLM, p. 14). In this case, however, 
the “atoms” do not have an attractive force 
field, and, furthermore, these groupings are 
associated with an increase in the entropy of the 


% In Fig. 2 the curve is merely indicated by the short 
vertical lines. 
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system. It appears, therefore, that this is not 
the same phenomenon as the “‘prefreezing’’ phe- 
nomenon we are discussing here. We may specu- 
late that both types of grouping can occur in an 
actual liquid, the prefreezing phenomenon in- 
volving larger groups and being characterized 
also by a change in the interatomic distance to 
that of the solid, which is probably considerably 
larger than the average distance between closest 
neighbors in the liquid in the range in which 
prefreezing occurs. 


The phenomena that we have here described as pre- 
freezing were already noted in our earlier work.? The 
explanation given there, namely, that the decrease in 
entropy as the volume is decreased is associated with an 
increase in the average coordination number of the atoms 
of the liquid toward the coordinate number characteristic 
of solid argon, has some resemblances to and some points 
of difference from the explanation offered here. The process 
of prefreezing described in the above paragraphs will 
obviously be associated with an increase in average 
coordination number, but of a somewhat more specialized 
type than that envisaged in the earlier work. In particular, 
in the present discussion we conceive of the coordination 
number as varying appreciably in different regions which 
may be of appreciable size, whereas in the earlier work it 
was supposed that the average coordination number 
changed gradually with no appreciable local variations 
extending over considerable regions. 


Conclusion 


In concluding we may point out that perhaps 
the greatest advance recorded in this paper is 
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the rather good correlation of the properties of 
the liquid with those of the solid. The attempt 
to use values of © taken directly from the solid 
appears to be successful. It is particularly inter- 
esting that it is absolutely essential to take into 
account the relaxation of the forces of the side- 
atoms in the liquid, and, in fact, as we have noted 
in Section 3, the liquid appears to be slightly 
more relaxed than our first estimate indicated. 
The correlation with the solid appears not only 
in the fact that the values of © for the solid give 
a good description of the properties of the liquid 
when proper account is taken of the interatomic 
distance, but also in the examination of the 
potential energy of the liquid and its comparison 
with that of the solid. Considerable uncertainty, 
unfortunately, still remains in the interpretation 
of the properties of the solid, since it does not 
appear to be possible to correlate 9 theoretically 
with the lattice distance. In the liquid we have 
a somewhat similar difficulty, since it requires 
recourse to the experimental data to correlate a» 
with a. Further progress in the understanding 
of the liquid will require solution of both prob- 
lems. Since only one of them is actually charac- 
teristic of the liquid, the other belonging to the 
domain of the theory of the solid state, it now 
seems fair to say that our present understanding 
of the liquid has advanced to a point not too far 
behind our understanding of the solid. 
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The pressure dependence of accommodation coefficients of gases on metals has been explained 
by assuming that the accommodation coefficient varies linearly with the fraction of the surface 
covered with adsorbed gas. It is further assumed that the accommodation coefficient has a 
negligibly small value ap on a gas-free surface and an asymptotic value a. on a saturated 
surface. The assumptions lead to an accommodation coefficient isotherm which reproduces 
the pressure dependence of 119 accommodation coefficient values for ten gases on platinum 
with an average absolute deviation of 1.5 percent. 





T has long been recognized'~* that the accom- 

modation coefficient of a gas on a metal 
surface depends upon the specific nature of the 
surface. Recently, values have been reported® 
for room temperature accommodation coefficients 
on platinum of ten gases over a wide range of 
pressure, so that sufficient experimental results 
are now available to test any theory which 
relates the value of the accommodation coef- 
ficient to characteristics of the surface which 
depend upon the gas pressure. It is the purpose 
of this paper to show that the observed variation 
of the accommodation coefficients with pressure 
can be explained quantitatively by considering 
the effect of adsorption on the metal surface. 
Although it frequently has been suggested'—* 
that adsorption is responsible for the increase in 
the value of the accommodation coefficient in 
changing from a gas-free to a gas-covered surface, 
no general functional relation has been proposed 
to date. 

If we assume that the accommodation coef- 
ficient varies linearly with the fraction of the 
surface covered, the following relation holds: 


a=(1—O)aot+ Oa., (1) 


1(a) J. K. Roberts, Proc. Roy. Soc. A129, 146 (1930); 
(b) A135, 192 (1932); (c) A142, 519 (1933). 

2'W. B. Mann, Proc. Roy. Soc. A146, 776 (1934). 

*'W. B. Mann and W. C. Newell, Proc. Roy. Soc. A158, 
397 (1937). 

*(a) W. C. Michels, Phys. Rev. 40, 472 (1932); (b) 52, 
1067 (1937). 

°B. Raines, Phys. Rev. 56, 691 (1939). 

° P. Rolf, Phys. Rev. 65, 185 (1944). 

7L. B. Thomas and F. Olmer, J. Am. Chem. Soc. 65, 
1036 (1943). 

*I. Amdur, M. M. Jones, and H. Pearlman, J. Chem. 
Phys. 12, 159 (1944). 


where apo, a, and a, are the accommodation coef- 
ficients at zero pressure, at the pressure at which 
the fraction © of the surface is covered, and at 
infinite pressure, respectively. Since we are 
dealing with gases adsorbed on relatively plane 
metallic surfaces at very low pressures and at 
temperatures far above the critical, we shall 
assume that the fraction of the surface covered 
is given by the simple Langmuir relation 


O=bp/(1+5p), (2) 


where p is the gas pressure and d is a constant 
for a given gas and depends only upon the tem- 
perature. Equations (1) and (2) may be combined 
to give the following relation for the pressure 
dependence of the accommodation coefficient: 
b 
a a ee = : ago+ op Onn. 
1+bp 1+ bp 1+bp 

It will be assumed, for the present, that the 
magnitude of ap is sufficiently small so that it 
may be neglected in Eq. (3). This assumption 
is in accord with the available room temperature 
information shown in Table I which compares ap 
and a, for several gases and metals. Although the 





(3) 


TABLE I. Comparison of room temperature values 
of a and aa. ae 








Gas Metal 





Platinum 
Platinum 
Nickel 

Tungsten 
Tungsten 


Helium* 
Helium» 
Helium® 
Helium4 
Neon?4 








8 See reference 2 of text. 
b See reference 6 of text. 
¢ See reference 5 of text. 
4 See reference 4b of text. 
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TABLE II. Test of general relation a=[bp/(1+6)) Jaw. 


























Aa Aa | 

100 100|—— 

Gas @> 6(mm) P (mm) Gobs Qale Gobs Gas a. b(mm) P (mm) Gobs Gale Qobs 
He 0.409 276 0.009639 0.300* 0.298 0.7 | Xe 0.881 206 0.007116 0.469* 0.522 11.3 
0.01011 0.296 0.301 1.7 0.01127 0.736* 0.614 16.6 
0.04097 0.386 0.376 2.6 0.03784 0.794 0.779 1.9 
0.05250 0.376 0.383 1.9 0.06307 0.810 0.816 0.7 
0.06076 0.392 0.386 1.5 0.07540 0.848 0.826 2.6 
0.07714 0.391 0.391 0.0 0.1250 0.865 0.846 2.2 
0.09940 0.385 0.395 2.6 0.1725 0.869 0.854 1.7 
0.1127 0.398 0.397 0.3 0.2111 0.866 0.859 0.8 
0.1450 0.400 0.400 0.0 0.2524 0.849 0.862 1.5 
0.1740 0.398 0.401 0.8 0.2891 0.856 0.864 0.9 
0.1995 0.406 0.402 1.0 0.3352 0.835 0.866 3.7 
0.2210 0.394 0.403 2.3 0.3986 0.864 0.868 0.5 
0.2461 0.411 0.404 1.7 -—~ 
0.2596 0.408 0.404 1.0 Mean 3.7 
eet aie eam “f° H: 0.316 545 0.008717 0.253* 0.261 3.2 
0.3417 0.400 0.405 1.2 0.009725 0.270* 0.266 1.5 
0.3583 0.406 0.405 0.2 0.04916 0.309 0.304 1.6 
; = 0.1013 0.311 0.310 0.3 
, 0.1027 0.316 0.310 1.9 
atone 4.4 0.1536 0.316 0.312 1.3 
0.2039 0.312 0.313 0.3 
Ne 0.718 165 0.01104 0.466 0.462 0.9 0.2479 0.309 0.313 1.3 
0.01374 0.487 0.496 1.8 0.3111 0.309 0.314 1.6 
0.04714 0.642 0.634 1.2 0.3434 0.308 0.314 1.9 
0.06002 0.659 0.650 1.4 ~ 
0.1018 0.685 0.675 1.5 Mean 1.5 
0.1709 = 0.700 0.691 1.31 bh, 0.396 697 0.01162  0.350* 0.353 0.9 
0.2322 0.705 0.697 1.1 0.04888 0388 0.385 08 
0.2900 0.706 0.701 0.7 , “i — . 
0.3669 0:700 0.703 04 0.06222 0.386 0.387 0.3 
03040 0-703 0-704 01 0.1041 0.395 0.391 1.0 
04671 0691 0706 22 0.1056 0.392 0.391 0.3 
; ' : a 0.1439 0.394 0.392 0.5 
Wien de 0.1767 0.398 0.393 1.3 
a 0.2048 0.399 0.393 1.5 
0.2112 0.396 0.393 0.8 
A 0.870 195 0.01095 0.596* 0.596 0.0 0.2903 0.389 0.394 1.3 
0.01312 0.612 0.627 2.5 0.3470 0.389 0.394 1.3 
0.05125 0.807 0.793 1.7 0.3804 0.387 0.395 2.1 
0.1045 0.833 0.831 0.2 0.4148 0.385 0.395 2.6 
0.1589 0.851 0.844 0.8 a 
0.1643 0.847 0.845 0.2 Mean 1.1 
0.2157 0.848 0.852 0.5 | N, 0.783 392 0.009121 0.620* 0.612 1.3 
0.2252 0.856 0.852 0.5 0.01146 0.625 0.641 2.6 
0.3015 0.853 0.857 0.5 0.04950 0.758 0.745 17 
0.3497 0.845 0.859 Lz 0.09147 0.765 0.762 0.4 
0.3883 0.842 0.860 2.1 0.1008 0.761 0.764 0.4 
nae 0.1106 0.761 0.765 0.5 
Mean 1.0 0.1524 0.784 0.770 1.8 
0.2107 0.776 0.774 0.3 
Kr 0.878 126 0.009823 0.445* 0.487 9.4 0.2181 = 0.775 (0.774 ry 
0.2973 0.767 0.777 1. 
0.04863 0.776 0.757 2.7 0.3255 0.763 0.777 18 
0.06158 0.800 0.780 2.5 _ — bes 
0.1040 0.826 0.818 1.0 ve 
0.1077 0.842 0.820 2.6 en 
0.1451 0.850 0.836 1.6 | CO 0.778 702 0.008454 0.665* 0.666 0.2 
5 0.1818 0.857 0.844 1.5 0.05471 0.756 0.758 0.3 
0.2173 0.855 0.850 0.6 0.09192 0.768 0.766 0.3 
0.2451 0.843 0.854 1.7 0:1216 0.778 0.769 1.2 
0.2914 0.841 0.858 2.0 0.1778 0.776 0.772 0.5 
0.3484 0.841 0.861 2.4 0.2064 0.781 0.772 1.2 
0.3764 0.840 0.863 2.7 0.2509 0.765 0.774 1.2 
0.4165 0.840 0.864 2.9 0.3017 0.768 0.774 0.8 
0.4766 0.826 0.867 5.0 0.3311 0.767 0.775 1.0 
Mean 2:8 Mean 0.7 
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TaBLe I1.—Continued. 








8 
& 
Et. 


| 2 


Gas a. b(mm) P (mm) obs 


O, 0.800 242 0.008557 0.534* 
0.05029 0.742 
0.09046 0.787 
0.1226 0.782 
0.1520 0.784 
0.2080 0.792 
0.2306 0.779 
0.2579 0.782 
0.3054 0.770 
0.3259 0.777 





me SSONO- 


| : 


te OOO LO tne 


Mean 1.3 





a, values are taken from measurements using a 
platinum wire,* they are assumed to be essen- 
tially independent of the wire used because of 
the saturation of the metal surface. 

Inspection of Eq. (3) shows that the elimina- 
tion of a produces an extremely small effect on 
values of @ calculated at pressures where an 
appreciable fraction of the surface is covered, 
since the small ao term is, under these conditions, 
multiplied by the term 1/(1+0p) which ap- 
proaches zero as the surface approaches satura- 
tion. The simplified theoretical equation for the 
pressure variation of the accommodation coef- 
ficient will therefore be taken as 


bp 
= ——— Ae 
1-+bp 


1/a=1/a,+1/(bpa.). (5) 


Values of 6 and a, were computed from Eq. 
(5) for each gas listed in Table II from the ob- 
served accommodation coefficients at the tabu- 
lated pressures. The method of least squares was 
used with the standard weighting procedure® 
necessitated by the inversion of Eq. (4) to produce 
Eq. (5) which is linear in the reciprocals of the 
desired constants. The starred a}, values differ 
from those originally published* because of a 
correction required in the original meter readings 
from which these values were computed. These 
readings fell on a different scale than the others 
and were found, upon calibration of the meter, 
to be in error by about 1 percent. This relatively 


* J. B. Scarborough, Numerical Mathematical Analysis 
(The Johns Hopkins Press, Baltimore, 1930), p. 370ff. 


small error has a large effect on the accommoda- 
tion coefficients at the low pressures involved 
because of the small difference between the 
meter reading in vacuum and at the given 
pressure. The entries under aeaie Were computed 
from Eq. (4) using the b and a, values given in 
the table. 

The agreement between the observed and 
calculated accommodation coefficients, as ex- 
pressed by the absolute percentage deviations, 
seems quite satisfactory in all but three cases, the 
first value for krypton and the first two values 
for xenon. This may be caused by the difficulty 
of obtaining precise values of ap; in these cases, 
where the low pressure and the low free molecule 
conductivity combine to give a very small dif- 
ference between the energies required to heat 
the platinum wire to a given temperature in 
vacuum and in the gases at the recorded pres- 
sures. It is also possible, but less likely, in the 
opinion of the author, that these larger devia- 
tions result from the assumption that ap is 
negligible in comparison with a, since, in this 
case, it would be difficult to explain the opposite 
sign of the deviations in the first two entries for 
xenon. 

It is not possible to make a quantitative com- 
parison between the results of Table II and those 
of other investigators who have determined 
accommodation coefficients on platinum, since 


TABLE III. 











Investigator yas @obs Pobs (mm) Peaic (mm) 
Rolf : 0.170 
Thomas and Olmer 2 0.220 0 
Thomas and Olmer 2 0.295 0 
Thomas and Olmer 0.238 0 
Thomas and Olmer N 0.57 0 
Thomas and Olmer , 0.89 0 
Thomas and Olmer 2 0.74 0 





1 X10-% to 1 X10? 2.6 X10~5 
to 5 X10 5.0 X10-8 
to 5X10 4.2 X10-% 
to 5X10 4.2 X10- 
to 5 X1072 2.4 X1072 
to 5 X10 

to 5 X10 








only Rolf* and Thomas and Olmer’ have pub- 
lished details which show that the heat flow 
problem has been properly analyzed, and in both 
cases only single values of the accommodation 
coefficient are reported for a wide range of 
pressure. Furthermore, the results of Thomas 
and Olmer, particularly those near unity, are 
uncertain because of an error in one of their 
formulas.* In Table III the accommodation 
coefficients reported by Rolf and by Thomas and 
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Olmer have been combined with the constants of 
Table II to ascertain if the pressures calculated 
from Eq. (4) fell within the ranges cited by these 
authors. 

The author feels that the qualitative agree- 


DUFFEY 


ment in six of the seven cases in Table III tends 
to support the present simple adsorption theory 
but that more reliable confirmation must await 
additional experimental results of the type given 


in Table II. 
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The hybridization of s and p orbitals to give non-tetrahedral orbitals is discussed. The 


. results are applied to cyclic hydrocarbons. For cyclopentane the theory yields an H-C-H 


HE method of directed valence bonds as 

developed by Pauling! is based on the 
following postulates. A single bond arises from 
the interaction of two unpaired electrons of 
- Opposite spin derived from different atoms. Of 
two orbitals in an atom the one which can over- 
lap more with the orbital of another atom will 
form: the stronger bond with that atom. The 
direction of this bond will be that direction in 
which the orbital wave functions of the two 
atoms overlap as much as possible. In applying 
these postulates to the carbon atom Pauling 
assumed that orbitals were formed from a com- 
bination or hybridization of s and # orbitals. 
The hybridization was made in such a way as 
to yield the maximum bond Strengths. 


HYBRIDIZATION IN THE NON-TETRAHEDRAL 
CARBON ATOM 


Let us now consider how this theory is to be 
altered if we specify one of the valence angles of 
the carbon atom. Such a restriction is present in 
compounds such as cyclopropane, cyclobutane, 
and cyclopentane. It will be assumed that a 
certain amount of symmetry exists. Thus num- 
bering the bonds of carbon as 1, 2, 3, and 4, it 
will be assumed that bonds 1 and 2 are equivalent 
and that bonds 3 and 4 are equivalent. 

1L. Pauling, J. Am. Chem. Soc. 53, 1367 (1931). See 


also W. Heitler and F. London, Zeits. f. Physik 44, 455 
(1927); J. C. Slater, Phys. Rev. 37, 481 (1931). 


angle of 111° 1’, and for cyclobutane an H-C-H angle of 180°. 





It is assumed that the radial parts of the wave 
functions y,, pz, Wey, Wz are so similar that their 
differences are negligible. The angular parts, 
normalized to 47, are given by the equations 


s=1, (1) 
p.=V3 sin 3 cos ¢, (2) 
py =V3 sin # sin ¢, (3) 
p.= V3 cos #, (4) 


where # and ¢ are the angles used in polar 
coordinates. 

The angular part of each of the bond-forming 
eigenfunctions is approximated by functions of 
the form 

fi=ast+bip.tcip,+dipz, (5) 
where 7 is 1, 2, 3, or 4. The coefficients a;, 0;, ¢;, 
and d; are determined by the normalization 
condition 


a?+6?+c?+d?7=1, (6) 


and by the orthogonality condition 


a,a;+6,b;+cc;+dd;=0, tj. (7) 


The direction of the maximum of the first 
bond eigenfunction is arbitrary. It is chosen to 
lie along the x axis. Hence 


fi=ais+(1—a,’)'p,, (8) 


where }; has been replaced by (1—a,’)! so that 
Eq. (6) might be satisfied. 













ar 





(1) 
(2) 
(3) 
(4) 
olar 


1ing 
s of 


(8) 


that 
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The direction of the maximum of the second 
bond eigenfunction can without loss of generality 
be chosen to lie in the xz plane. Thus 


fe=a25 + bop, +dop.. (9) 
From Eq. (6) 
a2? +b.?+d.? = 1. (10) 
From Eq. (7) 
2;d2+0);b2.=0. (11) 
Since bonds 1 and 2 are assumed to be equivalent 
bonds 
d2=d). (12) 
Hence 


be = —[ay?/(1—a,’)?], (13) 


and 


d,=[(1—2a,”)/(1—a,’) }}. (14) 
In a similar way it is found that 
d3=a4=(4—a,’)?, (15) 
bs=bs= —ai{[(1—2a;’)/(2(1—a,’)) ]'}, (16) 
C3= —Cy=1/Vv2, (17) 
d3=d4= —a,/[2(1—a,’) }}. (18) 


Now let us relate the direction of the maximum 
of the second bond eigenfunction to the value 
of a;. First substitute Eqs. (1), (2), (4), (12), 
(13), and (14) in Eq. (9). Then let ¢=180° and 
differentiate fz with respect to &%. Setting this 
derivative equal to zero then gives the result 


tan d=a,"/(1—2a,’)}. (19) 


Since the direction of the maximum of the first 
bond eigenfunction lies along the x axis, the 
angle between the two maxima equals 90°+#. 


CYCLOPENTANE 


It is assumed that all the carbon atoms in 
cyclopentane lie in a plane and form a sym- 
metrical ring. Also it is assumed that the direc- 
tions of the bonds coincide with the directions 
of the maxima of the bond eigenfunctions. 
Consider one carbon atom in the ring and choose 
the coordinate axes so that f; and fe represent 
the angular part of the eigenfunctions which 


form the bonds with the neighboring carbon 
atoms. The value of a; can be obtained using 
Eq. (19) and taking #=18°. The value of a; is 
then calculated using Eq. (15). Now consider 
the axes rotated so that the direction of the x 
axis coincides with the direction of a C-H bond. 
The new value of a; is equal to the value of a; 
in the old coordinate system. The value of the 
H-C-H angle can now be calculated using Eq. 
(19). The final result is 


HCH =111° 1’. (20) 


It is evident that this value differs appreciably 
from the tetrahedral angle 109° 28’. 


CYCLOBUTANE 


If a, is taken equal to zero, it is found by 
Eq. (19) that 3 equals zero and hence the angle 
between the directions of the maxima of f; and fe 
equals 90°. Proceeding as before one obtains 


<~ HCH = 180°. (21) 


CYCLOPROPANE 


The calculations of Herzberg? seem to support 
the view that the HCH angle in cyclopropane 
is a straight angle. This indicates an electronic 
structure similar to that described for cyclo- 
butane. However now the directions of the 
maxima of f; and fz would not coincide with the 
C-C bond directions. Instead there would be a 
difference of 15° between each C-C bond. and 
the maximum of f; and fe, respectively. 


DISCUSSION 


The results obtained are to be regarded as 
rough approximations. The principal effect neg- 
lected is probably the repulsion between different 
groups in the molecule. In the case of cyclobutane 
the HCH angle is probably altered considerably 
by the repulsion between the hydrogen atoms 
and the neighboring carbon atoms. 

2 G. Herzberg, Infrared and Raman Spectra of Polyatomic 


Molecules (D. Van Nostrand Company, Inc., New York, 
1945). 
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The films of n-heptane formed by adsorption on the surface of graphite between 25° and 40°C 
exhibit a first- and a second-order phase transition. The first-order transition occurs between 
the gaseous and liquid expanded phases. The second-order transition occurs between the 
liquid expanded and liquid intermediate phases, and between the gaseous and liquid inter- 
mediate phases above the critical temperature. Since the behavior of the film is the same as ' 
that observed with insoluble films on aqueous subphases, the nomenclature for the latter is 























used. The critical constants for the gaseous film were determined as: temperature, 31°C; 
area, 400A? per molecule; and film pressure, 1.05 dyne cm™. These values are different from t 
those found for x-heptane on ferric oxide, which shows that the solid as well as the gas plays an [ 
important role in the determination of the adsorption isotherm. The heat of transition for the ] 
first-order change is estimated as 13,000+5000 cal. mole“ at 25°C. It is found at corresponding | 
temperatures, within this large experimental error, that the heat evolved in the first-order ’ 
transition is the same on graphite and ferric oxide. 
C 
, 
I. INTRODUCTION solids may undergo first-order phase transitions 
ECENTLY, the writers have presented ©” increase of pressure from the gas phase. These 
experimental evidence which shows that are analogous to those which occur in insoluble 
films formed by adsorption on the surfaces of films on an aqueous subphase. Data were pre- 8 
09 sented for two systems: n-heptane on silver! and t 
J e e 
| | n-heptane on ferric oxide.” Z 
ost Of the two systems, that which involves ferric 3 
oxide was the more extensively investigated. The 2 
-— ™ system exhibited all the critical temperature 0 
rc phenomena in the film pressure (m)-area (¢) E 
a relationships associated in general with a first- Is 
+ OO Fic. 1. The isotherms bet ws 
“ 25° and 40°C of aia a order transition. It was also shown that the \ 
2 05-— oxygen-free graphite with less compressibility x of the film formed by adsorp- V 
o than 0.004 percent ash to pres- tion on the surface of the solid is 
= sures of 0.025 mm. At 25°, 7. 
= 94 28.5°, and 30°C the isotherms Vaud dv . 
3 40%--a-- exhibit a discontinuity in the c= 0: 
| volume adsorbed as a function of RTv?\a In p/ 7, . 9 
a —| the pressure. This indicates a 
= —— phase transition of the where Vy is molar gas volume at standard tem- st 
8 o2t— is perature and pressure; 2, the specific area of the th 
solid; R, the gas constant; 7, the absolute tem- cl 
01 aa perature; v, the volume of gas adsorbed per of 
gram of solid; and p, the equilibrium pressure. 
| | From the criteria for phase changes in films fo 
0 oF 02 presented by Dervichian® and by Harkins and rs 
PRESSURE .MM HG —_—— ri 
naan 1G. Jura, E. H. Loeser, P. R. Basford, and W. D. 

* The writers wish to acknowledge the helpful coopera- Harkins, J. Chem. Phys. 13, 535 (1945). 13 
tion in this work of Dr. Orlan M. Arnold of the Chrysler 2G. Jura, E. H. Loeser, P. R. Basford, and W. D. 
Corporation. This is a part of a larger research project Harkins, J. Chem. Phys. 14, 117 (1946). CI 
under the general direction of William D. Harkins. § D. G. Dervichian, J. Chem. Phys. 7, 931 (1939). 
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Boyd‘ it also was shown that for a first-order 
transition (dv/d In p)7,» becomes infinite, while 
for a second-order transition there is a finite 
discontinuity in the derivative. 

In the present paper, which considers the ad- 
sorption of n-heptane on graphite at low relative 
pressures, the experimental data indicate that 
both a first- and a second-order transition occur 
at pressures below 0.040 mm at temperatures 
between 25° and 40°C. 


Il. EXPERIMENTAL METHODS AND MATERIALS 


The experimental procedures and materials 
used have been fully described in previous 
papers.*® The graphite used is that which con- 
tains less than 0.004 percent ash. It was furnished 
by the courtesy of Dr. Lester L. Winter of the 
National Carbon Company. As in the previous 
paper, 4.22 m*g! was used as the area of the 
sample, as obtained from an absorption isotherm 
of nitrogen at —195.6°C, and calculated by the 
relative method of Harkins and Jura. 


III. EXPERIMENTAL RESULTS 


The adsorption isotherms of n-heptane on 
graphite are exhibited in Fig. 1. These were ob- 
tained at pressures up to 0.025 mm at 25.00° 
+0.02, 27.00°+0.02, 28.50°+0.02, 30.00°+0.02, 
35.00°+0.05, and 40.00°+0.05C. At 25°, 27°, 
28.5°, and 30°C, the isotherms exhibit a first- 
order phase change at a pressure of 0.014 mm 
Hg. The behavior of the isotherms in this region 
is exactly analogous to that of those observed 
when ferric oxide is the subphase. The best 
values for the critical constants on graphite are: 
T=31°C, o=400A* per molecule, and r=1.05 
dyne cm~!. The corresponding values for ferric 
oxide are: T=29°C, «=900A? per molecule, and 
™=(0.45 dyne cm~!. These two sets of values 
show the dependence of the critical constants on 
the nature of the subphase. Obviously, the 
critical constants must depend also on the nature 
of the gas. 

The heat evolved in the first-order trans- 
formation of the film of u-heptane on graphite 
neat D. Harkins and G. E. Boyd, J. Phys. Chem. 45, 20 

°C Jura and W. D. Harkins, J. Am. Chem. Soc. 66, 
1356 (1944). 


*W. D. Harkins, G. Jura, and E. H. Loeser, J. Am. 
Chem. Soc. 68, 554 (1946). 
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Fic. 2. The same as Fig. 1, but to pressures of 0.100 mm 
Hg. These curves exhibit at pressures above 0.025 mm a 
second-order phase transition of the film. 


at 25°C is estimated at 13,000+5000 cal. mole-, 
which is practically the same as that on ferric 
oxide at the same temperature: 12,000+5000 
cal. mole. If the comparison were made at 
corresponding temperatures for the film on 
graphite and ferric oxide, the heats evolved in 
the phase transformation would be the same. The 
experimental facts seem to indicate that the gas- 
gas interaction cannot be neglected; even when 
the adsorbed molecules are separated by com- 
paratively large distances. 

Figure 2 exhibits the isotherms to a pressure 
of 0.100 mm. The experimental points for the 
discontinuity and lower pressures have been 
omitted, since they are shown in Fig. 1. Each of 
these curves (except that at 27°C, where the data 
are insufficient) exhibits a sudden increase in 
(dv/0p)7,z at pressures decreasing from 0.038 mm 
at 25°C to 0.025 mm at 30°C and then again 
increasing to 0.032 mm at 40°C. The data in 
this region are completely reproducible, and the 
values do not change with time. In several of 
the determinations in this region two consecutive 
points are found at the same pressure. It is 
apparent that this rapid change in (dv/dp)7,z 
represents a second-phase change in the film. 
The isotherms below the critical temperature 
(T=30°C) cannot be used to decide this point, 
since too few experimental points can be ob- 
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Fic. 3. A cabinet projection of the r—o—T surface of 
the film of n-heptane on graphite. In this projection the 
correct —o relations are maintained. The temperature 
axis is drawn at 45°. Thus, an increase of 5°C causes an 
apparent increase of molecular area of 100A? per molecule 
and an apparent increase of 1 dyne cm™ in the film 
pressure. The surface shows four regions: (1) the gas, G; 
(2) the “liquid expanded,” Lz; (3) the “liquid inter- 
mediate,’ L;; (4) the two-phase region, gas+“‘liquid ex- 
panded.” 


tained between the identified first-order transi- 
tion at 0.014 mm and the possible second 
transition. At 35° and 40°C, however, a definite 
decision can be made. At these temperatures a 
study of the data shows that a transition does 
occur, the question being solely one of order, i.e., 
first or second. The fact that in some cases two 
points are observed at the same pressure might 
indicate a first-order change with an exceedingly 
small change of area. A study of the pressure (7)- 
area (a) curves indicates, however, that the only 
reasonable assumption is that a second-order 
change occurs. If the transition were first order, 
the compressibility of the phase formed by the 
transition would at first decrease and then 
increase as the pressure is increased. This be- 
havior is unreasonable from both the two- and 
three-dimensional behavior of matter under 
pressure. This difficulty is obviated if it is 
assumed that order of the change is second rather 
than first. The —o curves are exhibited in Fig. 3. 

If the second transition is considered as first 
order, the heat relations on this basis are peculiar. 
The AH for the second transition on the assump- 
tion that the change is first order, becomes posi- 
tive below the critical temperature and negative 
above this temperature. Also, at 25° AH=300 
cal. mole, which increases to 27,000 cal. mole! 
at 30°C and decreases to —700 cal. mole at 35° 
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Fic. 4. The —T phase diagram of the film of n-heptane 
on graphite. The triple point is placed at the critical tem- 
perature and critical film pressure. 


and 40°. The AH for the process of transforming 
the vapor to the third phase that appears would 
be —12,000 cal. mole! at 25°C, which would 
increase to 26,000 cal. mole! at 30°C and then at 
above the critical temperature change to — 700 
cal. mole. This behavior is also unlikely. Thus, 
it is almost certain that the phase change under 
consideration is of the second order. The fact 
that this change exhibits an appearance some- 
what analogous to a change of the first order is 
presumably due to an insufficient sensitivity of 
the pressure measurements. 

By analogy with the behavior of insoluble 
films on water, the phase formed from the gas 
by the first-order change would be the liquid 
expanded (L,.), while that formed by the second- 
order change from the gas and the analog of the 
liquid expanded would be considered the liquid 
intermediate (L;) phase. The phase diagram of 
this system is given in Fig. 4. The diagram is in 
complete accord with our present knowledge of 
the behavior of films and also of matter in three 
dimensions, when the pressure, temperature, area 
(or volume) relations are considered. 

Figure 4 does not represent all of the phases 
which appear in this system. At 25°C, where the 
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entire isotherm has been determined, the “‘liquid 
intermediate’’ phase transforms into a condensed 
phase at a film pressure of about 20 dyne cm™. 
The nature of this transition has not as yet been 
determined. 

The second-order change reported here has 
also been found with other solids when n-heptane 
is the adsorbate, and also with n-butane and 
triptane at 0°C on several solids. The work done 
to date with n-butane indicates that the ob- 
served phenomena are the same as those found 
for graphite above the critical temperature. The 
work of Orr’ strongly suggests that similar 
transformations take place when argon, nitrogen, 
and oxygen are adsorbed on ionic solids at tem- 
peratures in the neighborhood of 77°K. The 
writers are planning to repeat his work and, by 


TW, J.C. Orr, Proc. Roy. Soc. London 173A, 349 (1939). 


taking more points, to determine whether or not 
these are phase transitions of the second order 
in the low pressure region. First-order transitions 
with nitrogen at these temperatures have been 
observed. 

The results with graphite show, even when the 
adsorbed film is above the critical temperature 
of the liquid expanded film, that the increase 
of the pressure of the film transforms the gas 
into a phase which, below the critical tempera- 
ture, does not appear until much higher pressures 
are attained. This behavior has its counterpart 
in the solidification of helium and other gases 
above their critical temperatures. Thus, in films 
formed on the surfaces of solids by adsorption, 
phase transformations from the gas are possible, 
even above the critical temperature of the liquid 
expanded. 





Errata: Statistical Mechanics of Transport Processes 


I. General Theory 


[J. Chem. Phys. 14, 180 (1946)] 


Joun G. KirkKwoop 
Baker Laboratory, Cornell University 


ip this article an important distinction was 
made between the distribution function f” 
(p, q; ¢) at time ¢ and its time average f™ 
(p, q; ¢) over an interval 7 subsequent to the 
instant t. The two functions are related in the 
following manner, 


i 17 
J (p, a; 4) ~~ free. q, t+s)ds. 
0 


Unfortunately, serious confusion was caused by 
a series of typographical errors introduced be- 
tween galley and page proof, in which the bar 
Was omitted at random from the functions 
J™(p, q;¢). These errors are to be corrected as 
follows: 


1. f™ to replace f® on the left side of the second 
of Eqs. (14). 

2. f™ to replace f™: 

(a) On left side of second of Eqs. (20). 

(b) In integrand of Eq. (21). 

(c) On left side of the third of Eqs. (46). 


(d) On left side of the differential equation, 
Eq. (47). 

(e) On right side of the first of Eqs. (48). 

3. f® to replace f™: 

(a) In the integrands on the right sides of 
Eqs. (23), (24), (25), (26), and the second 
of Eqs. (60). 

(b) On the left side of the differential equa- 
tion, Eq. (49). 

(c) On both sides of the differential equation, 
Eq. (59). 

(d) First sentence, last paragraph, page 194. 


Two additional errors of a minor character 
are also to be noted: 


4. In the first of Eqs. (70), Ci, the concentrations 
in moles/cm’ should replace C;, the concen- 
trations in molecules/cm‘*. 

5. The expression 62np;/m; should replace 
67nmp; at the end of the eighteenth line from 
the bottom of the second column, page 198. 
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On the Weak Associations in Water 
and in Aniline 
R: PARSHAD 


Physical Laboratory, Council of Scientific and Industrial Research, 
University Buildings, Delhi, India 
January 22, 1946 


| peng concept of the hydrogen bond has ex- 
plained a vast field of experimental data.' In this 
note we wish to show that his concept can be further used 
to explain the very interesting property of the comparative 
weakness of associations in water. For example, taking 
the energy of the hydrogen bond as the thermal energy 
required to rupture it, the energy of the hydrogen bond in 
water is 4.5-5 kcal./mole, compared to that of 6 kcal./mole 
for alcohols and 7-8 kcal./mole for carboxylic acids. It is 
probably not well realized that the “abnormal” properties 
of water are not only caused by its vast network of asso- 
ciations, but also because these associations are weak. It 
is not the mere presence of hydrogen bonds but their pro- 
gressive breaking with increase of temperature, or pressure 
etc., which causes the open structure of water caused by 
extensive hydrogen bonding to crumble and bring in the 
abnormal properties. For that matter, carboxylic acids 
are hardly abnormal in the liquid states. Because of the 
nearly non-breaking of the strong dimers, the latter behave 
as single molecules until even vaporization. Indeed the 
vapor state rather than the liquid state of these acids may 
be called associated in the sense that water is in the liquid 
state. 

While the extensive nature of the association in water 
may be explained as caused by favorable steric factors, 
(viz. all atoms in a state to be engaged in hydrogen bonds, 
the donating power of O in terms of H’s to be equal to the 
accepting power, and the ability to be able to form hy- 
drogen bonds on all sides (tetrahedral bonding), etc., the 
weak nature of these associations may be explained as 
follows. 

According to Pauling, in the O—H valency link, there is 
resonance between the covalent state and the ionic state. 
A pair of electrons is at one time shared between O and H 
(covalent state) and in the next instant is all attracted on 
to O by its electronegativity, leaving the H a bare proton 
(ionic state). In the state of the ionic character of O—H, 
the proton H attracts and is attracted by an O of another 
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O—H in virtue of the latter’s total negative charge (when 
in the ionic state) and any unpaired electron pairs on it. 
Here we consider two factors: 

(1) As one O—H link of a water molecule is ionic, the O 
is negative, because of the pair of electrons that it has at- 
tracted on to itself. This negative charge on O will dis- 
courage the attraction of another pair of electrons from 
the other O—H link. In other words, the negative charge 
on O, formed because of the ionic character of O—H link, 
decreases the electronegativity of O and so the ionic 
character of the second O—H link. As both the O—H 


. links are of course indistinguishable, this means that the 


ionic character of both the O—H’s are reduced, and so 
the hydrogen bonds formed will be weak. The liquids like 
alcohols have one H to one O and so do not suffer this 
weakening of the ionic character. 

(2) Consider the effect of an H, on O, not of its own valency 
link but that with which it is involved in a hydrogen bond. 


In the above figure, the atoms concerned are in bold face. 
The H proton will attract electrons on O, and so the latter 
will increase its effective nuclear charge for its own H (on 
the right of the figure). The electronegativity of O is in- 
creased and so the ionic character of O—H, leading to in- 
crease of hydrogen bonding power. 

The effects of (1) and (2) run in counter directions but 
evidently the effect of factor (1) is greater than that of 
(2), as in the former the mutual effect of O and H in the 
same valency link is considered. We may note here that 
derivation of factor (2) is contrary to that of Pauling. 
According to him the O—H link concerned is strengthened 
(or its ionic character decreased). However, Pauling’s 
derivation all the more supports the inference of weak 
associations in water. 

The above considerations for water can also be used to 
explain why aniline containing an amino group forms no 
or very little hydrogen bonds, when it is also considered 
that N has lower electronegativity than O, and because 
of the extra phenyl group on N, lower number of unshared 
pair of electrons are left on N for attraction by a distant 
proton for forming a hydrogen bond. But the electron at- 
tracting property of the phenyl group itself will contribute 
some additional electronegative power to N, and because 
of this aniline may be more associated, if at all, than 
ammonia. 

1L. C. Pauling, Nature of the Chemical Bond (Cornell University 


Press, New York, 1940). 
2 Reference 1, p. 326. 





A Note on Communal Entropy. Remarks on a 
Paper by Henry S. Frank! 


O. K. RICE 
University of North Carolina, Chapel Hill, North Carolina 
April 8, 1946 


N a recent series of interesting papers, Frank has dis- 
cussed what he aptly calls the pictorial theory of the 
liquid state. These papers included, also, some considera- 
tion of the solid state, and the gas of hard elastic spheres. 
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In the latter discussion, Frank! has come into conflict with 
some statements previously made by me.? On page 487 of 
reference 1, Frank discusses the communal entropy of a 
closely packed three-dimensional gas of hard elastic spheres. 
This I had estimated as 3R per mole on the basis of Eq. 
(2) of reference 2. Frank rejects my Eq. (2) and sub- 
stitutes his Eq. (36), vut grants that the latter is not itself 
satisfactory. My equation, applied first to a simple cubic 
array, rests on the assumption that motions in the three 
directions of space are independent of each other, pro- 
vided the system is sufficiently closely packed. Thus if an 
atom A (hard elastic sphere) is displaced out of the line 
joining its two neighbors, B and C (Fig. 1), it is assumed 
that its motion along this line is unaffected. This is true if 
the displacement is small enough. But even if the atoms 
are densely packed, large fluctuations become appreciably 
probable if the number of atoms in the line is large. It is, 
therefore, necessary to place some arbitrary restriction on 
the total possible displacement of any atom, allowing it 
to be not greater than some finite multiple (which may, of 
course, be as large as we please, if the spheres are sufficiently 
densely packed) of the average distance between atoms. 
Some such condition is necessary in order to keep the sys- 
tem in the assumed simple cubic lattice; if left to itself, 
the system would presumably tend to go over to a lattice 
of the face-centered cubic, or some related type, which 
would have a greater entropy; but this is not the problem 
we are attempting to solve. Even though the simple cubic 
lattice is not a stable lattice unless some restrictions are 
placed on the motions of the atoms, it is still useful to 
consider it, because it may be used as a guide to the cases 
which have a greater resemblance to physical reality. 

Let us consider a sufficiently closely packed simple 
cubic array of N hard elastic spherical atoms packed in a 
cube of length J. Letting K=(2rmkT/h?)*%"?, the total 
partition function of this system will be given, assuming 
classical statistics, by* 


o=Kf{ {- ° + f ded: ° -dxydyidyo- e 

dyndz,dze-++dzy. (1) 
In evaluating the integral we allow the x’s, y’s, and 2’s to 
take all possible values, keeping, however, each numbered 
atom in its particular place in the array. This latter con- 


dition makes it unnecessary to divide by N!, since we do 
not allow exchange of atoms. Since motions in the three 
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directions of space are independent, we may separate the 
x, y, and z parts of the integration, and since these parts 
are exactly alike, we may write 


oren(ff- aes: dey) (2) 


The integral consists of integrations along N! exactly 


similar rows of N!} atoms each. Hence 
3 


ovan(SS- . » f dxdxy: . dy!) (3) 


Except for the restrictions on large displacements this 
integral is just that evaluated by Tonks.‘ It is clear from 
pages 4 and 5 of reference 2 that reasonable restrictions 
can be placed on the fluctuations without appreciably al- 
tering the value of this integral. We, therefore, insert 
Tonks’ value with the appropriate changes (N! for N), 
obtaining 


i 
QO’ = KE (I— Nia)" /NiNBN (4) 
where a is the diameter of the spheres. If N! is large enough 


(ee A » ° ‘ 
1/Ni!=e% (Nt), excluding a factor whose logarithm will 
be negligible. Letting v=/*/N be the volume per atom, 
and setting @=o*/v, the expression becomes 


QY = Ke®% yp ( 1—@! ) aN (5 ) 


entirely equivalent to Eq. (2) of reference 2. As noted in 
reference 2, this is greater by the factor e*” than the par- 
tition function which would be obtained if each sphere 
were confined to its own cell, the cells being all equal in 
size and completely filling the space. The molal communal 
entropy is, therefore, 3R. How to take care of the case 
where the spheres are packed in a cubic close-packed 
(cubic face-centered) type of arrangement, rather than 
simple cubic, is probably sufficiently indicated in reference 
2. This involves a certain estimate, rather than a rigorous 
straight-forward deduction, but it is undoubtedly a good 
estimate and, in any event, in itself, it in no way affects 
the communal entropy. The change to the cubic face- 
centered type may, on account of the blocking of fluctua- 
tions, have a slight effect on the communal entropy, some- 
what similar to the effect due to restricting the size of the 
fluctuations in the simple cubic array. 

While this seems to settle the problem from a formal 
mathematical point of view, it may be possible to bring 
out the difference between a closely packed aggregate of 
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spheres and a dilute gas, for which the communal entropy 
per mole is only R, more graphically. If the spheres are 
very densely packed, under the restrictions noted above a 
sphere cannot get out of its row and column. (For con- 
venience we restrict ourselves to two dimensions.) On the 
other hand, consider the situation which may arise in the 
case of point atoms, as illustrated by Fig. 2. This figure 
shows a square divided into cells by solid lines, one cell 
for each atom. If we assume that the atoms may pass the 
cell boundaries, we can get a situation such as shown in 
Fig. 2. This may be considered as produced in the main by 
fluctuations in the top row only. It can also be produced 
by another set of fluctuations as is indicated by the broken 
lines. Since the two sets produce indistinguishable con- 
figurations, they give rise to only one state of the system. 
In the case of closely packed spheres, however, the two 
sets of fluctuations, if they were of the same order of mag- 
nitude relative to the free space (therefore taking the mole- 
cules out of their individual cells), would, nevertheless, 
give only small displacements relative to the diameter of 
the spheres. The two sets of small displacements would 
give rise in this case to distinguishable configurations. It 
is, therefore, seen that the closely packed spheres can have 
a greater communal entropy than the dilute gas, because 
in the latter case row and column fluctuations interfere 
and overlap, giving rise to indistinguishable configurations. 

Frank has quite properly called attention to the arbi- 
trary nature of communal entropy. It depends upon the 
size of the “natural” cell which one chooses in the first 
place, when one assumes that each atom is assigned to its 
own cell. However, it is agreed that, insofar as possible, 
the cells should be assigned so that they are of equal size, 
and completely fill the space available. Then the assign- 
ment is easy in the case of a densely packed gas of hard 
elastic spheres. When the atoms exert forces on each other, 
it becomes a little more difficult. I made an estimate of 
communal entropy in this case because I wished later to 
combine the partition function for a gas of hard elastic 
spheres with that of a Debye solid in setting up the parti- 
tion function for a liquid, assuming each part carried its 
own communal entropy factor. It seemed reasonable to do 
this if the excitation of communal entropy was comparable 
in the two cases. Otherwise some special means of averaging 
might have been necessary. Though the Debye solid was 
found to have a communal entropy only 77 percent as 
excited as the dense gas of hard spheres, this did not seem 
sufficiently different to render any special computation 
necessary. 

In order to make a fair comparison of the communal 
entropy in the two cases, it seemed necessary to define 
the ‘‘natural’’ cell for a Debye solid in a manner as nearly 
analogous as possible with the definition for the hard elastic 
spheres. If, in the hard-sphere case, all atoms are fixed in 
position except one, this one has a free volume 2° times as 
large as if it were confined to its own cell. It was, therefore, 
assumed that if all the atoms of a Debye solid were held 
fixed except one, this one also would move in an effective 
volume equal to 2* times its natural cell. It is known that 
this atom would vibrate with the root-mean-square value 
of the frequencies of the normal modes. Therefore, it seemed 
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reasonable to say that a solid in which all atoms moved 
with this same mean frequency (Einstein model) was one 
in which the molal entropy was 3R In 2 times as great as 
it would be if each atom were confined to its “natural” 
cell, and the communal entropy of the Debye solid was 
calculated by comparison with the Einstein model, on 
this basis. In any event, it seems more ‘‘natural’’ to me to 
define the ‘‘natural’’ cell by means of an average frequency, 
Frank prefers the Debye maximum frequency, which seems 
“artificial” to me. However, nothing can be proved about 
this. I merely wish to emphasize that I had a real reason 
for choosing the natural cell as I did in view of the purpose 
I had in mind. 

Errata to reference 2. On page 3 of reference 2 (last sen- 
tence beginning in the first column) where it reads, “If 
this result is multiplied by the momentum integral and 
divided by h,’’ replace h by h*. This gives what we have 
called Q¥ (though in this case it is for a one-dimensional 
array) and Q (correctly given by Eq. (1) of reference 2) is 
the partition function referred to a single atom. 

At the end of page 6 and beginning of page 7, it was in- 
advertently stated “that the partition function per mole- 
cule of a Debye solid is e!kT/hv».”’ This is the partition 
function per degree of freedom and was used as such. 

1H.S. Frank, J. Chem. Phys. 13, 478 (1945). 

20. K. Rice, J. Chem. Phys. 12, 1 (1944). For some errata see ibid. 
12, 521 (1944) and the end of this article. 

3Q% is called Q by Frank; to be consistent with my own Eq. (2) 


I retain O for the partition function referred to a single atom. 
4L. Tonks, Phys. Rev. 50, 955 (1936). 





Postscript to the Foregoing Note by O. K. Rice 


HENRY S. FRANK 
Lingnan University, Canton, China 
April 8, 1946 


GREATLY appreciate the opportunity to see Professor 

Rice’s ‘‘Note on Communal Entropy” in advance of 
publication, and to place on record as promptly as possible 
my agreement with the derivation he has given of the 
limiting partition function for a simple cubic array of hard 
attractionless spheres compressed toward the limit of com- 
plete mutual contact, and with the conclusion which it 
establishes, that 3R is the correct theoretical fluctuation 
entropy for this case, appropriate to the choice of /*(1 
—6)3/N as the elementary cell. Rice’s derivation also 
makes it possible to see where my “superficially straight- 
forward” generalization of the one-dimensional case to 
three dimensions fails in this application. The latter corre- 
sponds, for the present model and for distinguishable par- 
ticles, to letting each atom permute only with other atoms 
in the same linear arrays with it in the x, y, and = direc- 
tions. A correct method of counting, permitting all 1! 
permutations, reproduces both Rice’s mathematics and 
his result. 

When the system becomes more expanded, the motions 
become too complicated for simple analytical discussion. 
There are doubtless several ways of “breaking down” the 
process by which the fluctuation entropy decreases as @ 
decreases, and the picture given here by Rice is very in- 
structive in throwing light on one of them. 
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Rice and I have never disagreed 5n the physical occur- 
rence of a considerable amount of fluctuation entropy in a 
real solid for which any intuitively “natural” choice of 
cell size is made. I agree with him completely as to the 
“naturalness” of choosing a definition of cell size which 
will be most appropriate to the purpose one has in view. 
I now agree with him also in accepting as significant the 
closeness of the agreement between the fluctuation entropy 
found, for his cell size, in solid argon and the theoretical 
value 3R which he has established. 





The Tetrahedral P,; Molecule 


JAMES R. ARNOLD 
Institute for Nuclear Studies, University of Chicago, Chicago, Illinois 
April 15, 1946 


HE methods developed by Pauling! in his treatment 

of the tetrahedral bond orbitals of carbon (and ap- 
plied by him and others to a number of individual cases) 
have been extended and systematized by Kimball? into 
the theory of directed valence. This scheme serves to pre- 
dict the bond configurations for covalent molecules with 
remarkable success. 

One of the few remaining anomalies is the tetrahedral 
P, molecule (with its analog Ass). The bond angles for a 
trivalent P atom might be expected to be somewhat over 
90°, as typical p bonds, and, in fact, all other compounds 
which have been investigated show angles in the neigh- 
borhood of 100°%, Pauling* has suggested that the P, mole- 
cule is merely a strained structure, but the great amount 
of distortion involved seems to render this explanation 
improbable. The purpose of this note is to point out a set 
of orbitals with more suitable bond angles. 

The symmetry group of the three bonds extending from 
one P atom is C3,, and from group theoretical considera- 
tions we can reduce the possible orbitals to p*, sp’, pd?, 
sd’, The bond angles for p* are obviously 90°, while for 
sp? and sd* they are 120°. The combination p,4,4,,, as for 
example in the set: 


¥1=1/V3yp.+ V2/V3Ydz2 
y2= 1/vV3yp.— 1/(6)bvd,2+ 1 /V2pdyz 
¥3= 1/vV3yp.— 1/(6)4var2— 1/v2pdyz 


yields bond angles of 66°26’. This is the smallest possible 
angle for orbitals involving s, p, and d states only. The 
bond-forming strength (in Pauling’s terminology) is 2.29, 
compared to 1.73 for pure orbitals. 

Some distortion is still required to compress these angles 
to 60°, but the strain is not now excessive. The advantage 
gained in replacing a triple bond by three single bonds 
may be used to account for the promotional energy re- 
quired to produce a pd? structure from the p* ground state, 
and the stability of P, relative to Ps. 

The above explanation sheds no light on the 60° bond 
angles found in carbon compounds, but it does account 
for the absence of a tetrahedral nitrogen molecule. 

'L. Pauling, Proc. Nat. Acad. Sci. 14, 359 (1928). 

3G. Kimball, J. Chem. Phys. 8, 188 (1940). 


- Pauling, The Nature of the Chemical Bond (Cornell University 
Press, New York, 1944), second edition, p. 80. 
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Ultra-High Frequency Denaturation of Proteins 


GEorGE A. BoyD 
University of Rochester, Rochester, New York 
April 18, 1946 


LECTROMAGNETIC waves in the ultra-high fre- 
quency range offers a new method for the study of 
protein-type molecules in the following suggested manner. 
From dielectic studies it is known that protein molecules 
are not rotated by frequencies above 10’ cycles per second 
(below 30 meters). However, atomic and electronic polar- 
ization of the proteins do occur at the higher frequencies. 
During each half-cycle, polarization on the molecule 
builds up to a peak and decays to zero. This produces a 
pulsating strain of twice the frequency of the electromag- 
netic wave. The pulsating strain will place the molecule 
in forced vibrations except at a resonant frequency, fr, 
where the molecular energy absorption will rise sharply. 
The frequency fr will depend upon whether the entire 
molecule vibrates as a two-body oscillator with each half 
of the mass of the molecule acting as one body or whether 
a prosthetic group and the remaining portion act as the 
two bodies. Because of the low strength of the hydrogen 
cross bonds and the small number required to be broken 
for denaturation, it is probably that at polarization reson- 
ance sufficient energy will be absorbed to ‘denature the 
protein for the first type of vibration. Denaturation will 
result for the second type of vibration if the prosthetic 
group is broken off. 

Now let us assume that each protein molecule of our 
laboratory experience is characterized by a unique energy 
state. For example, let us suppose that egg albumin has the 
energy Ei, and lactoglobulin, the energy E2 and that each 
energy state will resonate with one particular frequency. 
Thus, there will be a frequency fr: which will denature egg 
albumin and another frequency fre which will denature 
lactoglobulin. It then follows that each of the proteins, 
enzymes, viruses, genes, and other protein-like molecules 
can be denatured by frequencies which are related to their 
energy states. This would lead to a systematic identifica- 
tion of the molecules by their resonant frequencies. 

In this connection it seems that we should not conceive 
of a large number of lactoglobulin molecules, for example, 
as having identical energy states. It is perhaps more prob- 
ably that these molecules have a distribution of energy 
states, say a normal distribution with the mean state 
having the energy E2. Thus, the question arises as to the 
sharpness of the effect of the resonant frequency, i.e., how, 
broad is the band of energy states about the mean which 
will be denatured by fre. 

We do not know how to determine this band width at 
present but the concept is of value in setting up biological 
experiments. Let us assume two proteins, for example, 
egg albumin and lactoglobulin with mean energy states 
E, and E», respectively, with distributions of states such 
that the denaturable band of states do not overlap. If we 
irradiate the two proteins in a common solution with a 
frequency fri, which will resonate with the band of states 
about £;, but does not resonate with the band about Eo, 
we can visualize the denaturation of egg albumin without 
the denaturation of the lactoglobulin. There follows the 
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implications that a virus could be denatured without de- 
naturing the proteins of the host, or that the energy state 
of one protein, enzyme, or gene could be changed without 
changing the states of neighboring molecules having dif- 
ferent energies. 





The Recognition of Phase Transitions in 
Adsorbed Films on Solids 


D. H. BANGHAM 


British Coal Utilization Research Association Laboratories, 
Grosvenor Gardens, London 


April 19, 1946 


N a recent note Jura, Loeser, Basford, and Harkins! 
have given evidence of a first-order phase transition 
taking place in a film adsorbed at a solid/vapor interface 
(n-heptane on silver). This is not, however, the first case 
of its kind to come to light, earlier instances having been 
reported from the University of Egypt where a study was 


made of the swelling effects which accompany adsorption’ 


on charcoal; water (23.5°C)? and methanol (0°C)® are 
cases in point. 

The purpose of this note is to draw attention to the ad- 
vantages which accrue in studies of the molecular struc- 
ture of films on solids when the adsorbent is used in the 
form of a block which responds elastically by dimensional 
changes to the pressure F in the adsorbed films. The form 
of relation between F and the quantity adsorbed per g 
(s) can then be determined in two ways, viz: By direct 
measurements of the dimensional changes which accom- 
pany adsorption (Method I); and by use of Gibb’s equa- 


tion in its integrated form F=(RT/M2){ sd loge p 


(Method II); here p is the pressure of vapor, M its molecu- 
lar weight, 2 the area per g of adsorbent. In Method | it is 
necessary first to establish that the adsorbent block does 
in fact respond elastically to changes of F. This done, 
however, the method gives far more reliable information 
than Method II as to phase relationships. 

The statement of Brunauer* that in the experiments 
referred to my colleagues and-I “unfortunately did not 
measure the adsorption isotherms of the vapours they used 
in the expansion experiments” is contrary to fact. We did, 
however, consider that Method I gave more reliable re- 
sults than II because: (1) Unless a region of obedience to 
Henry's law can be established experimentally, the inte- 
gration of the isotherms in the region of low pressures is 
not merely difficult but hazardous, especially if transitions 
occur. (2) We were unable to satisfy ourselves that Gibb’s 
equation was obeyed over regions of phase transition, even 
when “time” effects and hysteresis were absent.” ° 

By Method I it was established that: (1) films adsorbed 
from vapors on the charcoal blocks used were mobile, 
and obeyed similar equations of state to films at liquid/ 
gas interfaces ;5 (2) definite phase transitions occurred ;?:* 
(3) metastable films were of frequent occurrence, transi- 
tion to a stable state manifesting itself by a contraction 
signifying increase of interfacial energy at constant ad- 
sorption ;** (4) while some adsorbed phases could be 
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classed as “‘gaseous’™” and some “liquid,” in other cases ° 


the behavior was intermediate;?*.* (5) that a two-di- 
mensional condensed phase does not necessarily confer 
complete wettability nor render capillary condensation in 
cavities inevitable; (6) that film thickening occurs with 
some adsorbates at quite small relative pressures.* 7 
These conclusions—perplexing though they appeared 
when first published—have been generally confirmed by 
later work, notably by Jura, Harkins, and collaborators in 
the course of their elaborate investigation. The general 
similarity of phase relationships at solid/gas and liquid /gas 


_ interfaces has been stressed by Gregg,* who also emphasizes 


some of the difficulties attaching to Method II. 


* In the papers quoted these were referred to as “equilibrium” or 
“annealed” films to avoid implications as to their classification. Their 
properties were, however, fully described (see D. H. Bangham, Proc. 
Roy. Soc. A 147, 175 (1934)). 

1 Jura, Loeser, Basford, and Harkins, J. Chem. Phys., 13, 535 (1945). 
ware Fakhoury, and Mohamed, Proc. Roy. Soc. A 138, 162 
one Fakhoury, and Mohamed, Proc. Roy. Soc. A 147, 152 

4). 

4S. Brunauer, Physical Adsorption (Princeton University Press, 
New Jersey, 1943), p. 433. 

5 Bangham and Fakhoury, J. Chem. Soc. 1324 (1931). 

6 Bangham and Razouk, Trans. Faraday Soc. 33, 1463 (1937). 

7 Bangham, Proc. Roy. Soc. A 147, 175 (1934). 

8S. J. Gregg, J. Chem. Soc. 696 (1942). 





Aggregational States in Adsorbed Films on 
Incompletely Wettable Solid Surfaces 


D. H. BANGHAM 


British Coal Utilization Research Association Laboratories, 
Grosvenor Gardens, London 


April 19, 1946 


F a solid surface fails to provide nuclei of condensation 
for a saturated vapor (i.e., supersaturation can occur), 
the adsorbed film must necessarily be in a state of aggrega- 
tion different from the bulk liquid. If the bulk liquid is 
brought into contact with the solid there will be a finite 
contact angle. 

It is not generally recognized that such saturated films 
adsorbed at incompletely wettable surfaces are usually of 
more than unimolecular thickness. Charcoal and meth- 
anol? charcoal and water,'* graphite and water,‘ mica 
and benzene® are cases in point. De Boer has recorded 
other instances of essentially the same phenomena but in 
which the bulk adsorbate was solid (iodine on calcium 
fluoride ;* caesium on calcium fluoride’). Earlier instances 
in which the adsorbent was liquid (water on mercury’*) 
have been reported by Frumkin and co-workers. 

It is remarkable that, where full data are available, the 
films at saturation frequently appear to be limited to 
about three molecules in thickness. Even where the iso- 
therms are sigmoid and rise steeply near saturation (this 
being the common type for non-porous solids) they con- 
verge to an asymptote beyond saturation pressure and the 
saturation adsorption appears limited to about four molec- 
ular layers. The data of Palmer® for various organic vapors 
adsorbed on vitreous silica (the surface area of which was 
measured independently) provide numerous examples. 

For explanation of such steeply rising isotherms (convex 
to pressure axis) we must clearly think in terms of a pal- 
terned assembly of adsorbed molecules, the interaction be- 
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tween which and the solid exceeds the sum of its interac- 
tions with an equal number of isolated molecules. The 
contact angle measures the incongruence of this pattern 
with that in the normal liquid and at the liquid/vapor 
interface. Substances completely wettable by water, there- 
fore, are those upon which the pattern of adsorbed water 
molecules resembles that of the molecules of ordinary 
water or can fade imperceptibly into it. Incompletely 
wettable surfaces are not necessarily those which do not 
adsorb water, but those on which the pattern is incompatible 
with that of normal water. It is highly improbable that this 
pattern can be identified with that of ordinary ice, since 
this already exists in water at room temperature, and there 
is no evidence that contact angles decrease markedly 
below 4°C. Hysteresis of the contact angle (where real, 
and not due to contamination) would imply that the ad- 
sorbed film can exist in either of two configurations of 
equal stability. 

By the use of a simple technique!’ it has been found 
possible to build up adsorbed layers of water and of or- 
ganic vapors on mica to visible thickness (from super- 
saturated vapors) and demonstrate their incongruity with 
the normal liquids. The differences in behavior between 
water and such substances as benzene and carbon tetra- 
chloride were too slight to suggest that the more general 
aspects of the phenomena depend on the hydrogen-bond 
structure peculiar to water. 

Though the heat of adsorption in the outer layers is 
not necessarily equal to the normal latent heat (with methyl 
alcohol on charcoal it exceeds it) there are other cases 
where the difference is small, indicating that the ‘“‘ab- 
normality” lies in the entropy rather than in the energy 
of the film. 

These observations may be helpful in connection with 
the. theory of solutions since they offer suggestions as to 
the nature of the configurational changes which take place 
in the near neighborhood of solute molecules. 


“ae H. Bangham and R. I. Razouk, Trans. Faraday Soc. 33, 1463 
(1937). 

?D. H. Bangham, Proc. Roy. Soc. A 147, 176 (1934). 

3D. H. Bangham, N. Fakhoury, and A. F. Mohamed, Proc. Roy. 
Soc. A 138, 162 (1932). 

‘F. A. P. Maggs, unpublished. 

’D. H. Bangham and S. ea bag Roy. Soc. A 165, 552 (1938). 

‘J. H. De Boer, Physica, 8, 145 (1928). 

7 J. H. De Boer, private rihanna es 

* A. Frumkin, A. Gorodetzskaja, B. Kabanow, and N. Nekrassow, 
Physik Zeits. Sowjetunion 1, 255 (1932). 

°W.G. Palmer, Proc. Roy. Soc. A 160, 254 (1937). 
ommn H. Bangham and Z. Saweris, Trans. Faraday Soc. 34, 554 

38). 





Sedimentation Equilibria of Polydisperse 
Non-Ideal Solutes 
MICHAEL WALES, MARGARET BENDER, and J. W. WILLIAMS 
Department of Chemistry, University of Wisconsin, Madison, Wisconsin 
AND 
ROSWELL H. Ewart 
General Laboratories, United States Rubber Company, 


Passaic, New Jersey 
_ interpretation of equilibrium ultracentrifuge data 
for solutions of extended long chain molecules has 
been rendered difficult in the past by pronounced devia- 
tions from Henry’s law at very low concentrations. This 
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difficulty can be overcome by taking account of the devia- 
tions from ideality, using thermodynamic data from os- 
motic pressure measurements. 

Consider a column of solution of unit cross section and 
length dx, in a centrifugal field. At equilibrium, the change 
in potential energy in moving c, grams of material through 
the distance dx to a point where the concentration is 
¢z+dc, will be balanced by the work done against the os- 
motic pressure difference (0p/0x)dx. 

Let 

p=(RT/ Mazer +b'c2 (1) 


where =osmotic pressure, in consistent units and M,, is 
the number average molecular weight of material at dis- 
tance x. 


Then 
C2(1— Vp )w®xdx = (0p/dx)dx. (2) 
Now differentiating (1) and remembering that 
~ 
mili 
Mn: = Dici/M? 
we obtain 
Liciw'x(1—V ’p)=RTZ; Fa. - +20’ oe ici (3) 
This can be considered to be the sum of 7 equations, so that 
= RT dc; , dz 
ciwx(1— Ve) =F a7 +5; (4) 


where 
2 "Dic, = Zid,’ Ci. 


It is to be noted that the sedimentation of each species is 
affected by the presence of the others. It is assumed that 
b;’ =2b’, as supported by the theory of Scott and Magat! 
and in the absence of evidence to the contrary. Proceeding 
from this point in a manner similar to that of Lansing and 
Kraemer? expressions for My: and M., are obtained, as 





follows: 
a'BZ 
a’Z 2A «x 
cn Sip <eusinmnninesdinian ES ren 5 
Moz cz(2Ax—Ba'Z) Be (1 _a’B =) ( ) 
. 2A x 
and 
in 22%! Z2X1 ie a(; Sire 
a Zix2 1—a’BZ,/2Ax (6) 





~ A(x?—x)—Ble2—c1) 


In these expressions, 








4a V0) 

. 2RT 

,_ 1 _de/dx 
“Ge 2 
_ 20gb’ 
ae. 


when c, is measured in g/100 g solution and 


mm solvent 


b’ in ( g 
100 ¢g 
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TABLE I. 











Ex- Uncor- Concen- Concentration 
peri- rected tration ratio between Corrected 
ment M: g/100g R.P.M. ends of cell Mw M: 
23 95,000 0.108 4700 4 97,000 153,000 
24 86,000 0.108 6270 10 107,000 141,000 
25 120,000 0.072 6440 10 108,000 150,000 
Osmotic pressure Sedimentation data, corrected Theoretical 
Mn Mw M: M:|Mw M:|Mw 
50,000 104,000 148,000 1.42 1.50 








The averaging formula for M, becomes 





» M,,Zdx 
_ Ja 1a’ BZ/2Ax 
M.=—— Zdx ; 


a 1—a’BZ/2Ax 


The averaging formula for My, remains unchanged. Some 
results obtained with the equilibrium ultracentrifuge for 
the system polystyrene-carbon tetrachloride are given be- 
low in Table I. 
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The constant B in c.g.s. units for carbon tetrachloride- 
polystyrene (c, in g/100 g) was 0.1705 X 10-4. The molecu- 
lar weight distribution as calculated from the assumed 
reaction mechanism gives 


M,: Mo: M,=1:2:3. 


Also it has been found mathematically possible to cal- 
culate the differential weight distribution curve of the 
polymer from one experiment in the equilibrium ultra- 
centrifuge, provided };’=2b’. The method is similar to 
that of Rinde.* 

However, the calculations are extremely tedious and they 
are of doubtful accuracy when based on only one experi- 
ment. In addition to this the curve will be skewed to a 
certain extent depending on the difference of b;’ from 2)’. 

The authors are indebted to Dr. R. A. Gregg of the 
General Laboratories, United States Rubber Company, 
for the preparation of the polystyrene used. 

1R. L. Scott and M. Magat, J. Chem. Phys. 13, 172 (1945). 

2W. D. Lansing and E. O. Kraemer, J. Am. Chem. Soc. 57, 1369 


(1935). 
3H. Rinde, Dissertation, Upsala, 1928. 
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